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1. — Introduction 

According to the present understanding of particle physics, the fundamental inter- 
actions among elementary particles arc described by the so-called Standard Model (SM) 
of the electroweak and strong interactions. It is a non abelian gauge theory based on the 
group SU{2) (g) U{1) for the description of the electromagnetic and weak forces in unified 
form (the so-called standard electroweak theory) [1] and the group SU{3) colour for the 
strong interactions (the so-called Quantum Chromo dynamics or QCD) [2]. 
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In the SM the fundamental constituents of matter are quarks and leptons (mat- 
ter fields), while the quanta mediating the interactions among the matter particles are 
the photon, the and bosons and the gluons (gauge fields). Both quarks and 
leptons interact via the electroweak force, while quarks in addition do participate also 
in strong interactions. In the electroweak sector of the SM, the fundamental fermions 
are grouped into throe different generations. Each generation contains a charged lep- 
ton (e,/Lt, t), an associated neutrino (j/g, i^^, i^x), an "up" quark {u,c,t) and a "down" 
quark {d, s,b). Fermions belonging to different generations have identical couplings but 
are distinguishable by their masses. Although the SM docs not predict the number of 
sequential generations, the present precision measurements at LEP/SLC rule out the ex- 
istence of a fourth generation with a neutrino of mass up to 45 GeV (see Sect. 4'1). All 
the matter particles assumed by the SM have been experimentally discovered. In fact, 
the long-awaited top quark has been recently observed by the CDF and DO Collabora- 
tions at the TEVATRON, and the value obtained for its mass {mt = 175.6 ± 5.5 GeV) 
is in good agreement with the mass range indirectly obtained from LEP/SLC precision 
measurements (see Sect. 4' 2.1). 

In order to realize the unification of electromagnetic and weak interactions, the 
standard electroweak theory predicts the existence of four vector bosons as carriers of the 
electroweak force: 7, W^, The electromagnetic interaction is described in terms of the 
massless photon 7. Both charged and neutral current weak interactions are predicted by 
the theory and described by the exchange of heavy and bosons. The intermediate 
vector bosons and and all the fermions, acquire mass by the mechanism of 
spontaneous symmetry breaking of the SU{2) U{1) gauge symmetry, the so-called 
Higgs mechanism [3]. Strictly speaking, no mass value is predicted by the theory; on the 
other hand, by making use of low-energy experimental data, it is possible to determine 
the vector boson masses and the masses of all the fermions with the exception of the 
top quark, for which high-energy data are necessary. As a consequence of the Higgs 
mechanism, the existence of a neutral scalar boson, the Higgs boson, is predicted by 
the theory. However, also the mass of this particle is a free parameter that has to be 
determined experimentally. At present, the Higgs boson has not been yet experimentally 
observed, and its discovery would certainly represent one of the most important tests 
of the SM. Present negative searches at LEP provide the mass limit mjj > 77 GeV at 
95% CL (see Sect. 4^2.2). 

The first experimental confirmation of the electroweak theory was the observation 
of the neutral currents in neutrino-electron scattering experiments [4]. Since then, im- 
portant experimental results, such as the discovery of the and bosons by the 
UA1/UA2 Collaborations in proton-antiproton collisions [5] and the evidence for the 7- 
Z interference in e+e" — > // processes at PEP/PETRA [6], provided further support 
to the SM. With the advent of the electron-positron accelerators LEP (Large Electron 
Positron collider at CERN) and SLC (Stanford Linear Collider at SLAC), the modern 
era of precision tests of the SM was started. At LEP in its first phase (LEPl) and 
SLC the bosons arc copiously produced in a very clean experimental environment, 
enabling the determination of the electroweak observables with impressive precision and 
offering the possibility to test the predictive power of the SM at the level of its quantum 
structure. 

LEPl and SLC colliders began their operation in the fall 1989. The energies of the 
electron (e~) and positron (e+) beams were chosen to be approximately equal to 45 GeV, 

in such a way that the available centre of mass energy was centered around the mass of 
the neutral vector boson of the weak interaction, the Z'^ boson (i?cm — 91 GeV). LEPl 
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was terminated in the fall 1995 in order to allow a second phase in the LEP experimental 
program (LEP2) at energies above the peak, while SLC will continue in data taking 
close to the resonance; still for a few years. The high-energy reactions studied at 
LEPl by the four experiments ALEPH, DELPHI, L3 and OPAL and at SLC by the SLD 
Collaboration are two-fermion production processes in e+e~ collision, i.e. 

with 

// = vv{ygV^.,v^v^i.,VTV^^ (invisible final states), 

Z"'"Z~(e"'"e~,/i"'"/i~,T"'"r~) (charged leptons final states), 

qq{uu, dd, ss, cc, bb) (hadron final states) . 

About 16 millions of Z'^ events have been detected at LEP and about 100 thousands at 
SLC, using in the latter polarized electron beams. Thanks to the very large data sample, 
a remarkable precision has been reached in the measurements of the Z-boson properties 
and observables, such as mass, total and partial decay widths, production cross sections 
and forward-backward asymmetries. For example, the mass of the Z° boson is at present 
known with a relative error of 2 x 10"^, Mz = 91186.7±2.0 MeV; the total Z° decay width 
has been measured with a relative error of 1 x 10^"^, i.e.Tz = 2494.8 ±2.5 MeV [7,8]. The 
evolution with time of the experimental uncertainties on Mz and Tz during the years 
1989-1997 is shown in Fig. 1. Three years of data taking at LEPl (1991, 1993 and 1995) 
were devoted to a precision scan around the resonance, doing measurements also off-peak 
in the range 88-95 GeV; the rest of data taking was performed exclusively on the Z° 
peak. Several technological ingredients contributed to reach the fantastic performances 
of the LEPl program summarized in Fig. 1. For instance, the method of resonant spin 
depolarization adopted in the LEP beam energy measurement allowed to reach a very 
precise calibration, of the order of 1 MeV. Quite curiously, many important unexpected 
phenomena affected LEP energy calibration during operation: tide effects and field rise 
due to the passage of the TGV on the nearby railway are probably the two most famous 
examples. The installation of precision luminometers also contributed significantly to 
the success of the LEP program, allowing a luminosity measurement at the level of 0.1% 
or better, that is unavoidable to make optimal use of the high statistics collected at 
LEP. In parallel with the technological and experimental progress, a large effort was 
undertaken by himdreds of theorists in the calculation of radiative corrections. The 
excellent performance of the machine, associated with the above theoretical effort, allowed 
to test the fine structure of the Standard Model of electroweak interactions with an 
unprecedented level of precision (precision physics at LEP). In this sense, precision tests 
at the Z^ peak can be compared with classical tests of Quantum Electrodynamics (QED), 
like g — 2 and Lamb shift experiments. Besides a highly nontrivial test of the standard 
theory at the level of quantum loops, the availability of both accurate measurements and 
calculations provided the opportunity of putting limitations on models predicting physics 
beyond the standard model. 

After the long data taking at the pole and a short run at intermediate energy 
in the 1995 fall (LEPl. 5 phase), the second phase of LEP (LEP2) started operating in 
1996. The main reason for the energy upgrade of the LEP machine was the precise 
measurement of the properties of the charged vector bosons of the weak interaction, the 

bosons, through the study of the reaction 

e+e- ^ W+W- ^ /1/2/3/4 , 
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Fig. 1. - The time evolution of the errors on Mz and Tz (from [7]). 



with the following possible four-fermion final-states 

Ivil'vi' (I = e, fi, t) (leptonic final states), 

li^iQiq2{l = e, fj,,T;q = u, d, c, s) (semi-leptonic final states), 
91923394(5 = u, d, c, s) (hadronic final states) . 

The 1996 LEP2 run was split into two parts. The first run was at the centre of mass 
energy = 161.3 GcV, i.e. 0.5 GeV above the nominal WW production threshold, while 
the second one at ^/s = 172 GeV. The current LEP2 schedule foresees running during 
the years 1997-1999 with centre of mass energies from 184 GeV in 1997 to 194 GeV in 
1998-1999. Compared to LEPl, the LEP2 physics potential will be characterized by a 
statistical error of the order of 1%, instead of 0.1% at LEPl. Nonetheless, also LEP2 
has to be considered as a machine for precision physics, since it will allow to measure 
the mass of the W boson with an envisaged precision of 40-50 MeV, that is better than 
the current best determination of Mw at hadron machines, and to study in detail the 
couplings of the Ws with the other gauge bosons 7 and Z^, thus directly probing the 
non-abelian nature of the electroweak theory. Therefore, the LEP2 experimental program 
does constitute a nice completion of the precision studies at LEPl and SLC. Needless to 
say, besides the precision measurement of the properties of the charged vector bosons, 
LEP2 could also provide important pieces of information on the Higgs boson and physics 
beyond the SM. 
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In the recent few years, several reviews on the status of precision measurements at 
LEPl/SLC and of the calculation of radiative corrections in the SM have been published 
in the literature. From an experimental point of view, the main aspects of precision tests 
at LEPl/SLC have been summarized for instance in ref. [9]. The concepts and results 
of the calculation of radiative corrections have been reviewed for instance in ref. [10]. 
Before the start of the LEPl/SLC operations, the status of theoretical calculations and 
computational tools was described in refs. [11-13]. In view of the final analysis of the 
precision data, the status of precision calculations was critically reviewed in ref. [14], 
with special emphasis on the uncertainty inherent to the theoretical predictions. The 
theoretical results and related software necessary for the experiments at LEP2 were 
summarized in ref. [15]. 

The LEP/SLC physics program is a very wide one, ranging from precision measure- 
ments of vector boson properties, to QCD studies, searches of new particles, determina- 
tion of dynamical properties of heavy flavours, and so on. The main goal of the present 
review is to give a comprehensive account of precision physics at LEP, meant as the whole 
of precision tests of the electroweak sector of the SM and its implications for new physics. 
More precisely, the theoretical apparatus needed for precision calculations at LEP/SLC 
will be described and its link with data analysis pointed out, aiming at emphasizing the 
intimate and fruitful connection between theoretical ideas and experimental results. Too 
technical details are inevitably omitted; however, particular care has been devoted to 
compiling an as complete as possible bibliography, where the interested reader can find 
additional information. The experimental data considered throughout the paper are the 
ones presented at the 1997 Summer Conferences [7,8,16]. 

The review is organized as follows. Section 2 is dedicated to small-angle Bhabha 
scattering and its relevance for the luminosity measurement. In Sect. 3 the large-angle 
processes at the resonance, relevant for the study of the Z-boson properties, are con- 
sidered, reviewing the theoretical ingredients necessary for high-precision calculations of 
the Z-boson observables. It is then shown how the theoretical tools developed are used 
to fit the experimental data, in order to determine the top-quark and Higgs-boson masses 
and the strong coupling constant , and eventually establish the existence of new physics 
beyond the SM (Sect. 4). The most important issues of electroweak physics at LEP2 are 
covered in Sect. 5, where two-fermion processes far from the resonance (Sect. 5T) and 
four-fermion processes at and above the W-boson pair production threshold (Sect. 5'2) 
are discussed. At last, the conclusions are drawn in Sect. 6. Three technical appen- 
dices are dedicated to relevant theoretical subjects common to various items of LEP 
physics, namely the imiversal photonic corrections (Appendix A), the vacuum polariza- 
tion correction (Appendix B) and the scalar integrals and dimensional regularization 
(Appendix C). 

2. — Small-Angle Bhabha Scattering and the Luminosity Measurement 

The process e~^e~ e~^e~ (Bhabha scattering) is a peculiar one. Actually, at 
a difference from all the other two-fermion production processes, that occur as e+e" 
annihilation, i.e. are s-channel processes, the amplitude for e+e~ — > e+e~ receives con- 
tributions both from s- and t-channel diagrams (see Fig. 2). In the energy region typical 
of LEPl/SLC (v^ ~ Mz, Mz being the Z-boson mass) the Bhabha scattering process 
exhibits completely different features, depending on the fact that the final state elec- 
trons arc detected at large or small scattering angles. In the first case, the amplitude is 
dominated by s-channel subprocesses, and in particular by the e+e^ Z'^ annihilation. 
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in such a way that the process is sensitive to the Z-boson properties. In the second 
one, the smaller is the scattering angle the larger is the contribution of the t-channel 
photon-exchange diagram, in such a way that for sufficiently small scattering angles 
Bhabha scattering is essentially a pure QED process, substantially independent of the 
physics. These dynamical features of Bhabha scattering allow to use it both as a tool 
for studying Z'^ physics (large-angle cross sections and forward-backward asymmetries) 
and as a tool for the high precision luminosity monitoring (small-angle cross section). 
The first case will be examined in Sect. 3, while the second one will be addressed in the 
following. Before starting the discussion on small-angle Bhabha scattering, it is worth 
noticing that a large amount of work has been dedicated to Bhabha scattering in general. 
The interested reader is referred for instance to refs. [17] and [18] for a detailed account 
and a comprehensive compilation of the relevant literature. 

2'1. Luminosity Monitoring. - The luminosity £ of a coUider is the proportionality 
constant between the event rate dN/dt and the corresponding cross section a for any 
given process, according to the relation 

(1) ^=Ca, N = a j dtC = aL. 

In the practice, an experiment measures the number of events for a given process and, 
by making use of the inverse of eq. (1), quotes the experimental cross sections as 

(2) a=jN. 

From eq. (2) it is clear that, in order to fully exploit the experimental information con- 
tained in TV, the error affecting the luminosity L must be smaller than the experimental 
error affecting N. The luminosity of a collider depends in a highly non-trivial way on ma- 
chine and beam parameters. It is of course possible to compute it given these parameters, 
but, in particular for LEP, such a determination is affected by an intrinsic uncertainty 
which is completely unsatisfactory in view of the extremely high experimental accuracy. 
This is why the luminosity is monitored adopting a different strategy, namely by identi- 
fying a process which is in principle not affected (or only slightly affected) by unknown 
physics, so that its cross section can be computed within a firmly established theory, and 
determining the luminosity via the relation 

(3) L = — - — N. 

^known 

The luminosity determined by eq. (3) is then used for the determination of all the ex- 
perimental cross sections via eq. (2) (for a review on the precision determinations of the 
accelerator luminosity in LEP experiments, the interested reader is referred to [19] and 
references therein). 

Such a process does exist, and is the Bhabha scattering e~^e~ e^e~ at small 
scattering angle. The main reasons why the small-angle Bhabha scattering (SABH) 
fulfils the requirements just described are the following: 

• the SABH process is substantially a QED process, dominated by photon exchange 
in the t-channel; this, in turn, implies that 

— its theoretical cross section is dominated by a contribution that is in principle 
calculable by means of perturbative QED at arbitrary precision; 
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— the ,Z-boson exchange contribution to its cross section, via Z-boson annihila- 
tion in the s channel, Z-boson exchange in the t channel and Z-j interferences, 
is very small; hence a detailed knowledge of the Z-boson properties, like the 
precise value of its mass and decay width, to be determined in large-angle 
processes, has a negligible influence on the luminosity monitoring; 

• the SABH cross section is large, and can be rendered much larger than the typical 
Z-boson annihilation peak cross sections provided the detection angular region 
for the SABH events is sufficiently close to the beam pipe, so that the statistical 
uncertainty of N in the r.h.s. of eq. (3) can be kept small. 




Fig. 2. - The Feynman diagrams for the tree-level Bhabha scattering cross section. 

The tree-level differential cross section corresponding to t-channel photon exchange 
is given by (see for instance ref. [17] and references therein) 

where t!) is the electron scattering angle, \/s is the total centre of mass (cm.) energy and 
a is the QED coupling constant. 

The calculation of the theoretical SABH cross section to be inserted into eq. (3) 
will involve in general the calculation of radiative corrections. These are numerically 
dominated by photonic corrections, i.e. by those corrections obtained by adding real or 
virtual photon lines to the tree-level amplitudes, which in turn depend very critically 
on the phase space integration region. Hence, in order to exploit eq. (3) for the high 
precision luminosity monitoring, the theoretical cross section must be computed by taking 
into account carefully the experimental definition of a SABH event. Skimming over the 
details of the various SABH Event Selections (ES's). which vary from experiment to 
experiment, the general features of SABH event are the following: 

• the forward and backward luminometers cover an angular region of few degrees, 
starting from, say, 1°; 

• a SABH event is detected as a signal coincidence of the forward luminometer with 

the backward luminometer; 
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• for technical reasons, asymmetric angular acceptances are adopted; for instance, 
a narrower (N) acceptance on the e+ side and a wider (W) acceptance on the e~ 
side; the results of the NW and WN ES's are averaged; 

• for technical reasons, a calorimetric measurement on the final-state electrons is 

performed; this means that experimentally a final-state electron cannot be distin- 
guished by a final-state electron accompanied by coUinear photons; 

• very mild cncrgy/acoUinearity cuts are imposed on the scattered electrons, com- 
patible with the previous items. 

From the general features of a SABH event, and taking into account that photonic 
radiation is mostly soft and/or coUinear to the charged lines, the following picture can 
be drawn: 

• events with the final-state electrons inside the luminometers, accompanied by soft 
photons or initial-state photons lost in the beam pipe are detected as good SABH 
events; 

• since initial-state radiation (ISR) causes a boost of the cm. of the reaction, and 
since the final-state electrons cannot be too acollinear, there is a natural cut-off on 
the amount of ISR allowed; 

• events with the final-state electrons inside the luminometers accompanied by col- 
linear photons are detected as good SABH events; this means that final-state ra- 
diation (FSR) must be integrated over a finite region surrounding the final-state 
electron; this in turn implies, for a sort of Kinoshita-Lee-Nauenberg (KLN) mech- 
anism [20], that for FSR the so-called "logarithmically enhanced" corrections are 
greatly suppressed (more on this in Sect. 2"3); 

• initial- or final-state photons can convert into additional fermionic pairs, which for 
most of the events are lost in the beam pipe or coUinear to the final-state electrons; 

• events accompanied by the radiation of virtual photons are degenerate with the 
elastic event; 

• since in general there is no particle identification in the luminometers, it can happen 
that, for instance, a final-state electron radiates a hard photon in such a way that 
the photon hits the luminometer, while the electron is lost; also such an event is in 
general registered as a good SABH event. 

Given the above picture, it is clear that in order to perform the phase space in- 
tegration over all the configurations that correspond to a good SABH event the most 
versatile tool, i.e. a Monte Carlo (MC) integrator and/or event generator, is mandatory. 
Nonetheless, as a matter of fact the calculations necessary to compute the theoretical 
SABH cross section arc so sophisticated that all kind of information, even analytical 
or semi-analytical, are precious. A detailed account of the status of all the theoretical 
calculations available at present will be given in Sect. 2'4. It is worth noting that the 
total error affecting the determination of the luminosity by means of eq. (3) receives 
contributions both from the experimental error of the number of events N and from the 
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theoretical error of the SABH process cross section. Actually, by performing the standard 
propagation of errors on eq. (3) one obtains 

6L 6N 6a 

where the first term in the r.h.s. is the experimental error, the last one is the relative 
error on the theoretical cross section and ® means the sum in quadrature. 

2'2. Sensitivity of Observables to Luminosity. - The LEP experiments determine the 
.Z-boson parameters by means of combined fits to the measured hadronic and Icptonic 
cross sections and leptonic forward-backward asymmetries (see Sects. 3 and 4 for more 
details). In particular, in order to facilitate the comparison and the combination of the 
results, the LEP Collaborations quote, besides other choices, the following set of nine 
parameters: 

(6) Mz, Tz, al Re,^,r, J^'fb\ 

where Mz and Vz are the Z-boson mass and total width, respectively, is the hadronic 
peak cross section for Z-boson exchange only, Ri are the ratios of the hadronic width 
to the i-th leptonic width and A^p^ are the leptonic forward-backward asymmetries at 
\fs = Mz for Z-boson exchange only. These parameters are quoted after having been 
corrected for the effect of ISR, final-state QED and QCD corrections, as well as for t- 
channel and s-t interference for e+e~ final states. By assuming lepton universality, the 
three values of i?, and combined into a single Ri and Ap^ , and the corresponding 

fit is a 5-paramctcr fit. In both cases, the so called derived parameters, such as for 
instance the leptonic width and the number of light neutrinos, are computed by expressing 
them as functions of the fundamental parameters. 



Table I. - Line shape and asymmetry parameters from 5-parameter fits to the data of the four 
LEPl experiments, made with a theoretical luminosity error of 0.16%, 0.11% and 0.06%. In the 
lower part of the Table also derived parameters are listed (from ref. [18]). 



theoretical luminosity error 


0.16% 


0.11% 


0.06% 


Mz [GeV] 


91.1884 ± 0.0022 


91.1884 ±0.0022 


91.1884 ± 0.0022 


Tz [GeV] 


2.4962 ± 0.0032 


2.4962 ± 0.0032 


2.4961 ± 0.0032 


al [nb] 


41.487 ± 0.075 


41.487 ±0.057 


41.487 ± 0.044 


Ri 


20.788 ± 0.032 


20.787 ± 0.032 


20.786 ± 0.032 




0.0173 ± 0.0012 


0.0173 ± 0.0012 


0.0173 ± 0012 


Th [GeV] 


1.7447 ± 0.0030 


1.7447 ± 0.0028 


1.7446 ± 0.0027 


r, [MeV] 


83.93 ± 0.13 


83.93 ± 0.13 


83.93 ±0.12 


a? [nb] 


1.9957 ±0.0044 


1.9958 ± 0.0038 


1.9959 ± 0.0034 


Th/Fz [%] 


69.90 ± 0.089 


69.90 ± 0.079 


69.89 ± 0.072 


Ti/Tz [%] 


3.362 ± 0.0037 


3.362 ± 0.0032 


3.362 ± 0.0028 


ri„„ [MeV] 


499.9 ± 2.4 


499.9 ±2.1 


499.9 ± 1.9 


Tinv/Tl [%] 


5.956 ± 0.030 


5.956 ± 0.024 


5.956 ± 0.020 




2.990 ± 0.015 


2.990 ±0.013 


2.990 ± 0.011 



In order to understand the effect of the luminosity error on the physical observables, 
one can assume, as a first approximation, that the luminosity error affects, among the 
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fundamental parameters, the hadronic peak cross section only, in a way that is determined 
from eq. (2) by performing the standard propagation of errors: 

(7) 5<jl = {5alU,®al^-^. 

In eq. (7) the first term in the r.h.s. is the experimental error on the hadronic peak cross 

section, while the second term is the error on the cross section as due to the luminosity 
uncertainty. The influence of the luminosity error on the derived parameters is then 
determined by their dependence on ct^. 

Actually, in the practice the LEP Collaborations combine their results, propagating 
both the systematic and statistical uncertainties, and taking into account all the corre- 
lations between the parameters by means of the full covariance matrix. Table I shows 
an example of the results of such a procedure, by assuming a relative luminosity error of 
0.16%, 0.11% and 0.06%, respectively. As can be noticed, the luminosity error affects, 
besides the hadronic peak cross section cr^, almost all the derived parameters. 

2'3. Radiative Corrections to Small-Angle Bhabha Scattering. - The radiative cor- 
rections to the SABH process are dominated by photonic corrections, i.e. by those cor- 
rections coming from graphs obtained from the tree- level ones by adding real and/or 
virtual photon lines, and by the vacuum polarization correction. This last effect is taken 
into account by simply using the running QED coupling constant a{—t), where t is the 
squared four- momentum transfer of the reaction, as shown in Appendix B. In the present 
section, some generalities on photonic corrections are reviewed. 

Once the infra-red (IR) divergence present in real and virtual corrections separately 
has been canceled by properly summing over all the degenerate states, the n-th order 
QED correction takes on the following form 

n 

(8) a(") = (^)";^ai")r^(g^), 

fe=0 

where cr(") is the n-th order contribution to the corrected cross section, C{Q'^) — Lq2 — 1 
and Lq2 = IndQ^I/rrig) is the so called collinear logarithm, being a typical scale 
involved in the process. The collinear logarithm Lq2 originates from the phase space 
integration of the emitted photon, and in particular from those configiirations in which 
the emitted photon is almost collinear to the radiating charged fermion line. In the 
SABH process, the relevant scale entering the collinear logarithm is = —t, t being 
the squared four-momentum transfer of the reaction. It is worth noticing that, due to 
the smallness of the electron mass nie, Lt is of the order of 15 for — i ~ 1 GeV^, i.e. for 
values of the four-momentum transfer typical of small-angle processes at LEP. This is 
the reason why QED radiative corrections are numerically very relevant. 

The coefficients a^"'' in eq. (8) are in turn given by the following general expression 

(9) 4"^ = EO^ 

where / is the so called IR logarithm, I ~ \n{E / E™"-'^) , E being the beam energy and 

j^max i^gjj^g gome cutoff on the maximum photon energy, respectively. As a property of 
QED corrections, the IR logarithm I can become huge when shrinking the photon phase 
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space. In general, however, such corrections depend critically on the details of the ES 
considered. 

The purpose of the present disc;ussion is to derive a sort of rule of thumb able 
to provide the order of magnitude of the various QED corrections. To this aim, an 
approximate calculation is a useful guideline. As a general warning, however, it is worth 
noting that the coefficients ak, due to the presence of the IR logarithm, are IR sensitive, 
so that they can become very large in pathological situations, i.e. when very tight cuts 
on the emitted photons are imposed. After this caveat, the corrected cross section a for 
the SABH process in the leading logarithmic (LL) approximation, obtained for instance 
by means of the Structure Function method, taking into account initial- and final-state 
radiation but neglecting convolution effects (see Appendix A for more details) can be 
written as [17] 

(10) cr ~ aoe^. 

The symbols in eq. (10) have the following meaning: ctq is the tree- level cross section, /3 
is given by 

(11) /3 = 2-£(Q2) 

TT 

and e — SiEf is the product of two, in principle different, cutoffs for ISR and FSR. 
The series expansion of eq. (10) up to third order in a reads: 

(12) a~ao (^l+ZJlne+^/J^ln^e+i/jSln^e^ +0{P^). 

It is worth noting that in the LL approximation the coefficients a„ are proportional to 
the n-th power of the IR logarithm I. 

Table II. - The canonical coefficients indicating the generic magnitude of various leading and 

subleading contributions up to tliird-ordor, for a non-calorimctric ES. The coUinoar logarithm 
Lt = L = ln(|f|/mg) is calculated for 'ffmin = 30 mrad and "dmin ~ 60 mrad and for two values 
of the cm. energy: at LEPl (y^ — Mz), where the corresponding \t\ = (s/4)i?^j„ are 1.86 
and 7.53 GeV^ and at LEP2 energy (y^ = 200 GeV), where the corresponding \t\ are 9 and 
36 GeV^, respectively (from ref. [18]). 
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In order to exploit the information contained in eq. (12), it is necessary to specify 
the main features of the ES one is considering. Let us begin with a non-calorimetric 

ES. It is not a realistic case; nonetheless, it can be considered as a very useful benching 
situation. Such an ES is characterized by the fact that in principle it is possible to 
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separate final-state fermions from the photon(s) radiated by tliem, in sucli a way that 
a maximum photon energy cutoff, irrespective of the fact that the photon comes from 
initial- or final-state particles, is meaningful. For such an ES, one can take in eq. (12) 

(13) £ ~ £^ ^ 4 ^ Sl„ 

in such a way that eq. (12) becomes 

(14) <7~(7o (^1 + 2/3 In £„e + ^ (2/3)' ln2£„, + 1(2/3)3 ln3£„,^ + 0(/3^). 

By defining now as canonical coefficients the coefficients appearing in front of the IR 
sensitive terms, the canonical coefficients of the 0{a"L") corrections, with 1 < n < 
3, can be directly read off eq. (14), and are 2/3, (2/3)^/2 and (2/3)^/3!, respectively. 
The algorithms described in Appendix A allow to take these corrections into account 
naturally. Going beyond the LL approximation, eq. (14) does not provide information 
any more. Actually, the non-leading corrections are typically process dependent, and can 
be computed only by means of a full diagrammatic calculation. On the other hand, a 
non-coUinear photon is known to produce a correction whose typical size is a/n, which 
sets the size of the non-leading 0{a) correction. Moreover, the 0{a^L) corrections come 
from configurations for which there is one collinear and one non-collinear photon, so 
that their typical size is 2/3Q;/7r. The situation is summarized in Tab. II, where also 
a numerical estimate of the canonical coefficients for a non-calorimetric measurement 
is given. It is worth noticing that both these non-leading corrections are relevant for 
obtaining a theoretical error of the order of 0.1%. 

Let us now consider the more realistic case of a calorimetric ES (see Sect. 2'1). 
For such an ES, it is not possible in principle to separate a final-state fermion from its 
accompanying radiation. This means that one is effectively almost inclusive on FSR, i.e. 
that £/ ~ 1 and £ ~ = Sc. Equation (14) is then modified as follows: 

(15) a~ao(^l + /31n£e + i/32ln2£, + i/33ln3£,^ +Oif3^). 

By comparing eq. (15) with eq. (14), one obtains the relation between the canonical 
coefficients for a non-calorimetric ES and those of a calorimetric ES. In particular, the 
0(a"i") coefficients are reduced by factor of 2, 4 and 8 for n = 1, 2, 3, respectively. One 
can not expect a priori a reduction in the coefficient of the 0{a) correction. Hence, there 
is a reduction of a factor of 2 in the coefficient of the 0{a^L) correction. The situation 
is summarized in Tab. III. 

At present, there is complete control of the 0{aL, a, a^LP', a^L^) corrections, for an 
arbitrary ES. (^) The 0{a^L) corrections are not fully under control, in the sense that are 
either analytically known for non-realistic (non-calorimetric) ES's or only approximated 
for realistic (calorimetric) ES's (see Sect. 2'4 for more details). As it will be seen in the 
following, the approximate knowledge of the O(a'L) corrections is the main source of 
theoretical error on the SABH cross section. An illustrative plot concerning photonic 
radiative corrections to the SABH process can be found in Fig. 3. The tree-level cross 
section corresponding to the ES adopted is ao = 140.02 nb. As can be seen in the first 



(-*-) Actually, at present a theoretical error of technical origin is attributed to the O(ct'L^) corrections, 
that however is negligible with respect to the one due to missing 0{a^L) corrections (see Sect. 2'5). 
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Table III. - 
Tab. II) and 


- The relation between the canonical 
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plot, the photonic corrections reduce the cross section by an amount ranging from around 
6% to aroimd 16%, depending on the Zmin cut (see the Figure for the definition of Zmin)- 
In the second plot, one can see that the non-leading 0{a) corrections reduce the cross 
section by around 1.5%, the higher-order LL corrections increase it by about 0.5%, and 
the non-leading 0(0^ L) corrections introduce a further increase of about 0.1-^0.2%. All 
these corrections are relevant for a theoretical prediction at the 0.1% level. 

There are also other radiative corrections to the SABH process, which however are 
much smaller than the photonic ones. They are the following: 

• light pairs: those corrections arise from photons converting into // pairs, and are 
dominated by e+e~ pairs; they give typically a contribution of the order of few 
10~^ with respect to the tree-level cross section (see [24] and references therein); 

• QED radiative corrections to Z-^ interference: they alter sizably the tree-level Z-j 
interference contribution, since it changes sign when crossing the resonance; they 
are under control (see [25], [26] and references therein). 

Non-QED corrections other than vacuum polarization are absolutely negligible. 

2'4. Computational Tools. In the present section the basic features of the computer 
codes available for the calculation of the SABH process cross section are briefly summa- 
rized. The aim of this discussion is to make an inventory of the theoretical approaches 
to the problem, and of their realizations in the form of FORTRAN codes, rather than to 
give an exhaustive description of the programs, which can be found in the literature. 
BHAGEN95 [27] — It is a MC integrator for both small- and large-angle Bhabha scattering. 
It is a structure function based program (see Appendix A'l) for all-orders resummation, 
including complete photonic 0{a) and leading logarithmic 0{a^L^) corrections in all the 
channels (long writeup in [18]). 

BHLUMI [28, 29] — It is a MC event generator for small-angle Bhabha scattering. It 
includes multi-photon radiation in the framework of YFS exclusive exponentiation (see 
Appendix A'3). Its matrix clement includes complete 0{a) and leading logarithmic 
0{a'^L'^) corrections. Some non-leading 0{a'^L) corrections are also taken into account. 
The program provides the full event in terms of particle flavors and four-momenta with 
an arbitrary number of additional radiative photons. It is the standard package used by 
the LEP Collaborations for the calculation of the theoretical SABH cross section. 
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Fig. 3. - The effect of photonic corrections to the SABH process at different perturbative 
order/accuracy as a function of Zmin = E+E- / E^^^^, at y/s = 92.3 GeV. are the energies 

deposited in the positron and electron clusters, respectively. The ES adopted is CAL02 (see 
ref. [18]). The tree-level cross section for this ES is (tq = 140.02 nb. The first plot shows the 
cross section including LL 0{a), exact 0(a), exact 0{a) plus higher-order LL and the same 
plus non-leading 0{a'^) corrections. In the second plot, Aa/a means the effect of non-leading 
0{a) corrections (triangles), of higher-order LL corrections (squares), and an estimate af the 
effect of non- leading 0{a^L) corrections (circles). The numerical results have been obtained by 
SABSPV [21-23]. 



LUMLOG — It is a MC event generator for the SABH process (part of BHLUMI, see [28]). 
Photonic corrections are treated at the leading logarithmic level, in the strictly coUincar 
approximation, and in inclusive way. Structure functions exponentiated up to 0{a^L^) 
are included. At the 0{a^) it includes the leading corrections ofthe kind 0{a^L^). It 
is used to improve OLDBIS (more on this later) in the sector of higher-order photonic 
corrections. 

NLLBHA [30] — It is the FORTRAN translation of the only available fully analytical non- 
leading second-order calculation. At the 0{a^) it includes all the next-to- leading cor- 
rections 0{a^L). It is also able to provide 0{a^L^) photonic corrections and light pair 
corrections, including the simultaneous emission of photon and light pair. It is a semi- 
analytical result, at present available for a bare ES. 

OLDBIS [31] — It is a classical MC event generator for the Bhabha process (the modern- 
ized version is part of BHLUMI, see [28]). It includes exact 0{a) photonic corrections. It 
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is used to improve LUMLOG in the sector of 0{a) photonic corrections. 

OLDBIS+LUMLOG — It is the "tandem" developed in parallel with BHLUMI in order to take 

into account higher-order corrections (LUMLOG) on top of the exact 0{a) result (OLDBIS). 
The matching between 0{a) and higher-order corrections is performed in additive form. 
This means that no 0{a'^L) corrections are present. 

SABSPV [21,22] It is a MC integrator, designed for small-angle Bhabha scattering. 
It is based on a proper matching of the exact 0{a) cross section for t-channel photon 
exchange [17, 32] and of the leading logarithmic results for the full Bhabha scattering 
cross section in the structure fimction approach [33]. The matching is performed both 
in additive and factorized forms, the first form being used for comparisons only. In its 
default mode (factorized cross section) it includes the bulk of the 0{a'^L) corrections. 



1 








aL 


a 








a^L 1 


a' 














aV 











Fig. 4. - The general structure and the present situation concerning QED corrections to the 
SABH process. The thick line isolates the corrections exactly known for every ES. The dashed 
line points out the O(a^L) corrections. They are at present the main source of theoretical error, 
since they are either approximately known for realistic ES's or exactly known for unrealistic 
ES's. 

The situation of the computer codes available for the SABH process can be summa- 
rized as follows: 

• with the exception of LUMLOG, all the other codes implement exact 0{a) photonic 
corrections for the SABH process; with the exception of OLDBIS, all the other codes 
implement some form of higher-order photonic corrections; 

• BHAGEN95, OLDBIS+LUMLOG and SABSPV in the additive form include higher-order 

corrections in the leading logarithmic approximation; 

• BHLUMI and SABSPV in its default form include also the bulk of the 0{a'^L) correc- 
tions, for every ES; 

• NLLBHA includes the full set of the 0{a'^L) corrections, but hmited to a bare ES. 
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The situation concerning the QED corrections at present under control for every ES 
is described in Fig. 4. 

2'5. The Total Theoretical Error. - Any theoretical prediction performed by means 
of a perturbative theory is affected by an intrinsic uncertainty because the perturbative 
scries must be truncated at some finite order (with the possible exception of those cor- 
rections that can be resummed to all the perturbative orders, such as, for instance, the 
universal photonic corrections). Hence the theoretical error of a perturbative prediction 
is dominated by the largest unknown corrections, i.e. by the lowest order corrections not 
under control. When applying the above definition to the SABH process, it turns out 
that there are several sources of theoretical error, namely the ones quoted in Tab. IV. 
Up to date, the most important one is due to the missing photonic 0{a^L) corrections, 
which at present are not fully under control as shown in Fig. 4. 

Table IV. - Summary of the total (physical+technical) theoretical uncertainty for a typical 
calorimetric detector. For LEPl, the above estimate is valid for the angular range within 
1° — 3°, and for LEP2 it covers energies up to 176 GeV, and angular ranges within 1° — 3° and 
3° -6° (fromref. [34]). 





LEPl 


LEPl 
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Type of correction/error 
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Present 


Present 


(a) Missing photonic C(a^L) 


0.15% 


0.10% 


0.20% 


(b) Missing photonic O(a^L^) 


0.008% 


0.015% 


0.03% 


(c) Vacuum polarization 


0.05% 


0.04% 


0.10% 


(d) Light pairs 


0.01% 


0.03% 


0.05% 


(e) Z-boson exchange 


0.03% 


0.015% 


0.0% 


Total 


0.16% 


0.11% 


0.25% 



Moreover, when considering the calculation of a cross section for a realistic ES, 
one has to implement the theoretical formulation of the problem in a computational 
tool, typically a FORTRAN code performing all the numerical integrations. In so doing, 
other sources of uncertainty are added, that are of technical origin, such as, for instance, 
approximations in the formulae, numerical algorithms performing the phase space inte- 
grations and so on. On the whole, also the actual implementation in itself of a given 
theoretical formulation gives rise to a finite contribution to the total theoretical error. 

From now on, the precision reached in principle by means of a given theoretical 
formulation will be referred to as the "physical precision" of the approach, as distinct 
from the precision reached in its actual implementation, which will be referred to as the 
"technical precision" . 

As a matter of fact, the LEP Collaborations use the MC BHLUMI to compute the 
SABH theoretical cross section. Hence, a key issue is determining the theoretical error 
of the prediction by BHLUMI. Since the codes described in the previous Section differ 
from one another in the treatment of higher-order next-to- leading corrections, a careful 
comparison of their predictions together with a deep understanding of their differences 
can be used to infer an estimate of the total theoretical error on the SABH cross section. 

The most extensive work in this directions has been performed in the context of 
the Workshop Physics at LEP2, held at CERN, Geneva, during 1995, and in particular 
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Fig. 5. - Geometry and acceptance of the simple (non-calorimetric) ES BAREl. This ES re- 
stricts the polar angles i?i in the forward/backward hemispheres and requires a certain minimum 
energy to be detected simultaneously in both hemispheres. Photon momentum is not constrained 



at all. The entire "fiducial" t?-rangc, i.e. the wide (W) range, is (i^Sin, 
and the narrow (N) range is (j^mim '^max); where i9^in = i^Sin + 



-- (0.024, 0.058) rad 
= '^max - ^-d and 



^ij = ('^max — '?min)/16- This ES cau be symmetric WW or NN, or asymmetric NW (see the 
description in the figure). The energy cut s' > ZminS involves the momenta of the outgoing 
[s' = {q+ +q-f) only (from ref. [18]). 



within the Working Group "Event Generators for Bhabha scattering" [18], whose main 
tasks were 

• to make an inventory of all the available MC event generators, developed by inde- 
pendent collaborations, for Bhabha processes at LEPl and LEP2, both at small 
and large scattering angles; 

• to improve the understanding of their theoretical uncertainties by means of sys- 
tematic comparisons of MC's between themselves and with non-MC approaches. 

The main emphasis was put on SABH processes, because of the pressing need to 
match the theoretical precision of the calculations with the much improved experimental 
accuracy (< 0.1%) of the luminosity measurement. In particular, the main achievement 
of the Working Group, which is the result of a combined effort by several collaborations 
addressing several theoretical and experimental issues, was the reduction of the theoret- 
ical error on the SABH cross section from 0.16% to 0.11% for typical ES's at LEPl, and 
a first estimate of the theoretical error on the SABH cross section at LEP2 [34]. 

The various components of the theoretical error on the SABH cross section are 
quoted in Tab. IV (see ref. [34]), where a summary of the past and present situation at 
LEPl together with the present estimate valid for LEP2 is given. The errors in the Table 
are understood to be attributed to the cross section for any typical (asymmetric) ES, 
for a LEPl experiment in the angular range 1° - 3°, calculated by BHLUMI 4.03 [28]. 
In the case of LEP2, the estimate extends to the angular range 3° --6°, and also to a 
possible narrower angular range (say 4° — 6°) that may be necessary due to the effect of 
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Fig. 6. - Geometry and acceptance of the calorimetric ES SICAL2. This ES restricts the polar 
angles i?i in the forward/backward hemispheres and requires a certain minimum energy to be 
detected simultaneously in both hemispheres. No restrictions on azimuthal angles tpi are there. 
The entire "fiducial" t?-rangc, (i?^;^, i?^iax) = (0.024,0.058) rad, includes the wide (W) range 
(''mini ''max) thc narrow (N) range (li'min! '''max) exactly as depicted in the figure. This 
ES can be symmetric WW or NN, or asymmetric NW. The energy cut and the i?-cuts involve 
the definition of the cluster. Each side detector consists of 16 x 32 equal plaquettes. A single 
plaquette registers the total energy of electrons and photons. The plaquette with the maximum 
energy, together with its 3 x 3 neighborhood, is called cluster. The total energy registered in 
the cluster is Ef and its angular position is {'df,ifif ), i — 1,2. More precisely the angular 
position of a cluster is the average position of the centers of all 3 x 3 plaquettes, weighted by 
their energies (the definitions of <^'s are adjusted in such a way that (fi = (p2 for back-to-back 
configuration). The plaquettes of the cluster which spill over the angular range (outside thick 
lines) are also used to determine the total energy and the average position of the cluster, as in 
the backward hemisphere (from ref. [18]). 



synchrotron radiation masks in the experiments. The entries include combined technical 
and physical precision. 

As can be seen in the Table, at the present stage the theoretical error is still domi- 
nated by the error on photonic corrections, quoted in entries (a) and (b), and in particular 
by the one due to missing O(a^L) corrections of entry (a). Since the error of entry (a) 
is by far dominant with respect to all the other ones, it is worth devoting some space to 
describe how it has been estimated, namely by adopting the following procedure: 

• only the photonic corrections to the dominant part of the SABH cross section, 
namely the one due to t-channel photon exchange, have been considered as a first 
step; 

• four families of ES's have been defined (for the details concerning the definitions of 
the ES's the reader is referred to [18]); the simplest one, BAREl (see Fig. 5), is an 
ES in which cuts are applied only to the "bare" final-state fermions; the other ones, 
CALOl, CAL02 and SICAL2, are calorimetric ES's, implementing more and more 
complex clustering algorithms; in particular, SICAL2 (see Fig. 6) is very similar 
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to a "real" experimental ES; since photonic corrections are very sensitive to the 
details of the ES, defining these four ES's allows to span in detail the photonic 
phase space; even if the ES BAREl is far from realistic, the presently available 
analytical calculation including the complete set of 0{a'^L) corrections refers to 
such an ES, and so provides a very important cross-check of the MC programs; 

• all the available codes have been run for all the ES's, varying inside any ES the 
threshold requirements for the final-state fermions/clusters; 

• a test concerning the technical precision has been performed, namely comparing 
the exact up to 0{a) cross sections provided by the various codes; agreement at 
the level of a few 10~^ relative deviation has been achieved; 

• finally, the results of the codes including the full higher-order photonic corrections 
have been compared, for all the situations explored (see Fig. 7). 
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Fig. 7. - Monte Carlo results for the symmetric Wide-Wide ES's BAREl, CALOl, CAL02 
and SICAL2, for matrix elements beyond first order. Z-boson exchange, up-down interference 
and vacuum polarization are switched off. The cm. energy is y/s = 92.3 GeV. Zmin is the cut 
condition on the final-state energies, defined as E+E-/E'^ > Zmin, E-^+ being the final-state 
energy of the bare electron and positron, respectively. The fiducial angular range is 0.024- 
0.058 rad (for more details on the ES, the reader is referred to [18]). In the plot, the 0{a^)^J^^ 
cross section from BHLUMI 4.03 is used as a reference cross section (from refs. [18,34]). 
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The result of this procedure aUowed the definition of "one-per-mill regions" , referring 
to reahstic threshold cuts, within which most of the predictions lie. Moreover, for those 
cases for which the predictions do not lie within the "onc-pcr-mill regions" , the reasons 
for the deviations involved have been carefully investigated and eventually understood. 
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Fig. 8. - Monte Carlo results for various symmetric/asymmetric versions of the CAL02 ES, 
for matrix elements beyond first order. Z-boson exchange, up-down interference and vacuum 
polarization are switched ON. The cm. energy is = 92.3 GcV. Zmin is the cut condition on 
the final-state energies, defined as E+E-/E^ > z,nin, E-^+ being the final-state energy of the 
electron and positron clusters, respectively. The fiducial angular range is 0.024-0.058 rad (for 
more details on the ES, the reader is referred to [18]). In the plot, the 0{a^)^xp^ cross section 
from BHLUMI 4.03 is used as a reference cross section (from refs. [18,34]). 



An analogous procedure has been followed after the inclusion of all the relevant 
radiative corrections (vacuum polarization, Z-boson exchange contributions and so on), 
and extending the comparisons also to asymmetric ES's, leading to the results shown as 
an example in Fig. 8. A similar analysis has also been performed for the first time in 
situations which are typical at the LEP2 experiments. The conclusion drawn at the end of 
all these comparisons is that now the theoretical uncertainty due to uncontrolled 0{a'^L) 
corrections is reduced from 0.15% to 0.10% for the LEPl situation, and estimated to be 
0.20% at LEP2. As far as entry (b), the "missing photonic 0{a'^L^y' uncertainty, is 
concerned, new estimates of the effect have resulted in a more conservative theoretical 
error, namely 0.015% to be compared with the old estimate of 0.008%. 

As far as all the other entries in Tab. IV are concerned, namely entries (c), the 
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"vacuum polarization" uncertainty, (d), the "light pairs" uncertainty, and (e), the "Z- 
boson exchange" uncertainty, two of them, (c) and (e), are reduced with respect to the 
prc'vioiis situation thanks to several new fits of the hadronic contribution to the vacuum 
polarization, and some additional original work on the Z-boson exchange contribution 
done during the workshop. New estimates, both MC and analytical, of the light pairs 
contribution, (d), featuring more complete calculations done during the workshop, have 
resulted in a more conservative estimate of the pairs effect uncertainty of 0.03%. 

In conclusion, the total (physical + technical) theoretical uncertainty on the SABH 
cross section for a typical calorimotric detector is at present 0.11% at LEPl and 0.25% 
at LEP2. While the theoretical precision reached in the LEP2 case is sufficient, further 
improvements in the LEPl case are desirable. 

2'6. Recent Developments and Perspectives. - After the completion of the Working 
Group "Event Generators for Bhabha scattering" [18], some additional work, relevant for 
a further reduction of the theoretical uncertainty on the SABH cross section, has been 
done by the BHLUMI [35-39] and SABSPV [23] groups. 




Fig. 9. - Examples of virtual corrections to single bremsstrahlung in the SABH process (see 
ref. [35,37]). Only graphs concerning the real emission by the electron line are shown. 

Concerning the BHLUMI group, the exact virtual one-loop corrections to the hard 
bremsstrahlung process in the SABH scattering (see Fig. 9) have been computed [35, 
37]. These results arc needed to complete the exact treatment of the 0{a'^L) photonic 
corrections, since the contributions from double bremsstrahlung [40] and the two-loop 
electron form factor [41] are known. None of these corrections is at present implemented 
in BHLUMI [35]. It has to be noticed that this is the first exact, completely differential, 
result for the virtual one-loop corrections to the hard bremsstrahlung process in the 
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SABH scattering. Thanks to this result, the authors are able to estimate the size of the 
missing 0{q^)L part in BHLUMI, for a realistic NW ALEPH SICAL luminometer, finding 
it below 2 X 10^"' of the Born cross section in the massless limit. Work is in progress in 
order to generalize the result to the massive case, and in order to compare its effect with 
the semi-analytical calculation of the NLLBHA group [30] . 

Concerning the SABSPV group, the theoretical formulation has been refined in order 
to eliminate a phase-space approximation in the 0{a^L) sector and the result has been 
analytically checked against the available complete 0(0^ L) results already present in the 
literature for an academic ES (namely a BARE ES) in the soft-photon approximation, 
both for the annihilation and scattering channels [23]. For the annihilation channel, the 
results have been compared in particular with the ones shown in ref. [42], finding that 

• the 0{a) perturbative result is exactly recovered, by construction; 

• all the IR-singular terms, namely the ones containing In^ £ and Ine, where e = 
AE/E, are exactly recovered at the level of O(a^i^), 0{a'^Ls) and 0{a'^), where 
Ls = ln(s/m2); 

• the difference between the two results starts at the level of {a/n)'^Ls times a con- 
stant; 



in particular, such a difference reads 



(16) 



5a 



3C(3)-^C(2) + ^ 



where 5a is the difference between the cross section of [42] and the cross section of 
ref. [23]. The difference numerically amounts to a relative deviation of about 1.7 x 10"''. 
The residual difference is at 0{a^) times a constant and is numerically irrelevant. For 
the scattering channel, the results have been compared with the ones of ref. [30]. The 
results of the comparison are the same as in the annihilation case up to the 0{a^Lt) 
corrections, namely the difference appears at the level of {af-K j^Lt times a constant and 
reads 



(17) 



5a 



3C(3) 



where 5a is the difference between the cross section of [30] in soft approximation and 
the cross section of ref. [23]. This difference numerically amounts to a relative deviation 

of about 2.2 x 10""*, since the overall factor of two is compensated by the fact that 
Lt ~ 2/3Ls. In this case, also an additional difference appears, namely at the level of 
the IR-sensitive truly 0{a'^) terms, which reads 



(18) 



5_a_ 



= -Q'[41n^. + 81n.], 



and is numerically irrelevant for a realistic situation. Therefore, the Oia^L) corrections 
taken into account by the method of ref. [23] represent the bulk of the complete set. By 
taking now into account that the canonical coefficient of the 0{a^L) corrections for a 
calorimetric (realistic) ES is one half of the corresponding one for the BARE ES (see 
Tab. Ill), and considering a safety factor of three, the authors of [23] estimate that the 
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Fig. 10. - Comparison of Monte Carlo's. The relative differences between the codes involved 
in the comparison and the cross section by BHLUMI 4.03 taken as a reference cross section are 
shown as functions of the cut Zmin = E+E- / E^^^^. are the energies deposited in the 

positron and electron clusters, respectively. The details of the clustering algorithms (BAREl 
and CAL02) are given in [18]. The cm. energy is ^/s = 92.3 GeV. For more details see ref. [23]. 



overall theoretical error of the approach, as far as QED corrections are concerned, is 
Sa/a ~ 3 X lO"". 

The situation is described in Fig. 10, where two ES's have been considered, namely 
a bare and a calorimetric ES. As can be seen, BHAGEN95 and OLDBIS+LUMLOG, which 
do not take into account 0{a^L) corrections, differ from BHLUMI, which is taken as the 
reference cross section, by about 0.1%. NLLBHA shows differences which are contained 
within 0.1%; since the NLLBHA formulae coincide with the SABSPV ones up to the 0{a^L) 
corrections in the soft-photon region, such differences in the soft-photon region arc pre- 
sumably due to lack of exponentiation in NLLBHA. BHLUMI and SABSPV differ at most 
by 0.025%, consistently with the independent estimates of accuracies provided by their 
authors, respectively. 

The results shown pave the way to a definite reduction of the theoretical luminosity 
error at the level of 0.05%. This is anyway not far from the smallest possible theoretical 
error, since the vacuum polarization uncertainty is 0.04% in its own (see Tab. IV). 

3. - Z° Physics 

The most important processes at the Z'^ resonance arc represented by two-fermion 
production. In the case in which there are no electrons in the final state, the tree-level 
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Feynman diagrams for the s-channel annihilation e+e" 7, Z' — » // (/ ^ e) are de- 
picted in Fig. 11. In the presence of electrons in the final state, there are additional 

^-channel amplitudes, which at large scattering angles arc essentially backgrounds to the 
dominant s-channel annihilation (see Fig. 2 in Sect. 2). The high degree of precision (at 




Fig. 11. - The tree-level Feynman diagrams for the process e^e ^, Z ff (/ 7^ e). 

the level of 0.1% or better) reached in the mcasmcmcnt of the clcctroweak parameters at 
the Z° pole requires the inclusion in the theoretical predictions of radiative corrections 
beyond the Born approximation, in order to perform a meaningful comparison between 
theory and experiment. Radiative corrections to e+e" ^ 7, Z ^ // processes arc due 
to electroweak and strong forces as described by the Standard Model (SM) of the fun- 
damental interactions, based on the gauge group SU{2) (g) {/(I) SU{3). At a difference 
from the small-angle Bhabha process, where the calculation of QED radiative corrections 
only is enough to reach the aimed 0.1% theoretical precision, the two-fermion produc- 
tion in e+e~ neutral current processes receives important contributions also from the 
electroweak and strong sector of the theory. Therefore, the fine structure of the SM 
can be probed by precision physics at LEP/SLC in a highly non trivial way at the level 
of quantum loops. The issue of radiative corrections has been studied at length in the 
past. For instance, the situation before the starting of LEP operations is reviewed in 
refs. [12], [13] and [43]. More recently, in view of the final analysis of LEPl data, the 
theoretical studies have been updated in ref. [14]. 

In the present section, the basic aspects concerning the theoretical treatment of 
radiative corrections in precision calculations for the Z° resonance are reviewed, through 
the analysis of the main ingredients of pure weak, QCD and QED corrections to two- 
fermion production amplitudes. It is also discussed how the theoretical formulae are 
realized in the form of computational tools used for data analysis. The uncertainties 
associated to the theoretical predictions and their impact on the precision calculation of 
the Z-boson observables are finally analyzed. 

3'1. Realistic Observables and Parameters. - The data published by the LEP 
Collaborations and used to extract information on the SM and its possible extensions 
refer to two main classes of observables: 

• realistic observables; 

• the Z° parameters or pseudo-observables. 
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The cross sections a{s) and the forward-backward asymmetries Apb{s) of the two- 
fermion production processes 

(19) e+e-^(7,^)^//(n7), 

quoted as functions of the centre of mass (cm.) energy and including real and vir- 
tual photonic corrections {i.e. "dressed by QED") are collectively referred to as realistic 
observables. They are defined as 

(20) a = aF + (TB, 

(21) = 

where (Jf,b are the cross sections in the forward and backward hemisphere, respectively, 
i.e. 

(22) ctf = 2it I rfcosi?-^, aB = 2TT rfcosi?-^. 

In the tree-level approximation, for unpolarized initial- and final-state fermions, 
the diflFerential cross section dao/dQ, as obtained from the calculation of the Feynman 
diagrams shown in Fig. 11, is given by 



z 



^ ' dfl dfl dn dn' 

where the 7, jZ and Z contributions read 

da^ a^Q)N, 

% = -'^^''^f^Mxis)) m (1 + cos^ ^) + 29l9i cos^] , 

^ ^W^^'^"^ m? + {9lf){{9if + {9if) (l + cos^^) 
(24) +^glglglgi cos^]. 

In cq. (24) the symbols arc defined as follows. {I is the fcrmion scattering angle. is 
the muon decay constant, a is the QED coupling constant and Mz is the Z-boson mass. 
Q/ is the final-state fermion charge, in units of the positron charge. Nc is the colour 
factor, A^c = 1, 3 for leptons and quarks respectively, g* and are the tree-level vector 
and axial-vector couplings of the iih fcrmion to the Z^ boson, given by 



(25) gl = Il-2QiSiu'dw, 5^ = 



3' 



If being the third wcak-isospin component of the ith fcrmion and sin^ the squared 
sine of the weak mixing angle. At last, x(s) is the resonating factor 

(26) X{s) 



(s - M|) + iVzMz ' 
Tz being the total Z-boson width. 
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From the explicit expression of the differential cross section given in eqs. (23) and 
(24) it is straightforward to derive the expressions for the total cross section and forward- 
backward asymmetry. Concerning the total cross section, at ^/s = Mz the ^-boson 
contribution can be written in terms of total and partial Z-boson widths (see eq. (31) 
below), while the term is exactly vanishing and the 7 contribution gives an effect 
contained below 1%. As far as the forward-backward asymmetry is concerned, still at 
•^i = Mz, it takes the form shown in eq. (32) below. 

The data presented by the LEP Collaborations for the realistic observables concern 
typically two different experimental configurations: 

• extrapolated set-up: a cut on the invariant mass of the final-state fermions or on 
the invariant mass of the event after initial-state radiation (ISR) alone is imposed; 
the data corresponding to this inclusive situation are also said perfect data; 

• "realistic" set-up: simple kinematical cuts, such as cuts on the energies or the 
invariant mass of the final-state products, angular acceptance and acollinearity 

angle of the outgoing fermions, are imposed. 

Owing to the critical dependence of the QED corrections on the applied cuts, the 
comparison of the theory with the data for both the above set-up necessarily requires the 
availability of formulations (and relative computational tools) as complete as possible 
in the treatment of QED effects. Moreover, when aiming to fit "realistic" observables, 
a special effort must be devoted to the development of compact analytical and semi- 
analytical formulae for QED corrections. 

The parameters, or pseudo- observables, are extracted from the measured cross 
sections and asymmetries after some deconvolution or unfolding procedure. Basically, 
these quantities are determined by the LEP experiments by means of combined fits to 
the hadronic and leptonic cross sections and leptonic asymmetries, after corrections for 
the effect of ISR. Besides ISR, other radiative corrections or specific uninteresting effects 
can be depurated from the measurements in the unfolding procedure. For example, 
the deconvoluted forward-backward asymmetry includes the Z-boson exchange only (i.e. 
after subtraction of the 7 and 7Z contributions) and also final-state QED and eventually 
QCD corrections are subtracted from the experimental data. For the quantities relative to 
the e+e" final state, the t and s-t interference contributions are subtracted as well. More 
details about the fitting procedure adopted by the LEP Collaborations for the extraction 
of the Z^ parameters are given in Sect. 4. Therefore, the Z'^ parameters are secondary 
quantities or, in this sense, pseudo-observables. They are anyway of utmost importance 
for the precision tests of the SM, since they depend on the details of the internal structure 
of the electrowcak theory. The Z° parameters considered in the literature are given in 
Tab. V. The effective sine is defined as 



where Q f is the electric charge of the; fermion / in units of the positron charge, gy and 17^ 
are the effective neutral current vec;tor and axial-vector couplings of the to a fermion 
pair //, respectively. By definition, the total and partial widths of the Z'^ include final- 
state QED and QCD radiation. The invisible Z-boson width Tinv, the ratios R and the 
hadronic peak cross section are defined as 



(27) 



4|Q/|sin2< 



eff 




(28) 



Fz-Fe 
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Table V. - The most relevant parameters. 



Observable 


Symbol 


hadronic peak cross-section 




partial leptonic and hadronic widths 


Ti {I = e,ix,T), Fc, Tt, 


total width 


Tz 


hadronic width 




invisible width 




ratios 


Ri, Rb, Rc 


forward- backward asymmetries 


aI Ab AC 
-^FBi -^FB-! -^FB 


polarization asymmetries 


pr^ pb 


left-right asymmetry (SLC) 




effective sine 


sin^ sinHl ff 



(29) Ri = ^, 

i I 

(30) Rb,c = 

i h 

(31) ^^"'^ = ^^'^^fri- 

By definition, (7^^^ includes only the Z-boson exchange. Forward-backward, left-right 
and polarization asymmetries, unlike the widths, arc depurated, as said above, of the 
effects of QED and QCD corrections and, like cr^ad' refer to pure Z-boson exchange. 
This allows to express them as simple combinations of the effective Z° couplings as 
follows: 

(32) 4b - l^eAf, 

(33) -^LR = 

(34) Pf = -Af, 

where 



(35) 



Analogously, the decay width of the boson into a // pair is given by the following 

expression: 

(36) Tf = 4N,To[{9^rRi + iaifRi], 

where Fq is given by Fq = G^M^^ /24A/27r, and Ry and i?;^ are factors taking into account 
QED and QCD final-state radiation (FSR) and mass effects. 

3'2. Electroweak Corrections. - The high precision reached in the measurements at 
LEP and SLC allows to test the SM at the level of its radiative corrections. This can 
be done by virtue of the renormalizability of the theory [44]. The one-loop radiative 
corrections to e~^e~ — > // can be divided into two classes, separately gauge invariant: 
the pure QED corrections and the electroweak ones. The former are obtained by adding 
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a virtual or real photon line to the lowest order Feynman diagrams, while the latter 
consist of all the remaining diagrams at one-loop order. As will be seen in the following 
the two classes of corrections arc very diffcircnt in many aspects and therefore different 
methods of calculation are required in order to reach the necessary theoretical accuracy. 
Contrary to the pure QED corrections, the electroweak corrections amount numerically 
to a few percent so that a one-loop approximation is already almost satisfactory. On the 
other hand, they depend on the fundamental parameters of the theory, so that particular 
attention must be devoted to the inclusion of potentially large higher-order effects in 
order to fully exploit the experimental accuracy to get information on the still unknown 
parameters. In the present Section the fundamental ingredients of electroweak corrections 
are illustrated and discussed. 

3'2.1. One-loop Feynman diagrams. The calculation of the electroweak corrections 
involve the evaluation of diagrams with closed loops, where an integration over the loop 
momentum is present which diverges for large momenta (ultraviolet divergences). For 
this reason a regularization procedure is needed in order to deal with finite and math- 
ematically well defined integrals, before implementing a renormalization program. The 
most used method in the case of gauge theories is the dimensional regularization [45 47] 
because it allows to maintain Lorentz and gauge invariance at any step of the calculation. 
Some details about this procedure are given in Appendix C. 

The electroweak one-loop corrections to e~^e~ // can be classified in vector boson 
self-energies, fermion self- energies, vertex and box corrections. The fermion self-energies 
are commonly accounted as part of the vertex corrections. The vector boson self-energies, 
also referred to as oblique corrections, include the following transitions: 7-7, 7-Z, Z-'y, 



Fig. 12. - One-loop fermionic contributions to vector-boson self-energies. 

All the above classes of corrections can be written in terms of general two-, three- 
and four-point functions introduced in ref . [48] , which in turn can be reduced to certain 
combinations of the basic scalar one-, two-, three- and four-point integrals introduced 
in ref. [49]. Some simple examples of the procedure are given in Appendix C. The 
general expressions of the scalar form factors for arbitrary momentum and masses are 
not simple functions of their arguments and have to be computed numerically by means 
of the relations given in refs. [48-50]. 

Within the set of the vector boson self-energies two different kinds of corrections can 
be distinguished for their different numerical relevance: the fermion-loop contribution 
which is the leading one (see Fig. 12) and the bosonic contribution (see Fig. 13). Being 
gauge invariant by themselves, the fermion-loop corrections can be resummed with some 
procedure in order to take into account of reducible higher-order effects. The complete 
expressions for the vector boson self-energies in terms of two-point functions can be found 



Z-Z and W-W. 




f 

w / \ w 
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Fig. 13. - One-loop bosonic contributions to vector-boson self-energies. 



in refs. [48,50,51]. Here only the fermion contributions are recalled because they deserve 
some important comments. The 7-7 transition is given by [50] 

(37) S^ip") = ^p2^Q2^c{8S2i(p',m,m) -4iJo(p^m,m)}, 

where B are two-point functions and the sum is extended to all charged fermions. This 
expression is the same as obtained in QED, since the SM coupling between photon and 
fermions recovers the QED case. The quantity 

(38) Ip-'ip') = 

is known as the photon vacuum polarization. It is interesting to work out the two-point 
functions for the following asymptotic cases, for a fermion of given flavour: 
• light fermions (jp^j » m^) 
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(39) n-(/) = (-A + inm^ - ^ + In g - i.) , 

• heavy fermions << m^) 

(40) n-(,^) = ^Q5.(-A + lnm^-^). 

Apart from the terms A — Inm^, which are related to the ultraviolet divergences and are 
removed by the renormalization procedure, it is worth noticing the term proportional 
to ln{\p^\/m'^) in the light fermions limit. At energies of the order of 100 GeV, as is 
the case at LEP, this terms is the origin of large corrections coming from the light 
fermion spectrum of the theory. While the leptonic contribution can be unambiguously 
calculated, the hadronic one is affected by a large uncertainty, because the light-quark 
masses can not be unambiguously defined. For this reason the hadronic contribution 
to the photon vacuum polarization is calculated by means of a dispersion relation as 
described in Appendix B. As far as the top-quark contribution to the vacuum polarization 
is concerned, taking into account the top-quark mass value of about 170 -r 180 GeV, 
the heavy fermion limit of 11''"'' can be taken, resulting in a small effect. This happens 
because the U{l)e.m. group of the SM is not spontaneously broken so that the decoupling 
theorem [52] can be applied. 

The combination of two-point form factors present in the fermionic contribution to 
the 7-Z transition is the same as the 77 case apart from the couplings. It is at this 
point worth noticing that, contrary to the case of the photon self-energy, the bosonic 
contribution to the mixed j-Z transition is not proportional to so that it is different 
from zero for ~ 0. 

The massive vector boson self-energies have the following expressions: 

5+(p2) = -^N^Y, [ 2 P^{B2l(p^md,m„) + Bi(p2,md,m„)} 

d 

(41) + {ml - ml)Bi{p'^ ,md,rriu) - mlBoip'^ ,md,rnu) 

where g is the SU{2) coupling constant and the sum is extended over the fermionic 
doublets; 

So{p^) = ^^cJ2 [P^ 4^2 {8^2l(p^ m, m) - 4Bo{p\m, m)} 

2 

(42) - -^Bo(p^,m,m) 

where the sum is extended to all fermions. As before it is interesting to consider the two 
limits of light and heavy fermions: 

• light fermions 

(43) 5+(/) = ^^iV,(-A + In l/l - i^), 



(44) 5o(/) = ^Ml!^±^p^(-A + lnb^| -i.); 
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• heavy fermions 

(45) S+{p') 

(46) Soip') 



127r2 
127r2 



N, 



id? ■ 



Inml) 

id)' 



Inm; 



4r2 



(A-lnm^), 



where the contribution of a single doublet in the charged vector boson case has been 
written. Some comments are in order here. The ultraviolet divergences are associated 
only with the real parts of the self-energies, as the imaginary parts are directly related 
to the widths of the gauge bosons. Furthermore it appears in the heavy fermions limit a 
term proportional to m^, that becomes numerically relevant in the case of the top quark. 




w,z 





Fig. 14. - One-loop vertex corrections. 



As far as the vertex corrections are concerned (see Fig. 14), neglecting the external 
fermion masses only two combinations of three-point functions C{p^ ,mf,My,mf) are 
present in the calculations [48,51]: 

C/j(p2) = - 2C•24(]3^m/,My,m/) - [Cii{p^ ,m},Mv,mf) 

(47) + C23(p',my, Af„mj)]p2 + i, 

when only one gauge boson is present in the loop, and 

C^^ip") = 6 C24(/,My,m/,My)+ [Co(p',My,m/,My) 

(48) + Cii(/,My,m/,My) + C23(/,My,m/,My)]p2_ 1, 
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when two gauge bosons circulate in the loop. In the case of the bb final state, due to 
the effect of internal top-qu&ik lines, five different combinations of three-point functions 
are needed, whose explicit expression can be found in [51]. As for the massive vector 
boson self-energies, terms of order mf originate from the vertex corrections to the bb final 
state [53]. In order to complete the vertex corrections, the wave function factors of the 
external fermions have to be added, which can be obtained from the fermion self-energy 
diagrams [50]. 

Last, the contribution of box diagrams for the process e+e~ ff (see Fig. 15) can 
be written in terms of two combinations of four-point functions. They arc ultraviolet 
finite, apart from the unitary gauge, and since they are not resonating at the peak, 
their effect is numerically very small, so that generally they are neglected for LEPl 
physics, although their inclusion is necessary for the gauge invariancc of the calculation. 
As will be discussed in Sect. 51.1, the box corrections become relevant far away from 
the Z° peak. 

e W,Z f e" f 

* sAAAAA^ 



« lAA/W\J « — 

e"" w,Z f e"" 

Fig. 15. - One-loop box corrections. 

It is worth noticing that, from the point of view of the determination of unknown 
SM parameters, the bosonic corrections involving the insertion of internal Higgs-boson 
lines are particularly interesting, since they introduce a dependence of the theoretical 
predictions on the unknown Higgs-boson mass niH- 

3'2.2. Calculational schemes. The tree-level Lagrangian of the SM contains in its 

electroweak part a certain number of free parameters which need to be fixed by compari- 
son with experimental data. The choice of the lagrangian parameters and their relations 
with a set of experimental data is the content of a renormalization scheme. 

The gauge sector of the minimal SM is characterized by three free parameters fixed 
by three input data points, which have to be known with high experimental accuracy, 
as they are the input for precision calculations. At present the three most precise ex- 
perimental quantities are the fine structure constant a, measured by means of one the 
methods quoted in Appendix B, the muon decay constant G^, measured through the 
muon lifetime, and the Z-boson mass measured at LEPl. Their precision makes them 
the most widely used experimental input data for the calculation of radiative correc- 
tions. Having defined the parameters, the calculation of any other physical quantity can 
be used as a test of the theory by direct comparison of the theoretical prediction for the 
observable with its experimental value. 

The relations between the lagrangian parameters and the measurable quantities de- 
pend on the order at which perturbation theory is carried on. In particular, going beyond 
the tree-level approximation, the situation is complicated by the appearance of ultraviolet 
divergences which make the differences between the tree-level parameters and the radia- 
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tively corrected ones infinitely large. Because of this, a commonly adopted procedure 
is the introduction of counter- terms: each bare parameter is split into the renormalized 
parameter and a counter-term which absorbs the ultraviolet divergences plus finite terms. 
The values of this finite parts depend on the renormalization prescriptions. Due to the 
renormalization group invariance, all renormalization schemes are equivalent. However, 
since the radiative corrections arc obtained from a finite order perturbative expansion 
(eventually supplemented by resummation to all orders of some class of gauge invariant 
contributions), the predictions for physical quantities are scheme dependent. The differ- 
ences between two calculations, performed at a fixed perturbative order in two different 
schemes, appear in higher-order terms. This fact allows to estimate the size of missing 
higher-order effects by comparing the results obtained with calculations based on different 
renormalization schemes (sec Sect. 3'6). Given the high precision of the experiments at 
the Z° peak, it is very important to keep under control these higher-order effects. In the 
literature several schemes have been adopted for the calculation of radiative corrections. 
With the exception of the low energy scheme [50,54 56], which has been used before the 
advent of LEP, it is worth mentioning the following calculational schemes adopted for 
various formulations of radiative corrections to observables at the Z° resonance: 

• the on-shell scheme [57]; 

• the JTS scheme [58]; 

• the scheme [59]; 

• the * scheme [60]. 

In the following a brief review of the various schemes is given. 

The basic idea of the on-shell scheme is the observation that Thomson scattering and 
the particle masses set natural scales where the parameters e, Mz, Mw, "niH and to/ can 
be defined. The finite parts of the mass counter-terms are fixed by the renormalization 
conditions that the particle propagators have poles at their physical masses. In the case 
of unstable particles like the and bosons the mass is not uniquely defined. For 
LEPl physics the vector boson mass is commonly defined as the zero of the real part 
of the inverse propagator. At one- loop order, after a Dyson resummation of the self- 
energy diagrams and neglecting external fermion masses, the propagator is given by the 
expression 



where the Lorentz structure has been neglected for simplicity and My is the bare mass. 
In the case of the Z-boson propagator, due to mixing eq. (49) is modified by reducible 

higher-order terms, which for simplicity are not considered here. Introducing the mass 
counter-term means the following relation between the bare mass My and the physical 
mass My: 

(50) M^^ = Mlr+6M^. 

Inserting the above expression for the bare mass in the propagator of eq. (49), the ex- 
pression of the counter-term can be easily identified as 



(51) 



5M^ = ReS'y(M^). 
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The charge counter-term is fixed by the condition that e equals the ee^ coupling con- 
stant in the Thomson limit of Compton scattering. At one-loop order, the ee^y coupling 
receives contributions only from the photon sclf-cncrgy and the mixed ^-Z transition, 
because the corrections related to the external particles cancel each other as a conse- 
quence of a generalization of the QED Ward identity. The expression for the charge 
counter-term reads 



(52) 



77 



(0) 



sw S-yz(O) 



CW 



Ml 



where sw and cw are the sine and cosine of the weak mixing angle defined in the SM by 



(53) 



sin 'dw = 1 — 



w 



Ml 



According to ref. [61], the above definition of sin^ 'd^ is assumed to be valid to all orders 
of perturbation theory. 

As a matter of fact, the precision of the PT-boson mass measurement is not adequate 
for precision physics at the resonance. For this reason the input W-boson mass is 
replaced by the more accurate value of obtained from the /^-lifetime r^: 



(54) 



1927r3 



m,. 



1 



2a rn, 
2iT V Stt rrie 



25 

T 



This equation, obtained within the effective four-fermion Fermi interaction with the 
inclusion of QED radiative corrections up to 0{a'^), is used as the definition of in 
terms of the experimental /i lifetime. Calculating the last one within the SM at one-loop 
order, the following relation can be established between and Mw'- 



(55) 



V2 



where Ar is a ultraviolet finite combination of counter-terms and loop diagrams. Since 
Ar is a function of e, M^/, Mz, mn and mj, cq. (55) can be solved itcratively for Mw- 
By inspection of the various contributions, Ar can be written in the following form: 



(56) 



Ar = Aa - 



A/)+(Ar), 



where Aa is the fermionic contribution to the photonic vacuum polarization, and Ap is 
the following ultraviolet finite combination [62,63]: 



(57) 



Ap = 



Sz(0) Evy(O) 



Ml 



which is quadratic in the top-quark mass 



(58) 



Ap = N, 



87r2\/2' 



Ap is the corrections to the p parameter, defined as p ■ 

to one in the minimal SM at the tree level [62-64]. 



M^/M| cos^ 'dw, which is equal 

takes into account non- 



(Ar) 



leading corrections, among which the most interesting ones, from the phenomenological 



36 



G. MONTAGNA, O. NICROSINI and F. PICCININI 



point of view, are due to Higgs-boson loops, yielding as leading contribution the following 
asymptotic logarithmic term (rn/f » Mw) [65-67] 



167r2 \ 3 V 6 
and to top-quark loops, yielding a logarithmic term given by 



At the one- loop order, no quadratic Higgs-boson mass corrections appear as a conse- 
quence of the so-called custodial symmetry of the Higgs-boson sector of the minimal 
SM. 

Due to the presence of large contributions to Ar of the kind a\n{Mz/mf), these 
terms are resummed to all orders by writing eq. (55) in the following form [68] : 

(61) ^ - 



V2 8s2^M2.(l-Ar)- 

This form takes into account to a good approximation also the terms of the order of 
\ia{mf /Mz) [68]. More details on higher-order terms are given in Sect. 3'2.4. 

The method outlined above of parameter renormalization is sufficient to obtain finite 
S'-matrix elements. However propagators and vertices by themselves arc not finite. To 
this aim also field renormalization has to be carried out. This allows to fulfill further 
renormalization conditions, in particular the vanishing of the mixed ■y-Z transition for 
real photons. The different ways of implementing field renormalization and the different 
gauges in which the calculations are performed are the main differences among the various 
realizations of the on-shell scheme present in the literature. 

In the MS scheme, the counter-terms are defined only through the divergent parts 
proportional to A — In of the radiative corrections to the bare parameters, as these 
are fixed only by the bare Lagrangian, and do not depend on particular renormalization 
conditions as is the case for the on-shcU scheme. The MS expressions of the parameters 
can be directly obtained replacing the bare parameters with the MS ones, and replacing 
at the same time the quantity A in the radiative corrections to the bare parameters with 
In/Lt^, fj, being the renormalization scale. In so doing, the MS parameters depend on the 
arbitrary scale /x, which can be naturally chosen to be Mz for electroweak calculations. 
The MS electroweak mixing angle can be defined as 

(62) §2^1-%, 

Mf 

where Mw,z are the MS vector boson masses. 

According to refs. [69] and [70], a way to link directly to the experimental inputs 
a, G/j,, and Mz is to calculate the radiative corrections to fj, decay within the MS 
framework introducing the correction Arw together with a generalized relation between 
the W- and Z-boson physical masses: 

(63) 5^(1 - Afw) = 



(64) M^ = c?pM\ 



V2G^Ml' 
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where = 1 — A/5 . 

Within the MS scheme the relation giving the Mw-Mz interdependence can be 
expressed as 



The advantage of adopting the MS scheme is twofold [71]. First, this calculational 

framework leads to effective couplings which absorb the largest part of the top-quark 
and Higgs-boson mass dependence of the radiative corrections. Secondly, the knowledge 
of the gauge coupling constants at the Z-boson mass can be naturally extrapolated to 
higher energies in order to test scenarios of Grand Unification. 

The idea of the scheme is to directly relate the bare parameters to the input 
data a, and Mz- This can be realized according to different methods, for example 
by assuming on-shcU parameters [72] or MS parameters [51,73]. Here a brief accoimt of 
the second realization is given, where the bare parameters are chosen to be the SU{2) 
coupling constant g, the T4^-boson mass Mw and the squared sine of the weak mixing 
angle s^. Starting from the observation that there is no one to one correspondence 
between lagrangian parameters and experimental data, in this scheme no attempt is 
made to define renormalized parameters (in particular s|) beyond lowest order. A set 
of three fitting equations [51,63] for the bare parameters is introduced, corresponding to 
the definition of a from Thomson scattering, from the ji lifetime and Mz from the 
zero of the real part of the inverse propagator: 



In eq. (66), df^ are the theoretical expressions of a, G^ and Mz in terms oig, Mw, and 
A, whereas d1^^ are the corresponding experimental values. The equations can be solved 
to first order in perturbation theory with respect to the bare parameters which contain 
ultraviolet divergences represented by the quantity A. They can be properly modified 
to account for resummation of relevant gauge invariant higher-order effects [51]. In the 
calculation of any physical quantity by means of one-loop diagrams calculated with tree- 
level parameters and tree diagrams calculated with up to one-loop order parameters the 
ultraviolet divergences cancel. This procedure is suitable for an order by order numerical 
renormalization. In the numerical calculation, the quantity A enters as an arbitrary 
parameter which can assume any numerical value without changing the final result. 

The condition = is fulfilled through a proper redefinition of the bare cou- 
pling , which automatically guaranties that the sum of all Zff vertices is ultraviolet 
finite [51]. 

It is worth noticing that in the scheme the PF-boson mass plays the same role as 

any other physical observable, contrary to the on-shell scheme where at the same time 
it appears as a prediction and as an intermediate parameter for the calculation of other 
observables. 

In the * scheme, an effective lagrangian to one loop is introduced, and the universal 
effects of oblique corrections resummed to all orders are absorbed in the following three 
running parameters (free of ultraviolet divergences): ei{q^), s1{<f') and G^^(q'^). The 
values of these three starred functions have to be fixed at a certain set of q a and then 
can be calculated at any other scale by means of a set of evolution equation. The choice 



(65) 




with A = (7ra/V2G/,) ^ . 



(66) 



dr^' = rff(,g,Mw,s^,A), 



(i = 1,2,3). 
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of ref. [74] is to fix el and G^^ at = by means of a and G^, and si at the ^-boson 
mass scale by means of the pole of the Z° propagator. In order to cure the problem 
S-yz(O) 7^ in the bosonic sector, a rcdiagonalization of the neutral current sector in the 
one-loop lagrangian is performed including in the self-energies the contribution of the 
universal and gauge dependent parts of the vertex corrections. 

Besides the calculational schemes described above, one should also mention the Z- 
peak subtracted representation of e+e" ff processes at one loop, recently suggested 
in [75] . In this approach, the input parameter is replaced by quantities measured on 
top of the Z'^ resonance. For instance, for leptonic final states the new input parameters 
are the leptonic Z-boson width F; and the effective sine sin^iJ^y^j, that reabsorb the 
bulk of loop effects, while the residual one-loop corrections are contained in quantities 
subtracted at = M"^. This approach has been proved to be particularly powerful for 
the description of e+e" annihilations above the Z° peak, especially if an investigation of 
models of new physics is the final goal. 

3'2.3. Transition amplitudes and effective couplings. The basic ingredients of loop 
calculations described in the previous subsections can be used to evaluate transition 
amplitudes for the physical process e+e~ ^ //. Before discussing the radiative corrected 
amplitude, it is worth recalling the general structure of the tree-level amplitude which, 
excluding for simplicity the case of Bhabha scattering, reads [76] : 

(67) M ~ i [QeQfla ® 7" + Xlaig', - 91ib) ® 7"(5^ - sfls))] , 

where % is the propagator ratio 

s 

^^^^ s-Ml + iTz{s)Mz' 

and the symbol F^ ® F" means product of spinorial currents. The s-dependent width 
Tz{s) ~ sF^/MI appearing in eq. (68) is due to the imaginary part of the Z self energy. 
In eq. (67), the photon and the Z-boson exchange diagrams are well separated and 
moreover the Z° contribution is written in a factorized form. By virtue of these two 
properties, simple relations between the pseudo-observables defined in Sect. 3'1 and the 
measured quantities can be established. 

However the two features are lost when the electroweak non-photonic corrections 
(self-energy diagrams, Zff vertex insertions and weak boxes) are considered. In this 
case the matrix element can be written in the following way: 

Mr^-^ [a(s)7a (8 7" + X {Klis, t)ic. ® 7" - Kli^, t)lo. ® 7^75 

(69) - F^i (s, f)7a75 ® 7" + is, i)7a75 ® 7^75)] , 

where a{s) is the QED running coupling constant with only fermionic contributions (in 
order to have a gauge invariant term) and the form factors F^j are complex valued 
functions of the input parameters and of the kinematical variables s and t. Usually the 
self-energies contributions are rcsummed in their fermionic component to take into ac- 
count large higher-order reducible diagrams. The bosonic parts are instead expanded at 
0{a), in order to preserve gauge invariance. The dependence of the form factors on t 
is only due to the presence of box diagrams. Being non-resonant at the Z° peak, their 
numerical effects are very small and can be neglected, even if from a theoretical point 
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of view the gauge invariance is no more respected. Other non-resonating contributions, 
such as the bosonic insertions to the photon propagator and photon-fermion vertex cor- 
rections, can be safely neglected. Another commonly used approximation is the Z'^ pole 
approximation, which means to fix the scale in the form factors at the value M\, since 
the weak corrections depend very mildly on the scale around the peak. After these 
approximations the factorization of the Z^ contribution is re-established and the form 
factors take on the form: 

(70) Fll{Ml) = G\{Ml)Gl{Ml). 

This means that the Z-boson part of the corrected amplitude can be obtained from the 
tree-level expression by replacing the vector and axial vector couplings with the corrected 
versions defined above by eq. (70). The imaginary parts of the couplings G'{(M|) are 
generally small with respect to the real parts, so that the following effective couplings 
are commonly introduced: 

(71) 9{r,A = ReG^,„(M|). 

As already noticed in Sect. 3"1, the decay width of the Z° boson into a // pair is 
given by the following expression in terms of the effective couplings: 

(72) = 4Ar,ro[(ff^)^4 + {g^fR^A- 

In the literature other sets of parameters replacing the effective couplings are present. 
For example according to ref. [66] the parameters p/ and fc/ are introduced related to 
Qy and by the following relations: 



(73) 


^ = l-4|g/|s2^A:/ 
9 a 


(74) 


with 




(75) 


Pf = 1 + %, 


(76) 


kf = l + 5kf. 



By means of these parameters, the Z^ partial widths can be written in the following way 
(neglecting for simplicity final-state QED and QCD interactions and mass effects): 

(77) Tf = T^NlpfWj - 2Qfsl,kff + 1]. 

The above described procedure of step-by-step approximations, leading from the exact 
one-loop amplitude to an approximate amplitude formally identical, or very similar, to 
the Born-approximation one, but written in terms of form factors evaluated at the Z^ 
peak, is known as Improved Born Approximation (IB A). Some different realizations of 
IBA's are available in the literature, and can be found for instance in [66,77,78]. 

3'2.4. Higher-order electroweak corrections. Given the high precision reached in 

the experimental measurements, particular attention must be paid to include in the 
electroweak corrections potentially large higher-order effects. These generally originate 
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from terms containing logarithms of the type ln(Mz/m/), where m/ stand for a generic 
fermion mass, or from contributions proportional to the top-qua,vk mass. In Sect. 3'2.2 
the prescription based on renormalization group arguments [68, 79] 

(78) 1 + 

has been introduced to include to all orders the leading logarithms of 0(a" In^ {Mz/ruf)), 
contained in Aa of eq. (55), giving the interdependence between Mw, Mz and Gf^. 
However, with this prescription the powers (Ap)" are not correctly resummed. In ref. [80] 
it has been shown how to take into account in the on-shell scheme the leading two-loop 
contributions proportional to mf by means of the replacement 

(79) - -r > - T • "2 h Arremt 

1-Ar 1-Aa i + f^Ap 
where Ap includes the contribution of two-loop one particle irreducible diagrams: 



(80) Ap = "ixt 



1 , (2) / mH 



mt 



Xt = 



87r2V2' 



The two-loop correction p^^) was first calculated in ref. [81] in the limit of light Higgs- 
boson mass tuh « mt, yielding p^^^ = 19 — 27r^. More recently the calculation has 
been carried out by several groups [82 84] for arbitrary values of the ratio mfz/m^ in 
the limit of vanishing gauge coupling constants. In ref. [68] it has been shown that the 
fermionic mass singularities of 0(a^ In (M^/m/)) are correctly taken into account by 
keeping Arrem in the denominator, i.e. 

1 1 

(81) 



1 - Ar 



(1 - Aa) f 1 + ^Ap) - Ar 



The effects of two-loop heavy Higgs-boson contributions have been investigated in 
refs. [85] and [86]. 

Very recently the two-loop sub-leading corrections of ©(a^m^/M^) have been cal- 
culated for the vector-boson masses correlation [87-89], and found to be potentially of 
the order of the leading contributions, depending on the Higgs-boson mass. More- 
over, the exact two- loop Higgs-boson mass dependence has also become available [90], 
introducing a shift in the M^-boson mass contained within 4 MeV [91] with respect to 
the calculation of ref. [88,89]. 

Another place where higher-order electroweak effects play an important role is the 
partial width Z //. The effective couplings contain leading universal parts, arising 
from self-energies and counter-terms, and flavour-dependent parts coming from vertex 
corrections. The leading universal contributions are given by [76] 

(82) ((5p/)„„i„ = Ap-h..., 

(83) (%)„„i„ = 

The higher-order terms can be incorporated by means of the following replacements: 

(84) + 

(85) fc/^l + 
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with Ap given by eq. (80). 

The partial decay width Z ^ bb contains an additional dependence due to vertex 
diagrams with virtual fop-quark lines, resulting in additional leading terms of the order 
G^m| for the effective couplings pb and kb- The additional leading two-loop electroweak 
effects of the order G^^mj have been computed in refs. [82] and [83], and amount to the 
following redefinition of the effective couplings: 



(86) 
(87) 

where Tf, is given by 
(88) 



n 



Pb 
kb ■ 

-2xt 



Pb{l+Tbf, 

kb 
l + n 



1 + xtr' -4- 



m 



H 



where (171^/171%) can be found in [82,83]. 

At the level of two-loop or higher-order corrections, other phcnomcnologically rele- 
vant contributions come from the mixed electroweak-QCD corrections, discussed in the 
next Section. 

3'3. QCD Corrections. - QCD corrections play a significant role as higher-order 
corrections due to virtual gluon insertions in electroweak loops (mixed electroweak-QCD 
corrections) and as final-state corrections in the hadronic decay channels. With the 
exception of very recent improvements, an account of perturbative QCD calculations for 
the Z-boson observables can be found in [14] and [92]. 

3'3.1. Mixed clcctroweak-QCD corrections. The need of very accurate predictions 
demanded by the high precision of the measurements on the Z'^ peak, motivated various 
authors to perform calculations of higher-order effects beyond the one-loop and higher- 
order electroweak corrections already discussed in Sect. 3'2. Great effort was spent in 
the calculation of corrections due to the insertion of virtual gluons in the electroweak 
quark loops, known as mixed electroweak-QCD corrections, giving rise to contributions 
simultaneously depending on the weak coupling and the coupling of the strong inter- 
action ag. Higher-order corrections of this kind have been calculated both for W- and 
Z-boson self-energies, thus affecting all Z ^ ff decay channels, and specifically for the 
peculiar Z ^ bb vertex. 

QCD corrections to the gauge-bosons self-energies have been computed in the SM 
considering the insertion of one and two gluons into the internal quark loops. Typical 
Feynman diagrams contributing to these corrections are depicted in Fig. 16. 





Fig. 16. - Example of Feynman diagrams for QCD corrections to gauge-bosons self-energies. 



The diagrams involving the exchange of a single gluon, described for instance by the 
first graph in Fig. 16, give rise to two-loop corrections of the order of aus, that are now 
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exactly known [93] and are dominated by terms of order G^asmf. The insertion of two 
virtual gluons, of the kind depicted in the second graph of Fig. 16, generates a three-loop 



contribution of the order of 



that is partially under control and introduces effects 



of the form G^a^mf [94]. In formulae, the effect of such corrections to the Ap factor 
introduced in Sect. 3'2 can be cast, in the heavy to;>-quark limit, as follows 



(89) 



Ap^ A/?-(1 + (5/3qcd) 



where dpqcD is the QCD correction to the leading G^mf term 



(90) 



SpQCD = -2.86as - 14.6a^, 



as = 



the running coupling constant being evaluated at the energy scale given by the top- 
quark mass. It is worth noticing that the QCD correction to the p parameter is negative, 
and tends to screen the electroweak contribution. (^) 

A second source of higher-order mixed effects comes from QCD corrections to the 
Z ^ bb electroweak vertex once corrected by the one-loop contributions depicted in 
Fig. 14 (see Sect. 3'2). Due to the presence of virtual top-quark lines in the one-loop bb 
vertex, significant non-universal mixed corrections, depending on as and mt, are addition- 
ally present for the bb final state. The effects induced by 0{as) gluon radiation introduce 
two-loop corrections to the leading electroweak term of the order of G^asWj [96] and to 
the hi{mt/Mw) term of the order of G^^as lii{mt/Mw), the latter with a very small nu- 
merical coefficient [97]. The missing ncxt-to- leading corrections of 0{aas) to the partial 
width Tb have been very recently calculated [98], thus improving the SM prediction for 
the Z ^ bb vertex. 

Finally, final-state mixed corrections to the Z ^ qq decay channels, due to the in- 
terplay between QED and QCD radiation off final-state quarks, are exactly available [99] 
and taken into account in the theoretical predictions. For example, the decay width 
into qq final states receives a contribution from these mixed corrections given by the 
factor 



(91) 



2 3 a{s) 



TT 



^ _ 1 as{s) 

3 TT 



with s 



Ml, 



that can be easily deduced from the QCD results for the Z 
in the next Section. 



qq decay width, discussed 



3'3.2. Final-state QCD corrections. In addition to the classes of radiative corrections 
discussed in the two previous Sections, an adequate theoretical description of the pro- 
cesses explored at LEP requires, as a further basic ingredient, the inclusion of the QCD 
corrections to the Z ^ qq decay channels, arising from the emission of real and virtual 
gluons off the final-state quarks (see Fig. 17). Among the observables measured at the 
Z° peak, these corrections affect the total hadronic cross section ahad, the FB asymmetry 
of heavy-quark {q = c, b) production A'p'^ , as well as the total hadronic width Th and 
the partial widths into c and b quarks r^b. When discussing final-state QCD effects 
in hadronic Z-boson decays, it is necessary to distinguish between corrections to the 
production of light, or massless, quarks and corrections to the heavy-quark production. 



(^) The SUSY-QCD corrections to Ap have been recently computed; they could introduce a shift of the 
order of 10 MeV on the W-boson mass [95]. 




Fig. 17. - Feynman diagrams for final-state 0{as) QCD corrections. 



In the case of light quarks (q = u, d, s) with iriq ~ 0, the QCD final-state correction 
to the cross section (or, equivalcntly, to the partial decay width) is completely known up 
to 0{al). This high perturbative accuracy is of course of the utmost importance in view 
of a precise determination of the strong coupling constant in Z-boson decays, as it will 
be shown in the following. The result for the Z° partial widths into massless quarks can 
be written as 

(92) r = To 
with 

(93) To = (A^c = colour factor) , 

gy and being the effective electroweak couplings. Kqcd is the up to three-loop 
correction factor for massless quarks given by [100, 101] 

(94) Kqcd = 1 + ^ + 1.41 (^Y - 12.8 (^Y . 

TT V TT / V TT / 

Actually. leading and next-to-leading second-order top-quarh mass corrections to the qq 
decay width in the massless limit are also available in the literature [102] and are taken 



QCD 



44 



G. MONTAGNA, O. NICROSINI and F. PICCININI 



into account in the more refined (but too lengthy to be shown here) QCD radiation 
factors usually implemented in standard computational tools for LEPl physics. The 
complete expressions, together with a discussion of the numerical effects of the various 
contributions, can be found in the recent review paper of ref. [92]. 

For massive quarks the situation is different, due to the presence of finite mass 
terms. In such a case, the vector and axial-vector contribution to the Z-boson decay 
width receive different corrections, known at different perturbative orders. Actually, the 
calculation of the axial- vector part of the hadronic .Z^-boson decay rate is more involved 
than that of the vector part. This is because the heavy quark does not decoiiple in the 
axial- vector part and hence one cannot avoid calculating massive diagrams. In particu- 
lar, the bb final state receives peculiar contributions from top-quavk dependent two-loop 
diagrams in the axial-vector part. Consequently, the formula of eq. (92) in the massless 
limit needs to be integrated with additional correction factors for massive quarks, of the 
form 



(95) 



r = ro 



The coefficients in the perturbative expansions 



(96) 
(97) 



Rv 
Ra 



TT \ TT / V TT / 



1^2 
TT 



+ ..., 



depending on the heavy-quark masses and m^, are calculated up to third order in the 
vector part and up to second order in the axial- vector part [92, 103 105]. Given the 
above formulae for the quark partial widths in the massless limit and massive case, 
the total hadronic width and cross section can be obtained as the sum over the partial 
contributions of each qq channel. 

Recently, the non-factorizable part of the two- loop 0{aas) correction to the hadronic 
Z-boson width has been calculated, introducing an extra negative contribution of about 
-0.6 MeV [106]. 

Concerning the forward-backward asymmetry for heavy quarks, the QCD corrections 
to ApB are available in the literature [107]. At the order as, they yield a correction of 
the kind 



(98) 



where k depends on the mass of the heavy quark, and its explicit expression can be found 
in [108]. 

Last, a few comments are in order. First, since the above QCD factors depend on 

the perturbative expansion parameter as and the latter depends on the energy scale 
according to the renormalization group equation for the /? function, a formula for the 
running of ag needs to be specified. Indicating with the renormalization scale and 
with nj the number of active quark flavours of mass rUq <C /x, the solution for ai'^\lj) 
typically used reads [92, 104] 



(99) 



"s(m) = Try 



1 - 



1 f3i\nL 



|(ln^.-ln.-l).| 
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with 

/3o = (ll - luf^ /4 , /3i = (^102 - y n/) /16 , 
« /^2857 5033 325 

In the above formulae L = ln(/i^/A|^). Actually eq. (99), together with the coefBcients 

of eq. (100), is valid in the modified minimal subtraction (MS*) scheme and corresponds 
to the three-loop expansion of the QCD (3 function. It is used with n/ = 5 for the num- 
ber of active flavours and it is technically called next-next-to- leading solution of as{p). 
At present, also the four-loop coefficient in the MS scheme is available [109], thanks 
to the recent calculation of the coefficient of the QCD /3 function [110]. A second 
comment concerns the mass correction terms entering the QCD factors for heavy-quarks. 
Analogously to the coupling constant, the quark masses are running parameters in QCD 
and obey the renormalization group equation controlled by the anomalous dimension 
7m,, that is known up to three-loop accuracy. The masses appearing in the formulae for 
c- and 6-quark production have thus to be understood as running (p, dependent) quark 
masses. The relations used to account for such effect read as follows [14]: 

m(/i) = m(m )exp ^ - / dx } , 

(101) m, = m,iml) |l + + K,^^ + Oial)], 

where rriq is the so-called pole mass (defined in quantum field theory as the position of pole 
of a renormalized quark propagator), rf^ is the MS running mass and Kc.t = 13.3, 12.4. 
Notice that the running mass depends on a^. Using, for instance, for the 6-quark pole 
mass nib = (4.7 ± 0.2) GeV and as{Mz) = 0.12, one obtains a running mass at the 
peak of rnb{Mz) — 2.8 GeV, leaving non-negligible mass effects. The running c-quark 
mass rfic is about a factor of five smaller than Mb, thus making c-quark mass corrections 
almost invisible in the Z-boson decays [92] . 

Before turning to QED corrections, it is worth quantifying the effect of the elec- 
troweak and QCD corrections discussed up to now. To this aim. Fig. 18 shows the 
comparison between the total cross section and the forward-backward asymmetry for the 
process e~^e~ 1^-^ fJ'~ at the tree level (dashed line) and the full prediction including all 
the one-loop and the relevant higher-order electroweak and QCD corrections (solid line). 
In Fig. 19 relative deviations for the cross section and absolute deviations for the asym- 
metry arc shown. The dotted lines represent the difference between strictly tree-level 
predictions and "tree-level like" results, obtained by running the QED coupling constant 
(and hence sin^ 'dw) and including final-state QED and QCD corrections in the total 
.Z-boson width in the propagator. The dashed lines represent the deviations of the last 
prediction from a full electroweak/QCD calculation. The full electroweak/QCD results 
have been produced by means of TOPAZO in its default mode, neglecting theoretical un- 
certainties (see Sect. 3'6). As can be seen from the figures, the bulk of the corrections is 
due to the running of the QED coupling constant and to final-state QCD corrections to 
the Z-boson width. However, there is an effect of 0.5 -j- 1% on the cross section and up 
to 0.002 on the asymmetry due to non-trivial corrections as computed by means of a full 
electroweak/QCD calculation. 
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Fig. 18. - Comparison between the cross section and the forward-backward asymmetry for 
e^e~ fj,^ , as computed in the Born approximation (dashed line) and including all the 
one- loop and the relevant higher-order electroweak and QCD corrections (solid line). Numerical 
results for this last case by TOPAZO [111]. 



3'4. QED Corrections. - As already said in Sect. 2'3, QED corrections originate from 

those diagrams with extra real and/or virtual photons added to the tree- level graphs (see 
Figs. 22, 23 and 24 below). Although these corrections are not particularly interesting 
with respect to the underlying theory, they need a special attention at LEP, as at any 
other e~^e~ collider, because the size of their effects strongly depends on the experimental 
cuts, and therefore their proper treatment constitutes the unavoidable link between the 
data taking and the physics analysis [43]. In fact, they are large corrections at high 
energies because, as already explained in Sect. 2'3, they are dominated by logarithmic 
contributions of the kind {a/Tr)L, L = ln(s/mg) being the so-called collinear logarithm. 
At LEPl, where s ~ M|, L is of the order of 25, and the effective expansion parameter 
in perturbation theory is (3 = 2a{L — Vj/tt ~ 0.1 rather than a. The Breit-Wigner 
line-shape of the Z° resonance in Born approximation is sensibly modified by the QED 
corrections for the following typical effects [112, 113]: 

• the peak height is lowered by around 25%; 

• the peak position is shifted towards higher energies by around 100 MeV; 

• a hard radiative tail, that increases the lowest-order cross section, appears above, 
say, 93 GeV. 

The significant distortion introduced by the QED corrections on the Z'^ line-shape 
of the process e~^e~ hadrons can be clearly seen in Fig. 20. 

Also the forward-backward asymmetry, ApB, which crosses the zero in the proximity 
of Mz, is significantly affected by the QED corrections, that shift the position of its zero 
and globally change its shape [114, 115]. The QED effects on ApB for the reaction 
e+e~ — > are illustrated in Fig. 21. 

The whole set of QED corrections to s-channel e+e" annihilations can be divided 
into three subsets, each separately gauge-invariant: 
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Fig. 19. - The left windows show the bulk of electroweak/QCD corrections as obtained by 
means of a "tree-level like" calculation (dotted line) and the residual effect of a complete elec- 
troweak/QCD calculation (dashed lino). The right windows show a blow-up of this last differ- 
ence. Numerical results for the full electroweak/QCD prediction by TOPAZO [111]. 



• initial-state corrections; 

• final-state corrections; 

• initial-final state interference. 

The initial-state corrections arc by far the dominant ones, because, at a difference 
from final-state and initial-final state interference contributions, they are responsible for 
a reduction of the cm. energy available for the hard-scattering reaction. 

3'4.1. Initial-state radiation. The Feynman diagrams for 0{a) initial-state QED 
corrections to a generic s-channel process e+e~ ^ 7, — > // are depicted in Fig. 22. 

Since at LEPl energies the electron-positron collision takes place at energies very 
close to the Z-boson mass, Mz, the radiation emitted by the initial-state (ISR) is es- 
sentially radiation of soft photons. This remarkable property is a direct consequence of 
the fact that the maximum energy of the photons emitted by the initial-state can not 
exceed the ratio Smax = ^z/Mz, with Smax — 0.03. In fact, the finite Z-boson width, 
Tz, acts as a natural cut-off that strongly inhibits the emission of those hard photons, of 
energy fraction larger than Emax: that would prevent the formation of the Z'^ resonance. 
This soft-photon dominance at the Z° peak is a very useful guideline when facing the 
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Fig. 20. - The effect of QED corrections on the Z'^ line-shape of e^e^ hadrons. The dashed 
hne is the QED-dcconvoluted cross section, the sohd hne is the QED corrected one. Numerical 
results by TOPAZO [111]. 



problem of including finite order perturbative corrections, beyond the leading logarith- 
mic (LL) approximation, that are actually demanded by the high-precision level of the 
LEPl measurements. 

As recalled in Sect. 2'3, the n-th order contribution to a QED corrected cross section 
a'"^ can be cast in the form 

n 

(102) -^"^ = (^"E4"^^^ 

fc=0 

where H = C{s) = L — 1. The coefficients a^"' in eq. (102) are given by the following 
expression 

k 

(103) 

3=0 

where / is the so called IR logarithm given by I = \n{E/AE), AE being the maximum 

energy of the emitted photons. This general expansion obviously applies to the IS QED 
corrected cross section. For example, the up to 0{a) cross section in the soft- photon 
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Fig. 21. - The effect of QED corrections on the forward-backward asymmetry Afb of e'^e~ 
jj,^ around tlic Z" peak. The dashed line is the QED-deconvoluted asymmetry, the soUd Une 
is the QED corrected one. Numerical results by TOPAZO [111]. 



approximation, as obtained by standard diagrammatic techniques, reads 

(104) a^^\s)=ao{s)l^l + ^ 

which shows a large 0{a) negative correction, clearly indicating that higher-order QED 
corrections are required for a high-precision reconstruction of the line-shape. These 
higher-order effects are kept under control at all orders by employing one of the ap- 
proaches described in Appendix A. For instance, in the structure function (SF) approach 
the LL-corrected cross section, in the extrapolated set-up, can be written as 

(105) cr(s) = / dxH{x, s)ao ((1 - x)s) , 

Jo 

H{x,s) being the QED radiator of eq. (186), and ao{s) being the lowest order cross 
section, possibly including all the short-distance process dependent corrections, as the 
electroweak and QCD ones previously discussed. This allows to exactly reproduce the 
LL content of eq. (102) at any order and, in particular, that of eq. (104) at 0{a) (see the 
series expansion in the LL approximation given by eq. (12) in Sect. 2'3 ). The sub- leading 
terms present in eq. (104) are not naturally reproduced by the perturbative expansion 



2l{L-l) + ^L+^-2 
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Fig. 22. - The Feynman diagrams for 0{a) initial-state QED corrections to s-channel process 



of the cross section obtained via the algorithms for the calculation of universal photonic 
corrections. It is worth noting that they are numerically important for a description of 
the Z'^ line-shape with a theoretical accuracy at the 0.1% level. For this reason, they are 
usually incorporated in the algorithms for the; resummation of LL corrcc;tions, by means 
of a proper matching with the perturbative diagrammatic results. For instance, in the SF 
approach, a standard procedure consists in replacing the Gribov-Lipatov form factor in 
front of the exponentiated term of the SF's {(3/2) (1 — x)^~^ by means of a soft-|-virtual 
iiT-factor A^_|_y given, up to the first order in a, by 



(106) 



A's+y(s) = 1 



a 
2^ 



-L 



in such a way that the up to 0{a) cross section of eq. (104) is exactly reproduced. 

By virtue of the soft-photon dominance in the ISR mechanism previously discussed, 
the hard-photon non-leading contributions, left-over in the above derivation, turn out 
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to be numerically unimportant at LEPl. For the applications to precision physics at 
the peak, the inclusion of sub-leading terms in the SF approach is pushed up to 
C(a^) [112], relying upon the exact second-order calculation of initial-state QED correc- 
tions to s-channel processes available in the literature [41,116,117]. By exploiting again 
the soft-photon dominance, one can extract from the second-order complete calculation 
the soft+virtual contributions only, so that the up to O{o?^ ii'-factor, to be placed in 
front of the electron SF, can be cast in the form 

(107) ^s^y{s) = 1 + A.^^Us) + A,f^{s) 
and is calculable through the relation 

(108) A5+y(,s) = (A'5+v-(.^))'-^/3'. 

The quantity A5_|_y(s) contains the non- leading non-IR sensitive diagrammatic correc- 
tions to the cross section corrected by the inclusion of soft -|- virtual photons up to C(a^), 
analogous to the one shown in cq. (104). Equation (108), that relates As+v{s) and 
^'s+vi^)' '^^^ derived from the explicit integration of eq. (172) in Appendix A in the 
soft-photon approximation. By inspection, A5_|_y(s) is given by 

(109) AsM^) = 1+0 ^slvis) + {^yAifyis), 

where the perturbative corrections Ag'^y{s) and A^g_^y{s) are explicitly given by 



Ai%{s) = ^-L + ^ 



3 

2" ' 3 
9 w 



(110) Aillis) = r--^] L'+ -^ + 11^^ + 3C(3) 



3 



45 11 



16 12 



L + constant terms. 



In eq. (110) C is the Riemann function, with ({3) ~ 1.202. With respect to the exact 
second-order calculation, 0{a^) sub-leading non-soft corrections are neglected, but they 
are negligible in the resonance region, their contribution to the cross section being at the 
level of 0.01% [113,118-120]. The above described procedure for matching LL results 
with exact finite-order calculations is valid in the soft-photon approximation. It is worth 
noticing, however, that a general prescription is also known, valid for any ES, i.e. taking 
into account also hard-photon contributions, as can be found in [23]. 

A complete treatment of IS C(a^) QED effects to the Z° line-shape requires, as a 
last ingredient, the inclusion of the corrections arising from the conversion of a photon 
into leptonic pairs and/or hadrons (pair corrections). These contributions, which are 
dominated by e+e^ pairs, can be included using the 0(0^) formulae available in the 
literature [121]. Their main effect is a reduction of the peak cross section of around 
0.3%, hence significant in the light of the remarkable experimental precision of the LEPl 
data. 

Equation (105) allows to understand the main effects of ISR on the Z° line shape. 
Actually, since AE/E is of the order oiTz/Mz around the resonance, the explicit calcu- 
lation of the integral appearing in eq. (105) in the soft-photon limit leads to the following 
reduction of the peak height: 
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On the other hand, eq. (105) is nothing but a weighted average of the lowest-order cross 
section for cm. energies smaller than the nominal one, the weight being the radiator 
H itself. This leads to a shift of the corrected peak position towards higher energies by 
about 100 MeV, and to the raising of a radiative tail above the resonance (see Fig. 20). 
It is worth noticing, in particular, that the contribution of soft-photon exponentiation 
to the shift of the peak position is of the order of 15 MeV, and thus phenomenologically 
relevant in view of the experimental precision of the determination of the Z-boson mass. 

Before concluding the discussion on ISR, it is worth noticing that the analogue of 
eq. (105) for the corrected forward-backward asymmetry in principle does not apply. 
Indeed, since ApB is a less inclusive quantity than the total cross section, kinematical 
effects can show up, thus introducing the need for improving eq. (105). However, at the 
resonance these effects are negligible due to soft-photon dominance, and anyway they 
can be naturally taken into account in a more appropriate SF formulation, as described 



3'4.2. Final-state radiation. The final-state QED corrections, although numerically 
smaller than the IS ones under many typical experimental conditions, are an essential 
ingredient for a high-precision phenomenology of the LEPl processes. They must be 
included in any formulation of physics that aims at a theoretical precision at the 
level of 0.1%. The Feynman diagrams for 0{a) final-state QED corrections to a generic 
s-channel process e+e" ^ 7, Z° — > // are depicted in Fig. 23. 

With respect to the emission of photons by the initial state, the mechanism of 
bremsstrahlung by the final-state particles is no more characterized by the soft-photon 
dominance discussed above, simply because the radiation process occurs after the forma- 
tion of the Z° resonance. This implies that the control on final-state QED corrections at 
the level of sub-leading and constant terms necessarily requires a full perturbative calcu- 
lation, including hard photon bremsstrahlung. A treatment of FSR completely flexible 
with respect to any kind of ES requires the usage of numerical (Monte Carlo) techniques 
in order to perform the phase space integrations in the presence of arbitrary cuts. How- 
ever, when considering particular experimental selection criteria actually adopted by the 
LEP Collaborations, it is possible to obtain the final-state 0{a) QED correction accord- 
ing to a completely analytical procedure [123, 124]. In this case, the final-state corrected 
cross section can be obtained by calculating the exact matrix element associated to the 
gauge-invariant set of Feynman diagrams describing final-state real photon radiation in 
e+e~ annihilation (see Fig. 23), integrating it over the phase space allowed by cuts and 
matching the obtained result with the soft-|-virtual correction. For example, if one as- 
sumes that only an invariant mass cut of the type M^(//) > sq is present, then the 
above strategy leads to the following analytical correction factors [78, 122, 124, 125] 



in [122]. 




x+-x'^ + 2 ln(l 



x) Lf 



(112) 




\nx + 2 ln(l 



x) (In X — 1) 




+ 2 ln(l - x) (In X- 1) - 2x 
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Fig. 23. - The Feynman diagrams for 0{a) final-state QED corrections to s-channel process 



(113) +2Li2(a;)}, 

where Lf = ln{s/m'j), x = sq/s and {F ± B) denotes the forward ± backward cross 
sections defined in eq. (22). In eq. (112) and (113), lA2{x) is the dilogarithm function. 
It is worth observing that when no cuts are apphcd {i.e. in the hmit .x ^ in the latter 
equations) the final-state correction factor for the total cross section reduces to 3a(5j/47r, 
that is the well-known inclusive final-state correction, while the asymmetry does not get 
any effect. The very small value (0.17% for fermions with unit charge) of the final- 
state correction to the fully inclusive cross section is a consequence of the cancellation of 
mass and collinear singularities established by the Kinoshita-Lee-Nauenberg theorem [20]. 
Besides an invariant mass cut, also the effect of acollinearity and/or energy threshold cuts 
can be treated analytically, and the expressions for the corresponding 0(a) correction 
factors can be found in the literature [124]. The availability of such results in analytical 
form is of utility for physics at LEPl in order to test Monte Carlo programs and to 
develop fast computational tools for fitting realistic observables. 

Whenever particularly tight cuts are imposed, the treatment of FSR at the 0{a) 
can become inadequate and a procedure of resummation of the LL contributions needs 
to be advocated [126]. The leading terms can be quite easily identified either analyzing 
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the explicit perturbative results or invoking one of the algorithms for the universal pho- 
tonic corrections described in Appendix A. This allows to see that final-state leading 
contributions are of the form 

(114) 2-Q^ln(l-so/5)(L/-l), 

TT 

so that they can be extracted from the finite order correction and summed up to all orders 

according to the; preferred resummation technique. By using, for instance, a "naive" 
exponentiation procedure, the final-state correction factor, including higher-order leading 
contributions and exact sub- leading terms at 0{a), can be cast in the form 

(115) = Ay%{s) exp {2 ^ Q} ln(l - so/s) {Lf - 1) } + 

where Ay"g(s) is the soft-l-virtual X-factor of eq. (106) (provided that the substitution 
m-e — > TOy: is performed in the coUinear logarithm) and (5^"''' is the residual correction 
factor including sub-leading terms as obtained after depuration of leading logarithms 
from the complete 0{a) results. 

All the above discussion and formulae concerning FSR are valid for a non-calorime- 
tric ES (see Sect. 2'3 for more details) according to which the energy of the final-state 
fermions measured by the experimental apparatus coincides with the energy of the "bare" 
particles, regardless of the emission of collinear final-state photons. However, when 
considering the process of large-angle Bhabha scattering, this assumption turns out to 
be unrealistic. Actually, what is detected in the real environment is an electromagnetic jet 
of half-opening angle 5c, where 5c is an experimental parameter describing the resolution 
power of the calorimeter. This electromagnetic effect can be analytically accounted for 
by adding to the 0{q) final-state correction for "bare" electrons the contribution due 
to a hard photon of energy fraction greater than 1 — x, where x is sq/s for an invariant 
mass cut and 2i?o/-\/s for an energy threshold cut, and collinear with the final fermion 
within an angle < i?-^ < Jc- For electrons in the energy regime of LEP the contribution 
reads [127] 

Fcoii = 2 — C, 

TT 



C = - ln(l - x) [in (1 -h r^a;^) - 1 

;i - (1 -1(1-0=))'] In (l + rV) 



4 

+ y + 4-2^^-^^ + 4^^-^) 

(116) -h21na;ln(l-a;)-h 2Li2(l-a;), 

where 



(117) 



5c\/s 



2me 

rrie being the electron mass. Equation (116) holds under the conditions 5c rad and 

r ':$> 1, that are both very well satisfied at LEP, where 5c is of the order of a few degrees. 
It is worth stressing that the effect of the calorimetric measurement, superimposed over 
the correction for "bare" final-state particles, depends very critically on the energy or 
invariant mass threshold, being at the level of 0.1% for low energy thresholds (around 
1 GeV) but raising to order 1 per cent at, say, 10 GeV. 
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3'4.3. Initial-final state interference. While the initial- and final-state QED correc- 
tions are functions of s only, the initial- final state interference correction depends also on 
the fcrmion scattering angle and therefore induces a modification of the dependence 
of the Born differential cross section. It receives contributions from the interferences be- 
tween the real radiation diagrams of Figs. 22 and 23 and from the interference between the 
tree-level amplitudes and the QED box diagrams shown in Fig. 24. The leading angular 




Fig. 24. - The Feynman diagrams for 0{a) QED box corrections to s-channel process e'^e 



dependent terms that enter the result are of the form ln{t/u) = In (tan(i?/2)) [128,129]. 
This explains why the initial-final state interference is numerically small under many 
realistic experimental conditions, even if its actual magnitude crucially depends on the 
applied cuts. A physical argument can also be advocated to understand the role of the 
initial-final state interference in the region of the resonance. In fact, when the is 
produced close to its mass shell, the wave functions for initial- and final-state radiation 
are separated in space-time chie to the finite lifetime, so that their overlap is small. 
This implies that for loose cuts the interference effect is typically at the level of 10~^ or 
even smaller. Only when particularly tight cuts are imposed or one moves away from the 
peak, the above argument is invalidated and the interference can become more sizeable. 
In order to keep under control all the situations, the correction due to the initial-final 
state interference is included in the theoretical predictions, generally using the existing 
exact 0{a) calculations (including hard bremsstrahlung) , that are sufficient for the realis- 
tic ES. For very tight, unrealistic, photon energy cuts, a proper exponentiation procedure 
can be introduced. 

3'5. Computational Tools. - In order to perform actual precision tests of the elec- 
troweak theory in e+e~ collisions at the Z^ pole, very precise measurements of cross 

sections and asymmetries are required, that in turn imply, together with a high statis- 
tics, a deep knowledge of systematics effects, such as the acceptance, selection efficiency 
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and backgrounds for a given reaction. The tools used by LEP Collaborations to determine 
experimental acceptances and efficiencies are unavoidably Monte Carlo event generators 
since they allow to simulate experimental cuts with the maximum flexibility. Just to 
give some examples, the LEP experiments commonly use the generators HERWIG [130] 
and PITHYA/JETSET [131] to simulate the process e+e" hadrons, KORALZ [132] for 
the production of /i and t pairs, BHAGENE3 [133], BHWIDE [134] and UNIBAB [135] for the 
large-angle Bhabha scattering. HERWIG and PITHYA/JETSET are generators for the study 
of hadronic final states in e+e~, ep and pp collisions, describing the phase of parton 
cascade in the framework of perturbativc QCD and the hadronization mechanism by 
using independent theoretical models. KORALZ is a Monte Carlo program developed for 
s-channel e+e~ — > // processes, including electroweak loops and using YFS exponen- 
tiation (sec Appendix A) for the treatment of photonic radiation. Large-angle Bhabha 
generators will be described in the following. 

However, in addition to Monte Carlo Generators, fast analytical and semi-analytical 
programs are of utmost importance for precision physics, because, although they can 
account for simple kinematical cuts only, they are necessary tools in order to extract the 
electroweak parameters from the measurements according to an iterative fitting proce- 
dure. These fitting programs arc also said electroweak libraries. As a matter of fact, 
the fitting programs mostly used by the LEP Collaborations are BHM [136], MIZA [137], 
TOPAZO [111] and ZFITTER [138]. However, the extreme complexity in the calculation 
of radiative corrections and the relevance of improvements and cross-checks for preci- 
sion physics at LEP, motivated a number of theoretical groups to develop completely 
independent and original electroweak libraries. The availability of various independent 
electroweak libraries turned out to be particularly useful in the context of the job of the 
Electroweak Working Group at CERN [76]. The aim of such a collaboration work was 
the estimate of the theoretical error inherent to the SM predictions for e+e" — > // pro- 
cesses, in view of the final analysis of LEP precision data (sec Sect. 3'6 for more details). 
The basic features of the codes that contributed to reach the above task are shortly de- 
scribed in the following, with particular emphasis on their main physics input as well as 
their numerical output. However, for more details on technical aspects of these programs 
and their underlying formulation, the reader is referred to the original literature and to 
ref. [76]. 

BHM [136] — It is a semi-analytical program for the calculations of the parameters 
and realistic observables for an extrapolated set-up only. It relies upon the on-shell 
renormalization scheme for the formulation of weak loops and the QED radiator for the 
treatment of ISR (see Appendix A). 

LEPTOP [139] — It is an analytical code, developed by the ITEP Moscow group, giving 
predictions for the the parameters in the so-called a-Born approximation. It consists 

in taking as "tree-level approximation" the one obtained by using a{Mz) overall, and 
computing siia^ 'dw from a{Mz), and Mz, with the residual electroweak corrections 
on top of this. 

TOPAZO [111] It is a semi-analytical program, developed by the Pavia and Torino 
groups, that can be used to calculate both pseudo-observables and realistic observables, 
the latter over both an extrapolated and a more realistic set-up. It is additionally able 
to calculate the full Bhabha scattering cross section at large angles. It employes the MS 
scheme for the treatment of the electroweak corrections. QED corrections are exactly 
included at 0{a) for s-channel processes, at the LL level for t and s-t contributions to 
the large-angle Bhabha scattering. On top of that, higher-order QED corrections are 
implemented using the SF approach (see Appendix A). 
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WOH [140] — It is an analytical code for the calculation of the pseudo-observables only. It 
is based on the on-shell renormalization scheme and, basically, it leads back to the same 
approach of BHM, so that BHM/WOH arc not completely independent of one another. 
ZFITTER [138] — It is a semi-analytical program, of the Dubna-Zeuthen group, that, as 
TOPAZO, allows to obtain predictions for both the pseudo-observables and the realistic 
ones, but excluding large-angle Bhabha scattering. The on-shcll renormalization scheme 
is used for electroweak loops; QED corrections are exactly treated at 0{a), together with 
soft-photon exponentiation for ISR and FSR. 

To summarize, all the above codes can provide predictions for the Z'^ parameters 
and therefore can all be used as fitting tools to these de-convoluted quantities. BHM, 
TOPAZO and ZFITTER can also calculate realistic observables and fit data for such observ- 
ables. TOPAZO is also additionally able to calculate the full Bhabha scattering observables. 
Concerning the treatment of electroweak loops, the codes employ completely indepen- 
dent calculational schemes, with the exception of BHM/WOH that essentially are based on 
the same approach. QED corrections in the three QED dressers are treated according 
to different theoretical methods, whereas QCD corrections are common to a large extent 
in all the codes. How these codes were used as tools to estimate the intrinsic theoretical 
uncertainties in precision calculations for the Z*' resonance is explained in Sect. 3'6. 

The special role played by large-angle Bhabha (LABH) scattering in Z° precision 
physics led the to development of dedicated semi-analytical and Monte Carlo programs 
for such a process. Two of them, BHAGEN95 and TOPAZO, have been previously described, 
the former in Sect. 2'4 and the latter here in the present Section. The other programs 
used by the LEP experiments for the study of e+e~ — > e+e~ at large angles are briefly 
described in the following. It is worth pointing out that the Bhabha process, when con- 
sidering large scattering angles and the energy region around the Z° peak, is basically a 
s-channel resonant process with "small" non-resonant contributions, so that its dynam- 
ics is completely different from that of the same process at small scattering angles, that 
has been already discussed in details in Sect. 2. In particular, whereas non-QED effects 
other than vacuum polarization are absolutely negligible in the SABH case, an accurate 
treatment of the electroweak loops is a necessary ingredient for a precision calculation of 
the LABH reaction [78, 141, 142]. As a consequence, the computational tools (and their 
underlying formulations) for the LABH process are very different with respect to the 
ones used for the luminosity monitoring and already described in Sect. 2'4. 
ALIBABA [142] — It is a semi-analytical program, implementing exact 0{a) weak and 
QED corrections to all Bhabha channels. The higher-order QED effects consist of LL 
0{a^) plus soft-photon exponentiation. 

BHAGENE3 [133] — It is a Monte Carlo event generator, including one-loop and the most 
relevant two-loop electroweak corrections. The 0{a) QED corrections uses the exact 
matrix element and are supplemented with higher-order soft-photon effects. 
BHWIDE [134] — It is a recent Monte Carlo event generator, based on the electroweak 
library of ALIBABA for the treatment of 0{a) weak and QED corrections. It includes 
multiphoton radiation in the framework of Yennie-Frautschi-Suura (YFS) exponentia- 
tion (see Appendix A) and can be considered as the extension of the code BHLUMI (see 
Sect. 2'4) to large angles. 

UMIBAB [135] It is a Monte Carlo event generator, relying upon an updated version 
of the ALIBABA electroweak library. The QED corrections are implemented assuming 
s-channel dominance and using photon shower algorithms for ISR and FSR (see Ap- 
pendix A). 

A more extensive description of the LABH programs, together with global compar- 
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isons between the semi-analytical and Monte Carlo programs, can be found in the recent 
work of ref. [18]. At LEPl, when analyzing LABH data, the common procedure employed 
is the so-called t-channcl subtraction, where f. and s-t contributions arc subtracted from 
the data. As a matter of fact, the two programs commonly used to perform this unfold- 
ing are ALIBABA and TOPAZO. Furthermore, besides the programs described above, one 
should also mention the codes ALISTAR and MIZA, that arc basically rearrangements of 
ALIBABA developed for fitting purposes and specific experimental needs. 

3'6. Theoretical Uncertainties. - The very high experimental precision reached by the 
LEP Collaborations in the measurement of the parameters and realistic observables 
(of the order of 0.1% or even below) , necessarily requires a careful estimate of the intrinsic 
uncertainties associated to the SM theoretical predictions for these observables. In fact, 
as already discussed in Sect. 2'5 when addressing the problem of the total theoretical error 
in the luminosity measurement, any prediction obtained via a perturbative expansion is 
intrinsically affected by an uncertainty that is mainly due to neglecting higher-order 
contributions. The evaluation of such an uncertainty turns out to be a particularly 
severe problem whenever considering precision calculations for the Z'^ resonance, since 
the higher-order contributions depend in a highly non-trivial way on the whole stuff of 
radiative corrections in the SM (weak, QCD, mixed ew/QCD, QED). For this reason, 
a combined effort by different groups of theorists is in principle the best strategy to 
address this delicate subject. Such an effort was pursued by the collaboration work of 
the Electroweak Working Group, held at CERN, Geneva, during 1994 and concluded at 
the beginning of 1995. The composition of the Electroweak Working Group consisted of 
the authors of those formulations and relative computational programs containing at that 
time the state-of-the-art of the radiative corrections to the electroweak Z'^ observables: 
BHM, LEPTOP, TOPAZO, WOH and ZFITTER. They have been already described in some detail 
in the previous Section. Because all these codes basically include the same content in 
the perturbative expansion ( i.e. exact 0{a) electroweak corrections plus higher-order 
leading contributions), but differ in the choice of the renormalization scheme and in the 
treatment of higher-order sub- leading effects, they are ideal tools to provide an estimate 
of the theoretical uncertainties. The latter can be classified as follows: 

• Parametric uncertainties: are related to the experimental precision of the input pa- 
rameters. Among them, the largest uncertainty comes from a{Mz): \Aa~^{Mz)\ 
= 0.09; other sources are the errors on as{Mz) and on the heavy-quark masses. 
These uncertainties can be estimated comparing the effects of variations in the 
input parameters within their errors [76]. 

• Scheme- dependence uncertainties: are associated to the choice of the calculational 
scheme. In fact, in a given renormahzation framework the truncation of the per- 
turbative series is realized in some specific way. An estimate of these uncertainties 
can be obtained by comparing the predictions for a given observable obtained with 
two codes based on different renormalization schemes. This is possible since, as 
discussed in Sect . 3'5 , BHM/WOH and ZFITTER employ the on-sheU scheme, TOPAZO 
makes use of the MS scheme and an original approach is implemented in LEPTOP . 

• Intrinsic uncertainties: are inherent in a given renormalization framework. They 
are a consequence of the still missing higher-order terms, reflecting in a certain 

degree of arbitrariness on how to combine the various theoretical ingredients in 
order to get the final expressions for the observables. This arbitrariness can be 



PRECISION PHYSICS AT LEP 



59 



quantified using the concept of option, i. e. a set of alternative but equally plausible 
implementations of the building blocks of the radiative corrections within a given 
calculational scheme [143]. 

To better illustrate the meaning of intrinsic uncertainties, it is worth presenting 
and discussing a few examples of option. A first one [144] refers to different possible 
implementations of QCD and QED final-state corrections to a given partial width 
Tf — T{Z ff). Let us suppose that 0{a) QED and 0{al) QCD corrections to the 
above width are known but that the mixed 0{aas) corrections have not yet been calcu- 
lated. Therefore, one can combine the two above final-state corrections either according 
to a factorized representation 



(118) 



or an expanded one 



T^EW 
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In the latter equations, denotes the elcctroweakly-corrected partial width. Owing 

to the lack of knowledge of the exact 0{aas) corrections, the two realizations are of 
course equally correct and could be both implemented into a program as two different 
options. They differ by a term of the order of 3aasQj/47r^, that can be seen as a naive 
estimate of the unknown mixed corrections. The explicit calculation of such corrections 
should reduce the uncertainty moving it to 0{aa^). Indeed, this correction is today 
available, giving the result 
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SO that the two above realizations can not be longer seen as options. 

Another theoretical uncertainty, which is related to the pure weak corrections to 
the pseudo-observables, is the loading-remainder splitting in the effective couplings dis- 
cussed in Sect. 3'2.3. These generally contain a leading universal part, which is usually 
resummed, and a non leading (remainder) part. The way of performing the separation 
between leading and remainder part and the way of treating the last arc not uniquely 
defined, so that they are a source of theoretical uncertainty. A typical example is pro- 
vided by Ar (the same reasoning, as the one shown in the following, also holds for the 
quantities Spf and Skf) introduced in eq. (55) and split into a leading and a remainder 
part according to eq. (56), where the leading term contains the light fermion mass sin- 
gularities and the terms proportional to . In this separation the remainder contains, 
among others, logarithmic contributions in the top-quark and Higgs-boson masses, which 
can be numerically important. In resumming Ar as in eq. (61), to take into account of 
the leading terms to all orders, different ways of treating Arrem are in principle possi- 
ble, i.e. they are different options which give an estimate of the associated theoretical 
uncertainty [76]: 



1 - ArL - Arrem 

1 / ^^T'rf^'m 
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1 -\- ^Vrem 

1-ArL ' 
1 

1 - Ari ^ ^r,em. 

In the specific case of Ar, the last two expansions are not vahd according to the arguments 
given in ref. [68]. 

The two examples given above should clarify the meaning of intrinsic theoretical 
error and how to estimate it by considering proper options. Actually, apart from the 
factorization of QCD and EW corrections and the leading-remainder splitting, other 
options have been carefully examined by the Elcctrowcak Working Group, such as the 
choice of the scale of a in the non-leading corrections (in particular the vertex correc- 
tions), the linearization of the radiatively corrected quantities like sin^i^^^^, different 
choices of implementing the resummation of the vector-bosons self-cnergics, in particular 
the terms related to the Higgs-boson contribution, and the choice of the scale of Ug in 
the mixed electroweak-QCD corrections. The results obtained by the work of the Elec- 
troweak Working Group, based on a careful analysis of both the pseudo-observablcs and 
realistic ones, allowed to draw at that time (beginning of 1995) the following conclusions 
on the theoretical accuracy of the precision calculations for the resonance [76] : 

• the differences between the predictions of different computational tools are small 
compared to the experimental errors; 

• the parametric uncertainty due to Auhad is the dominant source of error and the 

real "bottleneck" to improve the theoretical accuracy; only new accurate measure- 
ments of the cross section of e+e^ ^ hadrons at low energy could reduce it; 

• complete one-loop calculations, supplemented with highcr-or(l(T leading effects are 
adequate for Z° precision physics, but the control on two-loop electroweak sub- 
leading corrections would significantly reduce the still remaining uncertainty. 

The first two conclusions still remain valid today. Concerning the third point, 

progress in the calculation of important two- loop electroweak sub-leading corrections has 
been achieved after the conclusion of the work of the Electroweak Working Group. In 
fact, the two-loop next-to-leading heavy top-quark contributions of the order of Gi^rn^M^ 
to the Mw-Mz interdependence and to sin^ ^1// ^^^^ been completely calculated in the 
MS framework and for two different realizations of the on-shcU scheme [88, 89] (see 
also [90,91]). These newly calculated corrections give a nice reduction of the theoretical 
error of electroweak origin, especially for Mw and sin^i^^^j [89, 145], even if it should 
be noticed that they have not yet been implemented, at the time of writing, in standard 
computational tools for LEPl precision physics. (^) 

Concerning the theoretical uncertainty associated to the treatment of QED correc- 
tions, of utmost importance for the study of the realistic distributions, the conclusions 
of the analysis performed by the Electroweak Working Group, as obtained through a 
critical comparison of QED dressers over both an extrapolated and realistic set-up, can 
be summarized as follows: 



(•^) After the completion of this work, the two- loop electroweak sub-leading corrections of refs. [88,89], 
together with the QCD calculation of ref. [106] and the QED eflfects studied in ref. [120], have been 
implemented in a new version of the code TOPAZO [146]. 
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Fig. 25. - A comparison between theoretical predictions for the hadronic cross section. The 
squares, diamonds and crosses represent BHM [136], TOPAZO [111] and ZFITTER [138], respectively 
(from ref. [76]). The lower window shows the predictions of the codes together with their 
estimate of the theoretical error. 



• for s-channel annihilation processes, the theoretical error is of the order of 0.1%, 
almost independently of the considered cm. energy (see Figs. 25 and 26); 

• for full Bhabha scattering, the theoretical accuracy can be estimated, with due 
caution, to be at the level of 0.1-0.2% before and strictly around the peak, growing 
to about 1% (depending on the imposed experimental cuts) on the hard radiative 
tail. 

It is worth pointing out that the above uncertainties associated to the photonic 

corrections actually match the statistics and systematica of the LEPl data for the realistic 
distributions. In fact, whenever a somehow large theoretical error is present for full 
Bhabha scattering, this shows up in an energy regime (the hard radiative tail above the 
peak) where the experimental precision is much lower (i.e. around an order of magnitude 
larger) than the corresponding theoretical accuracy. 
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Fig. 26. - The same as in Fig. 25 for the fj, forward-backward asymmetry (from ref. [76]). 



As far as s-channel processes are considered, the estimate of the QED theoretical 
accuracy provided by the Electroweak Working Group has been recently reinforced by 
a novel investigation of the effects of LL 0{(3'^) photonic corrections to the QED radia- 
tor [120], usually neglected in standard computational tools and not taken into account 
in the analysis of the Electroweak Working Group. In fact, it has been shown that the 
0{f3^) corrections introduce, on an extrapolated cross section computed with a standard 
QED radiator with up to 0{l3'^) non-soft terms, a systematic shift of around -0.07%, 
that confirm a posteriori the estimate of the Electroweak Working Group. However, the 
size of the effect also indicates that these LL third-order QED corrections should be 
carefully takc'n into account in the theoretical predictions for the realistic obscrvablos of 
two-fermion production, in the light of the latest experimental data for some realistic 
observables, such as the cross section of e+e~ — > hadrons. 

Concerning the estimate of the theoretical error in the large-angle Bhabha scatter- 
ing, it substantially traces back to a certain disagreement between the predictions of 
ALIBABA and TOPAZO far from the peak and due, as discussed in the literature [78], to 
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a different implcnicntation of QED final-state corrections, yielding a difference in the 
treatment of 0{a^) sub-leading contributions. In fact, in ALIBABA final-state corrections 
arc implemented factorizing the leading terms only and summing up at 0{a) non- leading 
effects, whereas in TOPAZO a fully factorized prescription is followed, in order to preserve 
the so-called classical limit according to which the cross section of a real scattering pro- 
cess must vanish in the absence of electromagnetic radiation. This fundamental property 
is not guaranteed by an additive formulation. Furthermore, as shown and extensively 
discussed in [23] , a fully factorized form allows to keep under control in an effective way 
the bulk of 0(0^) next-to- leading QED corrections of the order of 0{a^L), that are on 
the contrary completely absent in the additive prescription (see also the discussion given 
in Sects. 2"3, 2'5 and 2"6 on the role of the O(a^L) corrections in the present theoretical 
error to the SABH process). Therefore, with respect to an additive formulation, a fac- 
torized formulation is more accurate from the point of view of the perturbative content 
as well as more theoretically founded. However, after the completion of the work per- 
formed in ref. [76] , a much more extensive analysis of the large-angle Bhabha process has 
been performed in ref. [18]. In this last analysis, several codes other than ALIBABA and 
TOPAZO, noticeably BHWIDE, have been used for the comparisons. Moreover, also more 
realistic ES's have been adopted, namely calorimetric ES's in addition to the unrealistic 
BARE ES adopted in ref. [76] , and this is a key point in order to obtain a reliable error 
estimate (see the discussion in Sect. 2 for more details). Thanks to these results, it may 
be concluded that the theoretical error of the large-angle Bhabha scattering is 0.3% on 
peak, and 0.5% off peak [18]. 

4. — Fits to Precision Data 

The application of the theoretical results discussed in Sect. 3 to the analysis of preci- 
sion data collected at LEPl/SLC leads to the indirect determination of the fundamental 
parameters of the Standard Model (SM). Actually, the high experimental accuracy makes 
the electroweak measurements on the peak sensitive to the mass of the particles cir- 
culating in the loops, although they arc not energetically accessible. Hence the precision 
data can be used to infer valuable constraints on the top-quark and the Higgs-boson 
masses, as well as on the value of the strong coupling constant ag- The level of agree- 
ment between theory and experiment allows in addition to derive hints on possible new 
physics scenarios beyond the standard description of fundamental particles and interac- 
tions. 

The present section is devoted to illustrate the procedure adopted by the experi- 
ments to extract the electroweak Z° parameters from the directly measured production 
cross sections and forward-backward asymmetries. The interpretation of the experimen- 
tal results in terms of the SM parameters is discussed. The most recent indirect limits 
obtained for the top-quark and Higgs-boson masses via the virtual effects due to radia- 
tive corrections to the precision observables are presented, and compared with present 
information from direct searches. The determination of from precision data at the Z'^ 
pole is compared with other measurements in different processes and at different energy 
scales. The issue of a more model-independent attitude towards the precision measure- 
ments is in conclusion addressed, discussing the S'-matrix approach as a framework for 
a model-independent determination of the Z-boson parameters, and the e parameteriza- 
tion of pure weak loops as a strategy for a model-independent analysis of the electroweak 
quantities. The implications of precision data for new physics beyond the SM are briefly 
examined. 
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41. Determination of the Electroweak Parameters. - The "primary" observables 
measured by the LEP experiments are the cross sections and forward-backward asymme- 
tries of the two-fermion processes e+e^ ^/Z f f quoted as a fimction of the centre 
of mass (cm.) energy (see Sect. 3'1). The channels considered are the hadronic (qq) and 
leptonic I = e,iJ,,T) final states. The cross section data are determined by making 

use of the relation of eq. (2) in Sect. 2'1, once the number of events for a given channel 
has been determined and corrected for the trigger efficiency, the geometrical acceptance 
and the eflaciency of the selection cuts. Analogously, the data for the forward-backward 
(FB) asymmetry ApB are obtained exploiting the definition 

where Np(Nb) is the number of events collected with a forward (backward) scattered 
fermion. This determination of the FB asymmetry is known as the counting method. 
Alternatively, the FB asymmetry is obtained from a max;imum likelihood fit of an asym- 
metric differential cross section of the form 

dcT 8 

(122) (X l + cos^'&+ -ApBCOs'd 

a cos 17 3 

to the experimental angular distribution. The expression given by eq. (122), where •& is 
the angle between the incoming electron and the outgoing fermion, is the form predicted 
by the electroweak theory in the lowest-order approximation and applies to s-channel 
annihilation only. For e+e~ — > e+e~, due to the presence of f-channel contributions, a 
similar but more complicated expression is used. This second determination of the FB 
asymmetry is known as the fitting or likelihood method. 

From the measured hadronic and leptonic cross sections and leptonic forward-back- 
ward asymmetries, the parameters of the Z'^ resonance are extracted by the LEP Collab- 
orations by means of a combined fit to the data [147]. The Z° parameters are obtained 
using a minimization procedure with the correlations among the data (common ex- 
perimental errors, theoretical luminosity error, uncertainties from the LEP energy cali- 
bration, etc.) taken into account using a covariance matrix. Several sets of parameters 
can be (and actually are) introduced to parameterize the measurements. However, in 
order to make possible a combination of the results obtained by each of the four LEP 
Collaborations as well as for comparison purposes, the LEP experiments use a standard 
set containing nine free Z^ parameters: 

• the Z-boson mass, Mz, and its total width Tz', 

• the hadronic pole cross section cr^, defined by eq. (31) in Sect. 3T; 

• the ratios Re,fj.,r = ^h/^e,fi,T; 

• the pole asymmetries Ap^, Ap^, Ap^, that can be expressed in terms of the effec- 
tive vector and axial-vector neutral current couplings of the fermions according to 
the relations given by eqs. (32) and (35) in Sect. 3'1. 

These parameters are chosen because they are most directly related to the experimental 
quantities and are weakly correlated. If the assumption of lepton universality is addi- 
tionally introduced, then the numbers of parameters reduces from nine to five, i.e. 

(123) Mz, Tz, al Ri=Th/Ti, A^^, 
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that actually is another standard set of parameters commonly used by the LEP Collab- 
orations. 

The extraction of the Z'^ parameters from a combined fit to the hadronic and leptonic 
cross sections and leptonic forward-backward asymmetries requires that an appropriate 
parameterization for these last realistic observables is introduced. By exploiting the 
relation between the effective couplings and the partial width, the cross section can 
be written in terms of mass, total and partial widths of the Z boson, without relying upon 
any particular assumption on the validity of the Standard Model (SM) [113, 148, 149]: 



where 



(125) I^-(^)=^Ha?|+^^ 

and the partial widths are understood expressed in terms of the Z° parameters. The 
term identifies the 7-exchange contribution known from QED, Rf and If describe 
the "f-Z interference and e is the correction due to finite final-state fermion mass effects. 
This procedure justifies the typical Breit-Wigner ansatz for the Z° contribution adopted 
by the LEP Collaborations, 

f^r,^^ f _ p f 

^ ' "^^^ M| ^^(s-M|)2 + s2r|/M|' 

for a given // channel, although the forms actually used by each experiment may dif- 
fer in some detail. The photon-exchange contribution is calculated from QED and the 
j-Z electroweak interference is taken into account as well, usually computing it within 
the SM. Although the interference contribution is small around the Z° peak, it is im- 
portant for the precise measurement of Mz and therefore some more general fit, where 
the SM dependence of the jZ effect is somehow relaxed, are actually performed by the 
experiments (see Sect. 4'3.1 below). 

It is worth pointing out that the above parameterization is introduced in order to 
describe just the short-distance behaviour of the cross section. Indeed, in the fit of 
the parameters to the measured data, the QED radiative corrections enter as an 
unavoidable theoretical ingredient, and they are taken into account according to the 
algorithms described in Sect. 3'4 and in Appendix A. A complication in the fitting 
procedure arises when considering the Bhabha channel e~^e~ e+e~, because of the 
presence of t-channel contributions. To minimize this contamination the experiments 
adopt a selection with large polar angle acceptance (typically 45° < < 135°) where the 
,s-channel Z-boson exchange largely dominates. Th(^ remaining t and s-t contributions 
arc generally calculated by using the programs ALIBABA and TOPAZO, in order to correct 
the cross section. The most recent LEP results for the Z-boson mass and width are 
shown in Figs. 27 and 28, respectively. 

In addition to the above parameters determined by means of a combined nine 
(five) parameters fit to the cross sections and asymmetries, other important electroweak 
observables are measiircd by the LEP experiments. Each of them is obtained according to 
a well-defined experimental strategy applied to a specific // final-state. A first example 
is given by the r polarization asymmetry P'^ that is defined as [150] 

(127) = 



s-Ml 



z 
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Fig. 27. - The Z-boson mass according to the most recent data analyses [7,8]. 



where <Jfi(^L) denotes the r-pair cross section for the production of a right-handed (left- 
handed) T~ . This quantity is experimentally determined by measuring the longitudinal 
polarization of r pairs produced in decays. The left-right asymmetry Alj^ is another 
important manifestation of parity violation in weak interactions, that, however, requires 
longitudinal polarization of the initial-state electrons. Owing to the difficulties in achiev- 
ing a high degree of polarization at LEP collider because of the presence of depolarizing 
effects, the left-right asymmetry has not been measured by the LEP experiments but 
by the SLD Collaboration at SLC, where the conditions for polarization measurements 
are much more favourable. In fact, an average degree of polarization of the incoming 
electrons < > reaching 77% has been obtained at SLC. The left-right asymmetry is 
determined by the SLD experiment exploiting the relation [151] 



(128) 



1 iVe, - N, 



<Pe> Ne, + iVe, 



where A^e; and A^e^ are the observed numbers of hadronic events using left-handed and 
right-handed electron beams, respectively. 

Further interesting electroweak observables in decays are derived from the anal- 
ysis of the final states containing c and b quarks {heavy flavours). These measurements 
require sophisticated tagging techniques of heavy flavours (such as mass or lifetime tag- 
ging), whose description is beyond the aim of the present review (see for instance ref. [152] 
for an overview of the tagging methods used at LEP). The Z'^ parameters that are de- 
termined from these studies are 



(129) Rt, = n/Th, Rc = rjTh, A'ps, A 



FB^ 



where F^fc) is the Z^ width for the decay into 6(c) quarks, A''pg,Apg are the b- and 
c-quark pole forward-backward asymmetry. As already remarked in Sect. 3T, the results 
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Fig. 28. - The Z-boson width according to the most recent data analyses [7,8]. 



for A''pg and Apg are quoted after having subtracted from the data, besides the effect 
of initial-state radiation (ISR), also the effects of QCD corrections, in order to deal with 
a pure electroweak observable. 

Once the above "primary" Z° parameters are extracted from the data, other impor- 
tant additional quantities, such as the partial widths F/i, F; and Tinv, are derived. 
^inv denotes the Z-boson invisible width that can be obtained via the relation (assuming 
lepton universality) 

(130) ^,r^v = Tz-Th-STi. 
Its present experimental value is [7, 8] 

(131) Tiny = 500.1 ± 1.8 MeV. 

This quantity is of particular interest for the extraction of the number of light neutrino 
species Ni, . In order to get rid of the bulk of the non-negligible top-quark mass depen- 
dence of the partial widths, iVj, can be conveniently derived from the comparison of the 
experimental ratio Tinv/Ti with the corresponding SM prediction, yielding the result [7,8] 



(132) 



= 2.993 ±0.011. 
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This determination, that firmly proves the existence of three standard hght neutrino 
famiUes, agrees well with the independent measurement obtained by means of the process 
e^e~ known as radiative neutrino counting reaction. It gives rise to a signature 

where only a single bremsstrahlung photon and nothing else is seen in the detector. The 
Feynman diagrams contributing in the SM to this process are depicted in Fig. 29. 




Fig. 29. - The tree-level Feynman diagrams for the process e^e — > vv^. 

The study of the radiative neutrino counting events offers the possibility to deter- 
mine in a rather direct way the partial width of the into invisible particles, since the 
.Z-boson exchange contribution, largely dominating the e"'"e~ vv^f cross section around 
the Z'^ peak, is proportional to the number of neutrinos. For that reason, this procedure 
is known as direct determination of the invisible width. In practice, the measured cross 
section as a function of the cm. energy is compared with the SM calculation for dif- 
ferent numbers of neutrino generations. A meaningful comparison between theory and 
experiments requires that electroweak radiative corrections, especially the large effects 
introduced by ISR, to e+e" vv^f are taken into accoimt, although the statistical preci- 
sion of the direct determination of the invisible width is smaller than the one achieved in 
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the measurements of the observables of two-fermion production processes (see ref. [153] 
and references therein). Yet, the study of radiative neutrino counting events clearly rules 
out the existence of a fourth family with light neutrinos in the SM. 

In addition to the number of light neutrinos, there are other derived electroweak pa- 
rameters that deserve mention for their role in precision tests of the electroweak theory. 
From the partial widths of the into charged Icptons, b and c quarks, the r polariza- 
tion asymmetry, the left-right asymmetry as well as the forward-backward asymmetry 
in leptonic and heavy- flavour channels, the LEP/SLC experiments determine the vector 
and axial-vector neutral current couplings of fcrmions. These results for the effective Z- 
boson couplings Qy and make use of the relations between the Z° parameters and the 
couplings themselves given by eqs. (32)-(34) and (35) for the asymmetries and polariza- 
tions, and eq. (36) for the partial widths given in Sect. 3'1. In particular, the comparison 
between electron, muon and tau couplings is in good agreement with lepton universality, 
with a precision of about 0.2% for gA. and 5 -t- 10% for gv [8]. Furthermore, since the 
forward-backward and polarization asymmetries depend only on the ratio gy/g^, it is 
possible to express the asymmetry measurements in terms of a single parameter given 
by the effective weak mixing angle, i.e. 

J 

(133) A\Qf\sir?^{ff = l-^-. 

9a 

The value of sin^ t^^^^ is flavor-dependent due to the effects of weak vertex corrections 
that are non-universal. However, since the most precise results are obtained for the 
charged leptons, the mostly quoted value for the effective weak mixing angle coincides 
with sin^ ^1//' ^^^^ allows to combine into a single electroweak parameter all the leptonic 
asymmetry measurements. 

It is worth pointing out in conclusion that, as a result of the; procedure described 
in the present Section, the Z° parameters that are extracted by the LEP experiments 
carry dependence on the whole stuff of electro-weak, mixed electroweak/ Q CD as well 
as, for the quantities referring to qq final states, QCD radiative corrections, reviewed 
in Sects. 3'2-3'3. Therefore, these observables can be conveniently used to measure 
quantum effects of the electroweak interactions, i.e. to derive interesting constraints on 
the SM parameters that enter the predictions via radiative corrections only, noticeably 
the top-quark mass rrit, the Higgs-boson mass tuh and the value of the strong coupling 
constant at the pole as{Mz)- Such an analysis assumes the validity of the SM ab 
initio. However, if an appropriate strategy is organized, a model-independent analysis 
of the Z° parameters can be also performed, exploiting the high level of precision of the 
data in order to constrain new physics predicted by possible extensions of the SM. The 
physical and theoretical implications of the precision measurements of the electroweak 
Z^ parameters are discussed in more detail in the following. 

4'2. Standard Model Fits. A summary of the most recent measurements of the Z'^ 
parameters obtained by the LEP and SLC experiments is given in Tab. VI. According to 
a standard presentation, the data are compared with the SM predictions (third column) 
corresponding to a fit of the electroweak data in terms of rrii, mu and ag{Mz), that 
will be described in more detail in the following. In the fourth column one can also see 
the pulls derived from the fit, i.e. the difference between each theoretical result and the 
corresponding experimental measurement, in units of the measurement error. A few but 
important comments about Tab. VI are in order here. 
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The first remark concerns the precision level reached in the experimental measure- 
ments at LEP and SLC. The Z-boson mass Mz is now known with a relative error of 
2 X 10""'' (absolute error of ±2 McV). The total Z° width and the hadronic peak cross 
section are determined with a relative experimental uncertainty of the order of 1 x 10~^ 
(with an absolute error for Tz of around ±3 MeV). All the asymmetry measurements 
have an absolute error smaller than 0.01, generally at the level of a few 0.001, and the 
effective weak mixing angle is derived with a precision of around ±0.0003. Therefore, in 
the light of these numbers, the LEPl/SLC experimental program can be considered as 
a great success, being the precision reached better than that expected at the beginning 
of operation. 



Table VI. - Electroweak data according to the most recent determination (see refs. [7] and [8]). 



Quantity 


Data (Jerusalem '97) 


Standard Model 


Pull 


Mz [GeV] 


91.1867(20) 


91.1866 


+0.04 


Tz [GeV] 


2.4948(25) 


2.4966 


-0.73 


al [nb] 


41.486(53) 


41.467 


+0.36 


Ri 


20.775(27) 


20.756 


+0.71 


Rb (LEP+SLC) 


0.2170(9) 


0.2158 


+1.38 


Rc (LEP+SLC) 


0.1734(48) 


0.1722 


+0.24 


-^FB 


0.0171(10) 


0.0162 


+0.89 


Ar 


0.1411(64) 


0.1471 


-0.93 




0.1399(73) 


0.1471 


-0.98 


A\b 


0.0984(24) 


0.1031 


-1.95 




0.0741(48) 


0.0737 


+0.09 


Ab (SLC direct) 


0.900(50) 


0.935 


-0.69 


Ac (SLC direct) 

sin^ Qej f (LEP-combined) 

ALR^sm^e^ff (SLC) 


0.650(58) 


0.668 


-0.31 


0.23199(28) 


0.23152 


+1.68 


0.23055(41) 


0.23152 


-2.37 


Mw [GcV] (CDF/D0+LEP2) 


80.43(8) 


80.375 


+0.69 




0.2254(37) 


0.2231 


+0.63 


mt [GeV] (CDF/DO) 


175.6(5.5) 


173.1 


+0.45 



A second comment regards the level of agreement between theory and experiment. 

No significant evidence for departures from the SM predictions is present. Close to the 
end of the LEPl/SLC experimental program, the data support the SM in a remarkable 
way. Looking at the data in more detail, one could point out that there is an around 
2(7 deviation of the measured value of with respect to the SM prediction, that 

Rh is 1.4(7 away from the SM expectation, that the determination of the effective weak 
mixing angle from A^ji at SLC differs significantly from the LEP average. However, 
one should also emphasize that the observed pulls just follow the pattern that can be 
expected from a normal distribution of experimental measurements [8,154]. Moreover, 
the disagreement between the measured values of Rb and Rc and their SM expectation, 
that received much attention in the past (sec for instance ref. [155] and references therein), 
has now substantially disappeared. Presumably, future improved measurements of the 
parameters, especially in the 6-quark sector at LEP and Alr at SLC, will clarify the 
above issues. 

Last but not least, the electroweak data contained in Tab. VI, by virtue of their 
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Fig. 30. Data-theory comparison for realistic cross section data. Experimental data from 
ref. [156]. Theoretical predictions by TOPAZO [111]. 



high precision, can be used to derive constraints on the top-quark and the Higgs-boson 
masses, as well as on the value of as{AIz), fully exploiting the predictive power of the 
SM, as a renormalizable quantum field theory, beyond the tree-level approximation (for 
instance, several fits performed by various authors, using older data sets, can be found in 
refs. [157 165]). Hence, the sensitivity of the electroweak data to quantum loops elevates 
the precision physics to the level of discovery physics. Moreover, the analysis can be 
generalized in order to constrain, disfavor or even rule out various extensions of the SM. 




Fig. 31. - Data-theory comparison for realistic forward-backward asymmetry data. Experimen- 
tal data from ref. [156]. Theoretical predictions by TOPAZO [111]. 



The procedure described above is the one that is commonly adopted by the LEP 
experiments. Its main advantage consists in the fact that it allows the extraction of 
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parameters (the Z° parameters) independent of the experimental details and hence easily 
comparable from experiment to experiment, both at different machines and energy scales. 
Moreover, and as a corollary, the procedure offers the possibility of fitting the unknown 
SM parameters by making use of a theoretical machinery which is completely independent 
of the experimental selection criteria. On the other hand, it is intrinsically a two-step 
procedure, requiring firstly the determination of the Z'^ parameters and secondly the 
fitting of the unknown SM parameters. It is worth noticing that an alternative but 
not antithetic procedure is possible, namely considering as fundamental experimental 
quantities the data for the realistic observables (sec Sect. 3'1), and fitting the unknown 
SM parameters directly on them. This last procedure requires the use of theoretical tools 
that are able to provide predictions for the realistic observables, taking into account 
all the details of the experimental set-up. As a consequence, the required theoretical 
framework is sensibly more involved, but on the other hand the procedure is a single- 
step one. This last strategy has been used in the past by the LEP Collaborations and 
also by other authors [78, 166], yielding results compatible with the outcomes of the 
standard procedure. Figures 30 and 31 show the comparison between the experimental 
data for realistic cross sections and forward-backward asymmetries of ref. [156] and the 
theoretical prediction at best fit performed by TOPAZO [111]. 

4'2.1. Determination of the iop-quark mass. One the most important recent achieve- 
ments in particle physics has been the discovery of the toy quark by the experiments 
CDF and DO at the vv colhder TEVATRON at FERMILAB in Chicago [167]. The 
present average of the mass values reported by the CDF and DO Collaborations is 
rrit = 175.6 ± 5.5 GeV [168], with a measured top-quark production cross section in 
fairly good agreement with the QCD prediction [169]. The main production mechanism 
for the to'p quark at the TEVATRON is the quark- antiquark annihilation into a gluon 
followed by the creation of a tt pair. Due to its short lifetime, the top quark can not 
be detected directly but only via its decay products, i.e. the weak decay t — > fol- 
lowed by the subsequent leptonic or hadronic VK-boson decays W li>i,ud. Hence, the 
production of the top quark at the hadron colliders gives rise to six-fermion final states 
{bbudWi,bbuddu), as depicted in Fig. 32. 




Fig. 32. - The main production mechanism for the top quark at the TEVATRON, yielding a 
six-fermion final state. 

The top-quark mass is reconstructed by CDF and DO measuring the energy and 

momenta of the final-state products in the dilepton (where both W decay leptonically) 
and lepton-plus-jets (where one W decays leptonically and the other one hadronically) 
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events. (^) The large value obtained for nit (the top quark is as heavy as a heavy nucleus) 
has as a direct consequence that the total width of this particle depends on the third 
power of the mass, reaching the value of about 1.5 GcV for rrit ~ 175 GcV. As said 
above, this corresponds to an extremely short lifetime of about 10~^^ s, so that the top 
quark decays immediately after its formation and tt bound states cannot be formed. 

The direct evidence for the top quark in proton collisions and the measurement of 
its mass in a range around 170-180 GeV, in agreement with the indirect determination 
performed at LEP/SLC, constitutes a striking success for the SM at the quantum level 
and a remarkable confirmation of the precision tests of the electrowcak interaction at 
LEP/SLC. Indeed, although the LEP energy is not sufficient to produce real top quarks, 
the effects of this particle and, in particular, of its high mass, can be felt in the precision 
measurements, via the radiative corrections induced by the loops containing the top 
quark as a virtual particle. In particular, thanks to the high statistics collected at the 
resonance and the leading quadratic dependence of the theoretical predictions on mt, 
the mass value of the top quark can be deduced from the virtual effects in elcctroweak 
observables with a precision comparable to the accuracy of the direct measurement. (^) 

Since the mass of the top quark is at present known from CDF and DO with a rather 
good precision, it is reasonable to adopt two different strategies in fitting the top-quark 
mass from the precision data, namely 

• consider all the available electroweak data with the exclusion of the CDF/DO mea- 
surement of rrit; 

• consider the whole set of data, including mt from CDF/DO. 

The former procedure allows to establish whether the estimate of mt from radiative 
corrections is in agreement with the direct measurement at the TEVATRON. The latter, 
being a more general type of fit, allows to test the overall consistency of the SM in a 
more complete way and to obtain the present best estimate of derived quantities such as 
the W-boson mass and sin^ ^e//- 

The results obtained according to the above strategies, when fitting the data with 
mt, mn and as{Mz) as free parameters, are summarized in Tab. VII. As can be 



Table VII. - Fits to as, mt and tuh (from ref. [8]). 



Parameter 


LEP (incl. Mw) 


All but Mw and mt 


All data 


mt [GeV] 
mn [GoV] 
a,(Mz) 
xVdof 


158+lt 

oq+168 
00_4<j 

0.121 ±0.003 
8/9 


157+1° 

4lt«t 

0.120 ±0.003 
14/12 


173.1 ±5.4 

0.120 ±0.003 
17/15 



seen from the table, the fitted mass of the top quark is somewhat lower than the direct 
measurement, even if in well agreement with it. This compatibility, as already remarked, 
is an impressive confirmation of the SM at the level of quantum effects and justifies 
combining the direct and indirect derivation of mt- The la errors returned by the fit 



(■*) Recently CDF employed also the fully hadronic sample, that however gives a less precise determina- 
tion of the fop-quark mass. 

(^) The values quoted for mt refer to the pole mass. 



74 



G. MONTAGNA, O. NICROSINI and F. PICCININI 



PRELIMINARY 



> 



g-180 



1 +2ct ZFITTER 
1CT BHM 
lo-TOPAZO 



direct 
measurements 




indirect 
measurements 
(LEP + SLC + i/N) 
within min. SM 
I 



[ GeV ] 



Fig. 33. - Contour plot showing the top-quark mass mt versus the VK-boson mass. The 
shaded area shows the SM prediction for the Higgs-boson mass {Muigga in the figure) between 
60 GeV and 1 TeV. The theoretical predictions axe performed by means of the codes BHM [136], 
TOPAZO [111] and ZFITTER [138] (see ref. [7]). 



include the uncertainties due to the error on a{Mz) (inducing a mt variation of about 
±4 GcV) and to missing higher-order electroweak corrections (affecting mt by about 
±1 GeV). 

By adopting the procedure of including the direct measurements of mt among the 
input data, then typical results obtained are [7, 8] 

mt = 173.1 ±5.4 GeV , 

mH = 115166^ (or '^H < 420 GeV at 95% CL), 

(134) as{Mz) = 0.120 ± 0.003. 

From this fit, an indirect measurement of the Vl^-boson mass Mw can be derived (see 
Fig. 33), i.e. 

(135) Mw = 80.375 ± 0.030 GeV. 



This small error of 30 MeV on Mw from radiative corrections is clearly a challenge for 
future direct precision measurements of Mw at LEP2 and at the TEVATRON. 
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4'2.2. Determination of the Higgs-boson mass. As already shown in Sect. 4'2.1, the 
electroweak precision data can be used to infer indirect Hmits on the mass of the yet 
elusive Higgs boson. As discussed in Sect. 3'2, the leading effect of the Higgs-boson 
mass on the observables is only logarithmic and correlated with rrit. These difficulties, 
associated with other problematic aspects that will be discussed later, naturally set the 
question how reliably the Higgs-boson mass mn can be predicted from electroweak data. 
Clearly, the answer to this question is of utmost importance for planning the search for 
the Higgs boson at LEP2 and future accelerators. 




The indirect bounds that can be put on mn from radiative corrections effects and 
their validity are matter of debate in the community of particle physicists and an unam- 
biguous, unique answer to the above question can not be given, especially whenever one 
takes into account the large amount of investigations and detailed analyses present in the 
literature, many of which appeared after the discovery of the top quark at the TEVA- 
TRON [170-181]. Nonetheless, there are aspects concerning the constraints on ttih that 
are common to most of the analyses of precision data and therefore deserve special men- 
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tion. Before entering the details of the discussion, it is worth recalhng that negative 
searches at LEP put the present lower limit on the Higgs-boson mass mjj > 77 GeV at 
95% CL, derived assuming the validity of the minimal SM [182, 183]. 

The first common feature concerns the dependence of the fitted central value for rriH 
and of the la errors on the set of input data considered. This rather strong dependence 
can be clearly seen, for example, from Tab. VII, looking at the variations of several 
tens of GeV on niH, when considering the LEP data alone, all the available the data 
with the exception of Mw and rrit, and the whole set of data. It has been in particular 
emphasized by several aiithors that the exclusion of the SLD data for Alr moves up rriH 
considerably, as a symptom of the clash between the sin^ 1?^^.^ determination from LEP 
and SLC. Actually, the SLD data alone leads to an unnatural bound on mjj at the la- 
level, in conflict with the limit from direct searches at LEP [161, 170, 184]. 




nriH [GeV] 

Fig. 35. - Ax^ versus the Higgs-boson mass, taking into account the parametric uncertainty on 
a{Mz) (from [185]). 



A second common feature concerns the effect of including in the fit the theoretical 
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uncertainties, both intrinsic and parametric, of the SM calculations discussed in Sect. 3'6. 
Contrary to rut, that is marginally influenced by missing higher-order electroweak cor- 
rections and parametric uncertainties, the upper boimds on mn arc sensibly affected 
by parametric and intrinsic theoretical errors [179]. For example, the upper bound at 
95% CL can vary by around 100 GeV as a result of different implementations of radiative 
corrections beyond the one-loop level, so that the ^ fimction of mn appears as a 

band, representing the intrinsic theoretical error of the SM predictions, rather than a 
single curve (see Fig. 34). (^) Similar considerations apply to the distribution as a 
function of mn when the input parameter of the fit a{Mz) changes within its error. 
This parametric variation is actually larger than the one induced by different treatments 
of higher-order effects, pointing out the limitation imposed by the uncertainty on a{Mz) 
in a sensible derivation of mjj from precision data (see Figs. 35 and 37). 



2 2.5 3 




Fig. 36. - Contour plot showing the top-quark mass {Mtop in the figure) versus the Higgs-boson 
mass {Mniggs in the figure). The theoretical predictions are performed by means of the codes 
BHM [136], TOPAZO [111] and ZFITTER [138] (see ref. [7]). 

In the light of the above caveat in establishing a quite precise upper bound on tuh, 



(®) The recently calculated two-loop next-to-leading electroweak corrections of the order of G'^^M^m^ are 
expected to reduce this theoretical uncertainty, provided tlicy arc included in the standard electroweak 
libraries. Actually, after the completion of this work, the program TOPAZO has been upgraded to include 
these corrections, together with the results of refs. [106] and [120] concerning QCD and QED corrections, 
respectively [146]. 
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it turns out to be difficult going beyond the conclusion that the electroweak precision 
data, even when accounting for the direct measurement of m*, imply an indirect upper 
bound at 95% CL of around 400-500 GeV, with a preference for a "light" Higgs boson of 
mass around 100 150 GeV (see Figs. 36 and 37). 

4'2.3. Determination of a^. Besides the mass of the top quark and of the Higgs 
particle, the precision measurements of the electroweak parameters at the pole provide 
the opportunity of an accurate determination of the coupling constant of the strong 
interactions a^. Actually, e+e" machines are an ideal laboratory for QCD studies. Being 
the hadronic activity in e+e~ collisions restricted to the final state, the measurement of 
as, as well as further tests of QCD, can be carried out in a particularly clean environment, 
where the experimental signatures of hadronic events are largely free of backgrounds (see 
ref. [186] for reviews on the subject). 




Fig. 37. - Contour plot showing as(Mz) versus the Higgs-boson mass {M Higgs in the figure). 
The theoretical predictions are performed by means of the codes BHM [136], TOPAZO [111] and 
ZFITTER [138] (see ref. [7]). 

Owing to the renormalization group dependence of the strong coupling constant 

on the hard energy scale Q (sec the discussion in Sect. 3'3.2), a reference energy scale 
needs to be specified when quoting a value for a^. Because of the large amount of data 
collected at LEPl and SLC, it has become conventional to use = Mf , Mz being the 
Z-boson mass, provided that the measured value of is run from the scale at which the 
measurement takes place to Mz (from now on, will be used as a short-hand notation 
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for as{Mzj). This is in turn a great advantage since it allows to compare independent 
determinations of in different reactions (e+e" annihilation, hadron-hadron collisions 
and dccp-inclastic Icpton-hadron scattering) at different energy scales, covering a range 
of from roughly a few to 10^ GeV^. In particular, "low-Q^" and "high-Q^" results 
can be compared, thus providing a non-trivial check of the running of Ug as predicted 
by QCD, and a world- average value of ag, as obtained by averaging all the available 
measurements once extrapolated to the Z-boson mass, can be derived. 



e T decays 
D S Bjorken SR 
D S GLS SR 
D S R2 (NMC) 
D S R2 (HERA) 
D S Fj (SLAC.BCDMS) 
pp direct 7 
Lattice QCD 
M', T decays 
DIS Fj, X Fj (CCFR) 
DIS jet rates 
DIS ev. shapes 
pp bb prod, 
pp VV + 1-;et 
e* e R below 

ev, shapes 22 GeV 

e" ev. shapes 35 GeV 
e* ev. shapes 44 GeV 
e* e" ev. shapes 58 GeV 
e* e" Z° ev. shapes 
e* fragment, fns. 
e* e" R at 

Z° line shape 

e" ev. shapes ' 33 GeV 
ev. shiupes ' 61 GeV 
e^ e" ev. shapes 1 72 GeV 



0.06 0.08 0.1 0.12 0.14 

«s (Mz) 

Fig. 38. - The as{Mz) values and total errors from measurements at e^e~ , lepton-hadron and 
hadron-hadron machines, and calculations from lattice QCD. The results are ordered vertically 
in ^/Q^. 

The data for the electroweak parameters can be used to measure ag according to 
two different procedures, i.e. 

• performing a fit, with as as a free parameter, to the single observable given by the 
ratio of hadronic to leptonic width Rf, 

• performing a global fit to the whole set of the electroweak data (usually including 
the direct measurement of nit), by allowing the three parameters mt, mn and ag 

or the two parameters mu and a,, to vary. 

For the inclusive ratio Ri, the QCD corrections are known up to 0{a^g) for massless 




80 



G. MONTAGNA, O. NICROSINI and F. PICCININI 



Table VIII. - World averages of as{Mz) 



Value Error Reference 

0.118 0.003 ref. 

0.118 0.005 ref. 

0.118 0.004 ref. 

0.119 0.006 ref. 

0.119 0.004 ref. 

0.119 0.005 ref. 

0.119 0.004 ref. 



quarks, as discussed in Sect. 3'3.2, and the electrowcak contributions are well under 
control since the main part of the weak radiative corrections, being common to the 
leptonic and hadronic decay channels, cancel in the ratio of partial widths. Prom the 
experimental point of view, i?/, being a ratio, is known with a precision higher than that 
of single partial widths and is slightly affected by the luminosity uncertainty. Using the 
combined LEP result for Ri the fitted value of ag is [7] 

(136) as = 0.124 ±0.004(exp.)± 0.002(mH), 

where the experimental error is essentially due to the limited data sample and the the- 
oretical error (mn) derives from the residual dependence on tuh- It is worth notic- 
ing that the extraction of q,, from Ri doesn't require any specific knowledge about the 
hadronization mechanism. This as determination at the Z'^ peak is exactly analogous 
to the procedure followed when considering the published measurements for the ratio 
R = (j(e+e~ hadrons)/(T(e"'"e~ — ^ fi^ at energies below the Z^ resonance. The 
value obtained in such a case when fitting as to R as measured in the energy range 
5 < \/s < 65 GeV is [188, 192] 

(137) Us = 0.1281^:^13 ±0.002(mH), 

where the second theoretical error is due to the variation of mn in the range 60-1000 GeV. 

An example of the second kind of determination of as from the Z^ line-shape data is 
given by the three-parameter fit discussed at the end of Sect. 4'2.1, yielding the result [7,8] 

(138) as = 0.120 ±0.003(exp.). 

Detailed analyses of the theoretical uncertainties underlying the above determinations 

suggest that they can be estimated to contribute about ±0.002 [193]. Actually, for all the 
input observables the non-perturbative effects are expected to be of the order of 1/Mz 
and are hence usually neglected. 

It should be kept in mind that the results for as from Z^ data are obtained under 
the assumption that the hadronic width is given by the SM. Therefore, a possible 
anomaly in Ff, is a potential source of bias in the ag derivation. 

Besides the above quoted determinations of a^, other measurements performed at 
e^e~ , lepton-hadron and hadron-hadron machines need to be mentioned for their rele- 
vance, together with the calculations from lattice QCD. The interested reader is referred 
to the detailed compilations that can be found in refs. [188], [189], [190] and [8], together 
with references to the original literature. 
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A compilation of as{Mz) determinations is given in Fig. 38, together with the cor- 
responding total error. It is worth noticing that the most precise determination of is 
obtained from a global fit to the whole set of precision data (Z" line-shape determina- 
tion). By inspection of Fig. 38 it can be seen that there is agreement, within the errors, 
between "low-Q^" and "high-Q^" measurements. Contrary to a few years ago, where a 
discrepancy between low- and high-energy determinations was observed, the present sta- 
tus of as measurements is quite satisfactory, basically as a consequence of the fact that, 
while the ag determination from deep-inelastic scattering and lattice QCD increased, the 
precise measurement from fitting the clcctroweak data decreased. The problem of de- 
termining a world average of ag has been considered by several authors. In Tab. VIII the 
most recent central values and corresponding errors are quoted, together with reference 
to the original literature. 

4'3. Model Independent Approaches and Physics Beyond the Standard Model. — 
The possibility of extracting predictions for the mass of the top quark and of the Higgs 
boson, as well as for as{Mz), via the analysis of the electroweak precision data, as 
described in the previous sections, relies upon the assumption of the validity of the SM. 
In spite of the great success of such a strategy, it is of course important to exploit the 
high precision of LEPl/SLC measurements in order to test the validity of the SM in 
a way as independent as possible of the details of the underlying theory and, possibly, 
to derive also constraints on models of new physics beyond the SM. In order to achieve 
these goals, a model independent strategy for the analysis of precision data has to be 
organized. In the following, a model independent way of extracting parameters (the 
S-matrix approach) is shortly described. Next, the s parameterization as a tool for a 
model independent analysis of the Z° parameters is considered. At last, the implications 
of precision data for physics beyond the SM are briefly discussed. 

4'3.1. ^-matrix approach. The original idea behind the parameterization of the 
line shape introduced in ref. [113,148,149] was to allow a model independent extraction of 
the fundamental parameters of the resonance, such as its mass, width and decay rates. 
In the following years this procedure has been posed on a firmer ground by appealing to 
the S'-matrix theory [194-197]. Within this framework, the scattering amplitude can be 
cast into the form of a Laurent expansion around the resonance as follows [198-200]: 

(139) A{s) = + J2 Bn{s - SpT = +Bo + 0{V%/Ml), 

S Sp S Sp 

where the quantities R, Sp and Bq are complex numbers. Taking into account that an 
overall phase in the amplitude is unobservable, the cross section depends on five real 
parameters, which are separately gauge invariant, because the amplitude is an analytical 
function of s. The position of the complex pole defines the mass M and the width F in 
a gauge invariant and process independent way [201-203] through the relation 

(140) Sp = - iMT. 

In the expansion given by eq. (139), the photonic contribution to the amplitude is part 
of the background denoted by Bq. Another approach followed in the literature [204] is to 
introduce explicitly the photon exchange contribution according to the parameterization 



(141) 



A{s) = + + B{s), 
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which has been implemented in the computer code SMATASY. In principle the parameter- 
ization of eq. (141) may lead to some difficulties related to the fact that the coefficients 
R~f and Bi arc not independent quantities [200,204], but they disappear if a truncated 
version is employed. Neglecting terms of the order of F^/M^, the cross section for a given 
final state can be written in the following form [205] : 



(142) 



a{s) = -ira^ 



r1 srT + {s-M^)jT 

S (S-M2)2 + M2r2 



where is the photon exchange term and the parameters and jr are related to 

the Z-boson exchange residuum and to the j-Z interference, respectively. On the same 
grounds as for the total cross sections, the asymmetries can be characterized around the 
resonance by two parameters [206] : 



(143) 

with 
(144) 

(145) 



Ao = — + o{ryM^), 



Ax 



IA_JT_ 



where and ja are the analogues of tt and jt for the numerators of the asymmetries. 
In order to allow for a realistic data analysis, the universal efii'ects of ISR need also to be 
taken into account, following the procedures described in Sect. 3'4 and in Appendix A. 

The Z-boson mass and width parameters extracted with the S'-matrix approach 
described above difi'er from the commonly defined on-shell Z-boson mass and running 
width by two-loop and higher order corrections, resulting in a shift of 34 MeV and 1 Mev 
respectively, according to the following expressions [198,207-209]: 



(146) 
(147) 



p2 

M = Mz 77- ^Mz- 34MeV, 

2Mz 



T = Tz 



^ Z 

2M| 



~ Tz - IMeV. 



Taking into account of the proper definitions, the .Z-boson mass and width obtained 
by fitting the data with the S'-matrix formalism are in agreement with those obtained 
within the SM, showing its internal consistency, even if the error on the Z-boson mass is 
larger in the S'-matrix approach, because the 7-Z interference is a free parameter. The 
fact that the 7-Z interference is sizeable off peak can be used to reduce the error on the 
parameter 7^"''. Actually, by adding to the peak results also data above and below the 
resonance, the values obtained from such an analysis are Mz — 91.1882 ± 0.0029 GeV 
and jl^"'^ = 0.14±0.12 [8]. These values are to be compared with Mz = 91.1867±0.0020 
from the standard analysis, and j^"'' = 0.22 as the SM expectation. 

4'3.2. Model-independent analysis of precision data. Since new physics effects can 
be more easily disentangled in the analysis of the Z° parameters if not obscured by our 
ignorance of SM parameters, it is convenient to introduce variables that are as free as 
possible from large mt effects and hence sensitive to new physics. This attitude is fol- 
lowed in the model independent approach proposed and further developed in ref. [210], 
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where the vahdity of the SM is not assumed ab initio and appropriate variables (the e 
parameters) are defined in order to provide an efficient parameterization of the most im- 
portant clcctrowcak data with respect to the sensitivity to new physics. More specifically, 
four independent quantities, indicated as ei,e2,£3 and £(,, are introduced in one to one 
correspondence with the observables Mw/Mz,ri,A''pg (assuming lepton universality) 
and Th, that arc chosen as primary defining measurements of the e parameters. From 
these input data, the strategy consists in isolating the key quantities Ap,Arw,^k,Sb, 
that, as discussed in Sect. 3'2, are the dominant effects in weak radiative corrections 
due to gauge bosons self-energies and vertex corrections to the Z ff coupling. Indeed, 
the four quantities Ap, Ar^y, Afc, £(,, whenever calculated in the SM, are dominated, for 
sufficiently large top-quark mass values, by quadratic terms in rrit of the order of Gf_i,m^. 
In particular, Ar^/ is connected to Ar defined by cq. (56) in Sect. 3'2 according to 
the definition (1 — Ar) = (1 — Aa)(l — Arw), in such a way that the running of a 
{a{Mz) = a/(l — Aa)) due to known physics is extracted from Ar. Furthermore, £b is 
introduced in order to take care of the important non-obliquc corrections to the Z ^ bb 
vertex discussed in Sect. 3 '2. The explicit relations between the primary defining observ- 
ables and the Ap, Arw, A.k,Sb factors can be found in the original literature [210,211]. 
Since the aim, as already stressed, is to provide a parameterization unaffected as much 
as possible by the relative ignorance of rrit, it is convenient to keep Ap {ei = Ap) and £(, 
and introduce, in place of Ar;^ and Afc, the two following linear combinations 

£2 = CqAp -h -2 -2 



Cq 5q '2iS^^.k 



(148) es = clAp+{cl-sl) Ak, 

where ,Sq is an effective weak mixing angle containing only photon vacuum polarization 
effects, given by 

2„2 _ MMz) „2_i 2 

(149) soco- ^^^^2' ^0-1 ^0- 

In this way, contributions of the order of G/imf do not enter £2 and £3, that contain 
only logarithmic terms in m*. Quadratic top-quark mass effects are confined in £1 and 
Sb- Furthermore, the leading logarithmic terms for large Higgs-boson mass arc contained 
in £1 and £3. For instance, the leading expressions for the e parameters within the SM 
can be written as [210,211] 

£1 = — ^— rfr „ JT tan vw m -n- + • • • 

G^^Mlr , ( mt 



^ G,M^ fmH\ G^M^ fmt\ 
' 12n^^/2 \MzJ 6n^^/2 \Mz J ^ ' 

Finally, the relations between the defining observables and the £j can be inverted in order 
to get the formulae for the Si in terms of the data. The explicit formulae read as follows 



£1 = -0.9882 -h 0.011963 Fj/MeV- 0.1511a;, 
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Fig. 39. - Electroweak precision data versus theory in the ei-es plane (from ref. [154]; see also 
refs. [210,211]). 

£3 = -0.7146 + 0.009181 r^/MeV- 0.69735 a;, 
£6 = -0.62 £i + 0.2463 + 0.436 {Thl^m - 1) , 
(151) £2 = l-43ei- 0.8653 + 0.43 ArvF, 

where X = gv/SA = 1 — 4(1 + Afc)so and r^o is the value of in the Umit when all the 
e's are neglected. 

The parameters si, S2 and £3 can be directly related to the variables S, T and U, that 
arc other quantities proposed in the litcratiue in order to parameterize the oblique correc- 
tions [212] and widely employed for model-independent studies of electroweak data [213]. 
However, while S, T and U are defined as deviations with respect to the SM predic- 
tions for specified values of m.t and mn, the s parameters are defined with respect to a 
reference approximation that is independent of mj, in such a way that they are exactly 
zero in the SM in the limit of neglecting all pure weak loops [210]. Moreover, it is worth 
noticing that the definitions for the e parameters are quite general because they do not 
refer to any particular model. Therefore, they are useful since they allow to perform a 
model-independent analysis of the electroweak precision data. 

As a first step of such an analysis, it is possible to derive Icr contours for the e 
parameters by using the minimal set of data given by Mw/Mz, F;, A'-p^ and F;,. To 
include additional observables in the e's analysis, further assumptions are required in 
order to maintain a consistent definition of the parameters. However, under appropriate 
hypotheses, the analysis can be generalized to include all the observables measured on 
the peak, supplemented with low-energy electroweak measurements. The results 
obtained from the data can be then compared with the SM predictions, as a function of 
rrit and m^r, obtained by using computational tools that contain the state-of-the-art of 
radiative corrections. Examples of such a procedure are given by Figs. 39-40. Prom the 
comparison between theory and experiment in a given Si — Ej plane, it is possible with 
the present data to draw a number of significant conclusions. First, a preferred range for 
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Fig. 40. - The same as in Fig. 39 in the et-ss plane (from ref. [154]; see also refs. [210,211]). 
The effects of the errors on a{Mz) and as{Mz) are also shown. 



the fop-quark and the Higgs-boson masses can be inferred, according to a strategy that is 
different from what has been described in previous sections. Referring to the most recent 
data (with the TEVATRON result for nit excluded from the analysis), one can see a 
preference of the precision data for a "light" Higgs boson and a rrit value somewhat lower 
than the CDF/DO determination. These conclusions on mt and tuh corroborate the 
results obtained from SM fits discussed in Sects. 4'2.1 and 4'2.2. Second, it is possible to 
observe a good agreement between the full SM predictions and the data. Noticeably, one 
can point out a strong evidence for pure weak radiative corrections, since the beautiful 
agreement between data and theory can not be simply explained in terms of an improved 
Born approximation based on tree-level SM plus QED and QCD corrections. This very 
important conclusion demonstrates by itself the constraining power of precision physics 
at LEP. The latter has been also emphasized by some authors that have studied the 
sensitivity of the data to a very peculiar subset of weak radiative corrections, namely the 
pure bosonic loops involving trilinear gauge-bosons and Higgs couplings [172]. As shown 
by these investigations, the LEPl/SLC data are so precise that they require the inclusion 
of such corrections in the SM calculation in order to precisely fit the data. Therefore, 
the present LEPl/SLC precision data feel the electroweak non-abelian couplings that are 
presently measured at LEP2 [184]. 

4'3.3. Physics beyond the Standard Model. As already stressed, the main virtue 
of the £ parameterization is that they can be rather easily calculated in extensions of 
the SM, and therefore can be used to analyze models of new physics in the light of 
precision data. Typical examples studied in the literature are technicolour models [214], 
models with an extended gauge group [215,216] and the Minimal Supersymmetric Standard 
Model [217]. Also specific analyses within specific models of new physics beyond the SM 
have been considered in the literature, for instance concerning models with an extended 
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Higgs sector [218-221] or models with extra fermionic generations [222,223]. The main 
conclusions of such analyses are summarized in the following. 

Concerning technicolour models, it can be said that they tend to produce, in their 
typical realizations, corrections to the s parameters (large and positive to £3, large and 
negative to £(,) that are disfavored by the data [224]. In short, simple technicolour 
models are basically disfavored by LEP experiments [225], even if they can not be ruled 
out completely since particularly sophisticated versions with non trivial behaviour under 
the electroweak group could avoid their typical bad consequences [226] . 

As far as models with an extended gauge group are concerned, the case of the simplest 
models with an extra U{1) (and with an associated new neutral vector boson) has been 
addressed in the literature [227]. For such a situation, two new parameters have to be 
introduced, the former being a mixing angle ^ in order to define the mixture of the 
standard and the new neutral vector bosons, the latter being the shift 6pM induced at 
the tree-level in the p parameter by the above mixing. Considering the class of models 
based on E6 and for the left-right symmetric model, the main implication of precision 
data is that very strong constraints on ^ and 6pM can be deduced. In particular, the 
amount of mixing allowed is very small, less than, say, 1%. It is worth observing that, 
in recent years, models containing an extra Z' boson with enhanced hadronic or almost 
vanishing leptonic couplings (hadrophilic or leptophobic Z') received particular attention 
as candidate models able to explain the simultaneous anomaly of Ri, and Rc data (see 
for instance rcf. [155]). At present, the Rt-Rc anomaly has substantially disappeared, 
but anyway these model have not been completely ruled out by the new data, and could 
be tested at LEP2 [228]. 

The Minimal Supersymmetric Standard Model (MSSM) deserves a special discussion 
since it is the most predictive framework beyond the SM. However, because of the very 
high number of free parameters, a convenient strategy when confronting theory and 
experiment is restricting to consider two limiting cases (the so-calk^d "heavy" and "light" 
MSSM) rather than attempting a direct fit to the data. The "heavy" MSSM implies all 
supersymmetric particles to be sufficiently massive. In this limit, the MSSM predictions 
for electroweak data essentially coincide with the results of the SM with a light Higgs 
boson, say, lighter than 100 GeV [229]. Therefore, this particular realization of the MSSM 
can nicely accomodate the LEP data essentially in the same way as the SM does. In the 
"light" MSSM the masses of some of the superparticles are nearby their experimental 
lower limit. For such realization, the pattern of radiative corrections differs from that 
of the SM, with peculiar effects in vacuum polarization diagrams and Z ^ bb vertex, 
that can be easily incorporated in the e's [230]. These peculiar effects are reviewed in 
some detail in [231]. In particular, it is possible to find mechanisms able to explain a 
possible departure (if real) of Rb from the SM [232]. In conclusion, the MSSM agrees 
well with precision electroweak data, since it doesn't give rise, both in its "heavy" and 
"light" realizations, to effects that are inconsistent with the LEP measurements. This 
is due to the particular nature of the MSSM as new physics model that doesn't alter 
significantly the structure of the SM. This positive outcome of precision physics at LEP 
pointing towards supersymmetric extensions of the SM is corroborated by independent 
investigations of the MSSM with respect to the electroweak data performed in [160]. 
In these papers, instead of analyzing the MSSM implications on the data in terms of e 
parameters, a more orthodox approach has been followed, namely a complete calculation 
of electroweak precision observables has been carried out, using one value for the Higgs- 
boson mass and tan/?, together with the set of SUSY soft breaking parameters fixing 
the chargino/neutralino and scalar fermion sectors. In this kind of analysis, an optimal 
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set of SUSY parameters can be extracted from a global fit to the data, showing how 
supersymmetry can consistently describe precision data without contradicting present 
limits from direct search on its mass spectrum (see Fig. 41). 
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Fig. 41. - Experimental data normalized to the best fit results in the SM and MSSM (from 
ref. [160]). 



Other specific models analyzed in the light of the LEP/SLC data concern the ex- 
tension of the Higgs-boson sector. In the minimal SM the tree- level masses of the vector 
bosons are linked by means of the weak mixing angle through the relation p = I, with 
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Since the Higgs-boson sector has not yet been tested with high accuracy, it is interesting 
to investigate a possible extension of the SM adopting an enlarged scalar sector. The 
simplest way to do this, by changing the trcc-lcvcl expression of p, is to introduce, in 
addition to the minimal doublet representation, a Higgs-boson triplet with a vacuum 
expectation value different from zero in the neutral sector. In this case the masses of 
the vector bosons arc no more related by the weak mixing angle, and the p parameter 
deviates from unity already at the tree level according to the following form: 



where Ii{l3i) is the weak isospin (third component) of the i— th Higgs-boson multiplet, 
and Vi the respective vacuum expectation value. In addition to the standard Higgs boson, 
other Higgs particles belong to the physical spectrum. The scenarios arising in non 
minimal standard models with Higgs-boson triplets has been studied in rcfs. [218 220]. 
The calculation of the radiative corrections to precision observables within these non 
minimal models requires not only the evaluation of the extra loop diagrams involving the 
non standard Higgs bosons, but also an extension of the rcnormalization procedure. Due 
to the fact that the bare W- and Z-boson masses are not linked through the weak mixing 
angle, a fourth input data is necessary in order to perform the relative subtraction. In 
other words the p parameter is not calculable as in the minimal standard model but 
need to be fixed by a rcnormalization condition. The procedure has been discussed in 
ref. [51] within the MS scheme using the physical W^-boson mass as the fourth data 
point, in addition to a, and Mz- Recently, the radiative corrections to the full set 
of LEPl electroweak observables have been calculated in a triplet model with a neutral 
Higgs boson and a pair of charged Higgs particles in addition to the standard Higgs 
boson [221]. The adopted rcnormalization framework is an extension of the on-shcU 
scheme, where the fourth input data has been chosen to be the effective leptonic mixing 
angle at the resonance. The interesting result of this study is that for a variation 
of the top-quark mass in the range = 175 ± 6 GeV, the predictions of the SM and 
of the triplet model are practically indistinguishable. As far as the Higgs-boson mass 
dependence is concerned, while the SM has a preference for a heavy Higgs boson from 
the observable Tz and for a light Higgs boson from the mixing angle measurement, in 
the triplet model a light Higgs particle is compatible with all precision observables. 

Models with extra fermion generations have been also investigated. As previously 
discussed, one of the achievements of the LEP experiments is the determination of the 
number of light neutrinos {N^, = 2.993 ± 0.011) and hence of the number of fermionic 
generations. However, the theory does not exclude the presence of new sequential quark- 
lepton generations with heavy neutral leptons, whose production would be kinematically 
forbidden at LEP. The experimental lower bounds on the masses of the hypothetical 
new fermions available in the literature (including LEP 1.5 data) are of about 60 GeV 
for the leptons and 100 GeV for the quarks [222]. In the literature, the indirect bounds 
obtained through the study of their effects on radiative corrections for LEPl electroweak 
observables [222,223] have been examined. A large mass splitting between the members 
of the doublets would spoil the agreement with the precision data, so the new heavy 
fermions should be almost degenerate in mass within a doublet. In order to reduce the 
number of parameters, a reasonable assumption is to work with quarks and leptons of 
equal masses. Within this scenario, the effect of one new generation can be compensated 
by increasing the fitted value of the iop-quark mass. An analysis based on older data 
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sets shows that only at most two new sequential generations are allowed [222] . However, 
as a general caveat, it should be noticed that new experimental data, both from LEP2 
and the TEVATRON, could modify this scenario. Considering the effects of elcctrowcak 
vacuum stability and the absence of a Landau pole in the Higgs-boson potential allows 
to put upper bounds on the masses of the new fermions, depending on the scale at 
which new physics will appear. In the case of absence of new physics up to the grand 
unification scale, the upper bounds on the masses of the fourth generation are of the 
order of 100 GeV, a value viable for LEP2 direct searches [222] . 

As a last comment, it is worth pointing out that the present precision elcctrowcak 
data coming from the LEP analysis can be used to derive constraints on the structure 
of a grand unified theory (GUT). In particular, within standard GUT's the evolution of 
the gauge couplings starting from the scale given by the Z-boson mass Mz does not lead 
to the unambiguous identification of a grand unification scale. Viceversa, such a grand 
unification scale can be identified within the framework of supersymmetric GUT's [233]. 
This is at present considered as an indication of new physics beyond the SM, and in 
particular of supersymmetry [154]. 

5. - Physics at LEP2 

As discussed in the Introduction, the ALEPH, DELPHI, L3 and OPAL experiments 
at LEP terminated in the fall 1995 their data taking at energies around the mass of the 
resonance. After a short run at intermediate energies (the so-called LEP1.5 phase), 
the collider enters a new phase of operation far away the Z° peak, namely in the range 
^/s ~ 161 — 192 GeV. The energy upgrade of the LEP machine offers the possibility of 
exploring a new energy regime for electron-positron collisions, beyond the threshold of 
H'^-pair production (see Fig. 42). The precision measurements of Z-boson properties at 
LEPl/SLC can be therefore completed by the experiments at LEP2 with the precision 
determination of the VF-boson properties. In addition to the precision physics of the W 
boson, experiments at LEP2 will continue to perform tests of the Standard Model (SM) in 
two-fcrmion production processes and in photonic final states, search for the Higgs boson, 
perform QCD studies as well as explore direct new physics signals [15]. Actually, as can be 
clearly seen from Fig. 42, there are at LEP2 many interesting processes with measurable 
cross sections, such as e+e^ ff, e+e^ 77, e+e" vv^, e+e" W~^W~ , ZZ, the 
latter yielding four-fermion (4/) final states. 

In the present section, the processes of main interest for precision studies of the 
elcctrowcak theory at LEP2 are examined. The two-fermion production processes and 
radiative events are discussed, pointing out new physics aspects with respect to LEPl. 
The physics of the W boson in 4/ production processes is then addressed, with particular 
emphasis on the status of theoretical predictions as well as first experimental results. The 
search for the SM Higgs boson in 4/ final states is finally discussed. 

5'1. Two-Fermion Processes. - Although LEP2 was mainly conceived to measure 
precisely the properties of the W boson through the production of W pairs in e+e~ 
collisions, the processes of quark and lepton pair production, still remaining among the 
most copious reactions at LEP2, is an important physics item. In this new energy 
regime, the two-fermion processes can be used to perform further tests of the SM, with 
enhanced sensitivity to the electroweak interference, as well as to obtain limits on possible 
extensions of the SM, such as the existence of new Z' bosons or effective 4/ contact 
interactions [235]. However, the de-convoluted cross sections for these processes are 
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Fig. 42. - Standard Model cross sections at LEP (from [234]). 



typically of the order of 10 pb in the LEP2 energy regime, i.e. three orders of magnitude 
smaller than the cross sections for the same processes at LEPl. Therefore, even when 
considering the highest final integrated luminosity expected at LEP2, i.e. 500 pb""'^ [236], 
it will be possible to measure the two-fermion observables with a statistical error of 
the order of 1%, to be compared with the corresponding 0.1% of LEPl. Nonetheless, 
in order to meaningfully compare theory and experiment, precise predictions arc still 
demanded and, therefore, radiative corrections need to be known with a good precision 
and incorporated in the theoretical calculations. To this aim, it is worth noticing that 
at LEP2 new aspects in the sector of radiative corrections to two-fermion processes do 
appear with respect to LEPl, noticeably [234] 

• weak boxes {i.e. diagrams containing WW and ZZ internal lines) are at the level of 
a few per-cent of lepton and hadron cross sections, and hence no longer negligible 
as at LEPl; 
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• initial-state radiation (ISR) of hard photons is not inhibited as around the Z° 
peak, where, as discussed in Sect. 3'4, the Z-boson width Tz acts a natural cut-off; 
consequently, hard-photon bremsstrahlung corrections become numerically relevant 
at LEP2. 

The phenomenon of the so-called Z'^ radiatkie relurn constitutes the main man- 
ifestation of the relevance of hard-photon radiation in two-fermion production above 
the peak. It corresponds to the emission of very energetic photons that reduce the 
two-fcrmion effective centre of mass (cm.) energy back to the Z-boson mass, and con- 
siderably enhance the cross section. Actually, one can observe an increase of about a 
factor of four in the total rate of two-fermion events at LEP2, between the tree-level and 
the QED corrected predictions. This effect directly arises from the dependence of the 
two-fermion cross section on the cm. energy, when going from LEPl to LEP2. Actually, 
the Breit-Wigner behaviour for ^/s ~ Mz, followed by a fast 1/s decrease after it, makes 
very likely the production of a // pair with invariant mass m^j clustered around Mz, 
accompanied by the emission of a real photon carrying away the residual energy. In fact, 
the photon energy spectrum shows a pronounced peak around the value (1 — M|/s) \/s/2, 
corresponding to those events (LEPl-like events) with an effective cm. energy after ISR 
of a/I ~ Mz {Z'^ radiative return) . 

The above comments are referred to s-channel annihilation processes. As far as 
large-angle Bhabha scattering is concerned, it is worth noticing that already some GeV 
off-resonance the cross sections is dominated by t-channel photon exchange. Hence, large- 
angle Bhabha scattering at LEP2 is much more similar to small-angle than to large-angle 
Bhabha scattering at LEPl (see ref. [18] for a detailed discussion). 

S'l.l. Weak corrections. As already discussed in Sect. 3'2, the one-loop pure weak 

corrections to e~^e~ jZ'^ f f arise from those diagrams that involve corrections to 
the vector boson propagators, from the set of vertex corrections (with the exclusion of 
the virtual photon contributions that are accounted among the QED effects) and from 
the box diagrams containing the exchange of two massive gauge bosons. In fact, the 
box diagrams where at least one boson is a photon are classified as belonging to the 
QED corrections simply because they are IR divergent and need to be combined with 
the photonic initial-final state interference in order to get a meaningful result. Therefore, 
the genuine weak boxes are those with two W- and Z-boson propagators. Differently 
from the propagators and vertex corrections, they are ultra-violet finite in any gauge 
but the unitary one, where they actually give rise to UV infinities. Moreover, weak 
boxes introduce, in addition to s, a dependence on the scattering angle, that makes 
the parameterization of weak corrections in terms of "running" (s-dependent) effective 
couplings, successfully employed at LEPl, no longer viable at LEP2. In fact, since the 
box diagrams are negligibly small around the Z° pole, any dependence on the scattering 
angle can be neglected in the calculation of pure weak effects around the peak. Off 
resonance, i.e. in the LEP2 energy regime, the weak boxes become relevant corrections, 
due to the natural appearance of WW and ZZ thresholds. Therefore, they are actually 
included as relevant short-distance effects in the theoretical predictions for two-fermion 
physics above the Z^ peak. 

In order to make some quantitative statements on the effects of weak boxes on two- 
fermion observables, it is necessary to remember that the contribution of weak boxes is 
not gauge invariant by itself. It follows that a sensible evaluation of these corrections, 
as well as a comparison of their size as obtained from different computational tools 
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based on different theoretical framework, can not be carried out unless a gauge-invariant 
procedure of subtraction of weak boxes (the so-called de-boxization) is introduced. A 
possible procedure has been proposed in rcf. [234]. Denoting the correction due to the 
WW box diagrams by Bww{0 computed in a general gauge, then B\yw{0 can 
be split as 

(154) BwwiO = Bwwil) + ie - 1) A(e). 

In ref. [234], in order to get an estimate of the effect of weak boxes at LEP2, it was 
agreed to subtract from the full one-loop theoretical prediction Bww{£, — !)• This 
contribution at LEPl is of the order of 10~^ and its evaluation far from the peak can 
give an idea of the numerical effects introduced by the weak boxes at LEP2. The effect 
of weak boxes, as defined above, has been evaluated for the most relevant two-fermion 
processes (e^e~ — > ^~ ,uu,dd,hh, hadrons) and for several cm. energies. The main 
conclusion of this analysis, based on the predictions of different tools such as TOPAZO, 
ZFITTER and WOH (see Sect. 3'5), is that the weak boxes introduce a significant relative 
correction of the order of a few per cent, especially around the WW threshold and at the 
highest LEP2 energies. In particular, for the effect is positive, of the order of 1-2% at 
i/i = 161 — 175 GeV and negative, of the order 1% at 205 GeV; for ahad the correction 
is positive at ^/s = 161 GoV, of the order of 1%, negative and around 3% and 4% at 
175 and 205 GeV, respectively. For the muon forward-backward asymmetry the effect of 
weak boxes is within 0.01. 

5T.2. QED corrections. Because of the overwhelming effect due to ISR in the LEP2 
energy regime discussed above, some contributions associated to the emission of hard 
photons that arc irrelevant at LEPl play a new role far above the peak. Actually, 
the soft-photon dominance in the ISR mechanism emphasized in Sect. 3'4 for the energy 
region around the resonance ceases to be valid at LEP2, and previously negligible 
hard-photon brcmsstrahlung effects, both at the leading and ncxt-to-lcading level, show 
up. This implies the need for an upgrade in the treatment of QED corrections with 
respect to LEPl, since the numerical impact of the new corrections is typically at the 
1% level and hence significant in view of the LEP2 experimental precision. 

The first example of such effects due to hard-photon bremsstrahlung is given by the 
0{a) NLO corrections arising from the emission of a hard acollinear photon by the ini- 
tial state, because hard collincar photons are already taken into account by the leading 
logarithmic (LL) prescriptions (see Appendix A). Since these corrections are strongly af- 
fected by the imposed experimental cuts and can be very difficultly cast in a fully analytic 
form, one can evaluate their size by performing an exact 0{a) perturbative calculation, 
integrating it numerically and depurating the final result of the universal LL as well as 
of the non-leading (process dependent) soft -|- virtual corrections [237]. This procedure 
allows to isolate and quantify in a clean way the effect of hard-photon contributions in 
0{a) ISR. The results of such an analysis can be found in ref. [237] for typical event 
selections (ES's) adopted by the LEP Collaborations, namely 

• only a cut on the invariant mass of the event after ISR; this matches the ES adopted 
by recent analyses of hadron and lepton-pair production by the LEP Collaborations 
at LEPl. 5 and LEP2 energies; 

• fermion acceptance (40° <"!?-< 140°) and acoUinearity (acoll=10°, 25°) cuts in 
association with an invariant mass one. 
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Fig. 43. The theoretical error of the "Z -peak recipe" for the fj, pair production cross section 
(relative deviation), as a function of the invariant mass cut Xcut defined by s' > XcutS, s' 
being the invariant mass of the final-state fermion pair after ISR. The results correspond to two 
acoUineaxity cuts, for the same angular acceptance cut. The arrows point at Xcut = M%/s. 



For the second ES, the results are shown in Figs. 43 and 44. They refer to four cm. 
energies: ^/s = Mz (as a test case), 140, 175 and 190 GeV. In all the cases, the final 
state considered is 

Whenever a s' cut alone is imposed in the data analysis (that in practice coincides 
with what is done by the LEP Collaborations above the peak), 0{a) non- leading hard- 
photon corrections are compatible with zero for hadronic and Icptonic cross sections, and 
negligible at the 10~^ level for the lepton asymmetry. This applies both to loose and 
tight s' cuts, i.e. including or excluding the events with radiative return. Therefore, 
for such typical ES's, those LEPl calculations accoimting for soft plus virtual NLO 
terms only (Z°-peak recipe) can be safely extrapolated at higher energies without loss 
of accuracy. For more complicated ES's (see Figs. 43 and 44), including an acceptance 
and acoUinearity cut, the theoretical error due to the neglect of the 0{a) non-leading 
hard-photon terms can grow up to about 2.5% for the cross section and to 1 x 10~^ 
for the forward-backward asymmetry. However, if the Z^ radiative return is cut away, 
the above theoretical error is limited to 0.3% relative deviation for the cross section and 
3 X 10~^ absolute deviation for the forward-backward asymmetry, that are acceptable in 
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Fig. 44. - The same as in Fig. 43 for the forward-backward asymmetry (absolute deviation) . 



the light of the LEP2 experimental precision [237] . 

The second example of hard-photon contributions damped at LEPl but no longer 

negligible at LEP2 comes from the IS second-order non-leading corrections and third- 
order LL contributions. The former arise from the configurations where a hard photon 
is radiated off in the direction of incoming electron or positron in association with a 
large-angle, acoUincar additional hard photon, and provide corrections of the order of 
a^L. The latter are due to the emission of three hard photons collinear to the colliding 
leptons, originating terms of the order of O{0^). The 0(a^i) non-soft effects are exactly 
known from explicit pcrturbativc calculation of the ISR spectrum [117]. They can be 
easily included in the radiator or flux function (see Appendix A) in such a way that 
the exact 0{a^) calculation for an inclusive cross section is reproduced. The LL O{0^) 
contributions, being universal photonic effects, can be kept under control employing one 
of the algorithms described in Appendix A. Referring in particular to the structure 
function (SF) method, these LL O{0^) corrections are known for the electron SF [119, 
238,239] as well as the radiator [118,120], both in factorizcd and additive form. It 
has to be noticed that the 0{a^) non- leading non-soft corrections are in principle of 
the same order of magnitude as the third-order LL 0(/3^) corrections. Therefore, a 
careful evaluation of the effects induced by higher-order hard-photon corrections in the 
LEP2 energy regime necessarily requires that both the contributions are simultaneously 
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Fig. 45. - The effect of 0(/3^) and 0{a^L) corrections on the e^e — > /x^/j total cross section 
as a function of the s' cut (see ref. [120] for more details). 



included in the theoretical predictions. The interplay between 0{a'^L) and 0(fi'^) effects 
for the QED corrected muon-pair cross section, as a function of the s' cut defined as 
s'/s > Xcut, is shown in Fig. 45. The relative deviations with respect to the cross section 
computed by means of a standard additive radiator with up to 0{a'^) LL hard-photon 
corrections (0-2) are shown for several cm. energies. As can be seen, both the 0{(3'^) 
(upper window) and NL £)(a^) (lower window) corrections amount to a contribution 
of several 0.1% when the Z'^ radiative return is included, but they tend to compensate 
one another, being of the same order of magnitude but of opposite sign. When the Z'^ 
radiative return is excluded, or close the Z^ resonance, the NL 0{a^) corrections are 
confined at the level of 0.01-0.02%, whereas the 0{l3'^) ones remain at the level of 0.05- 
0.1%. More precisely, the inclusion of the NL 0{a^) plus 0{P^) corrections above the Z° 
peak causes a reduction of the QED corrected muon-pair cross section of about —0.1% 
when the Z^ radiative return is excluded and an enhancement of about 0.25% when it is 
included. In the case of the hadronic cross section the effect of enhancement is around 
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Fig. 46. - Data-theory comparison for two-formion production cross sections at LEP. The 
experimental data can be found in [156] (on peak) and [240] (^i > 130 GeV). The theoretical 
prediction is performed by TOPAZO [111]. 



0.4%, i.e. comparable with the expected experimental precision for such an observable. 

The theoretical predictions for cross sections and forward-backward asymmetries of 
two-fermion processes, as obtained by means of the code TOPAZO including all the weak 
and QED corrections discussed above, are compared with experimental data and shown 
as functions of the cm. energy in Figs. 46 and 47, respectively. 

5' 1.3. Radiative events. As already discussed in Sect. 4, the study of single photon 

production in the process e+e~ i/Pj has been exploited at LEPl as an alternative 
method for the determination of the invisible Z-boson width. At LEP2, the radiative 
processes, both with one and more final-state visible photons, still remain an important 
physics item, but within a different framework [234,241]. For instance, they are used as 
QED or SM tests or, more generally, as tools to investigate physics beyond the SM. As a 
general feature, all these processes at LEP2 are affected by an error which is dominated 
by statistics and is of the order of some per cent. This in turn implies that a theoretical 
knowledge of the cross sections with a theoretical error of the order of one per cent is 
mandatory, thus requiring the inclusion of all the phenomenologically relevant radiative 
corrections. 

The process e+e~ wf, n > 2, is well suited for testing QED. Actually, it is 
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marginally affected by weak or strong corrections, whose amount is below 1%. Since it 
provides a particularly clean signature, it can also be used to probe models predicting 
the existence of excited leptons or contact interactions. At present, the data- theory 
comparison confirms pure QED [242]. 

The production of single radiative events in the SM is dominated by the process 
e+e~ — > vi^'y, with a sizeable background of radiative Bhabha events with the electrons 
lost. All these events manifest themselves as a single visible photon accompanied by 
missing transverse momentum {^±)- It is worth noticing that, at a difference from the 
LEPl case where the radiative photon is essentially soft, at LEP2, due to the phenomenon 
of the Z'^ radiative return discussed above, the observable photon is essentially a very 
energetic one, and thus easily taggable (see Fig. 48). Standard computational tools 
employed by the LEP Collaborations for the evaluation of the SM cross sections are 
described in refs. [132], [244] and [243]. These events can also be used to probe processes 
beyond the SM, such as production of neutralino-pairs or excited neutrinos. Moreover, 
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Fig. 48. - The energy spectrum of the observed photon at LEPl and LEP2 energies. Numerical 
results by NUNUGPV [243]. 



they are in principle sensitive to WWj anomalous couplings (see Sect. 5'3.2 for more 
details). Typical examples of new physics are the following: 

• in the framework of the Minimal Supersymmetric Standard Model (MSSM), with 
R-parity conservation and the neutralino Xi as the lightest supersymmetric particle 
(LSP), the process of interest is e+e~ — > XiXi7 [245]; 

• in the framework of the so-called "no scale supergravity model" [246] , the neutralino 

Xi should be visible through its radiative decay into a photon and a gravitino G, 

Xi 7G, so that the process of interest is e+e~ ^GG. 

Unfortunately, both the processes are hardly detectable because of their too low rate. 
Actually, at present the data analyses do not point out any anomaly with respect to the 
SM predictions [247]. 

A radiative process that received particular attention during the last months is 
e+e~ — > 77 + ^_L. Actually in the MSSM, such signatures could be produced by the 
process e~^e~ X2X2 ^ 77XiXij where X2 is the ncxt-to-lightest neutralino decaying 
radiatively into Xi7- In the framework of the "no scale supergravity model", instead, 
the process of interest is e+e~ — > XiXi ~* 77GG. At present there is no evidence for 
physics beyond the SM, and the data are used to put constraints on the masses of the 
supersymmetric particles [247]. 
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5'2. Four-Fermion Processes. - As discussed in Sect. 5'1, two-fermion production 
processes certainly are an important physics item at LEP2, but the true novelty is rep- 
resented by the processes involving the production of VF-boson pairs (and their back- 
grounds), that at LEP2 energies are kinematically accessible with cross sections com- 
parable to the ones typical of two-fermion production. The study of PF-boson pair 
production can add information concerning two important items of the SM, namely the 
determination of the M^-boson mass [248] and the structure of the triple gauge-boson 
couplings [249]. 

At LEPl, given the set of input parameters a, Gf^ and Mz (see Sects. 3 and 4 for 
more details) all the other observables can be derived within the SM. In particular, the 
W-boson mass Mw satisfies the relation (see eq. (61)) 



which is essentially the SM prediction for the muon decay. In cq. (155), Ar represents the 
radiative corrections to the tree-level matrix element, and hence depends, in particular, 
on the to;>-quark and Higgs-boson masses, m* and tuh- At LEP2, the H^-boson mass 
is measured and therefore eq. (155) becomes a constraint that, via Ar, correlates 
and niH to one another. Of course, the more precise is the knowledge of the VF-boson 
mass, the stronger is the resulting constraint, which can eventually shed light on the 
Higgs-boson mass. 

Concerning the triple gauge-boson couplings, within the SM they have a specific 
form determined by the underlying SU{2) (g) f7(l) symmetry. Any departure from such 
a specific form (anomalous couplings) is in general responsible of a bad high-energy 
behaviour of the cross sections, spoiling the renormalizability of the theory. As a matter 
of fact, at LEP2 energies an observable very sensitive to anomalous couplings is the 
angular distribution of the W bosons in VK-pair production. 

In order to exploit the experimental information as a stringent test of the SM, it 
is clear that one needs reliable theoretical calculations of the various W^-pair production 
observables. It is also clear that the best knowledge is mandatory, i.e. the observables 
should be computed taking also into account all the relevant radiative corrections. 



(155) 




air 1 



V2M^{1 - M^/M^) 1 - Ar 




W 



Fig. 49. - The tree-level Feynman diagrams for on-shell VK-pair production. 



A first step in this direction consists in calculating the cross section for T4^-pair 
production in the so called narrow- width (NW) approximation, i.e. assuming that the 
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M^-boson width is zero and hence the W bosons are stable particles (see Fig. 49). A huge 
amount of work along this direction has been performed in the past [250-252] , and today 
the process of on-shell PF-pair production is known together with its full 0{a) radiative 
corrections. 

On the other hand, the W bosons are, as well known, unstable particles, so that 
the effect of their finite width is in principle important, especially around the threshold 
region. Actually, whereas the cross section for the production of two real W bosons is 
exactly zero below ^/s = 2Mw, when the effect of the finite width is taken into account 
the cross section must become non-zero also below ^/s = 2Mw, due to the possibility 
of producing off-shell W^s which subsequently decay (see Fig. 50). Moreover, the W 
bosons, being unstable, can be only intermediate states of the reactions, the true final 
states being four-fermion (4/) states. This in turn implies that, for a given 4/ final 
state, besides the processes concerning the production and decay of a W-boson pair, one 
has to consider also those processes that lead to the same final state, but via different 
intermediate states (background processes). 




140 150 160 170 180 190 200 210 



Vs (GeV) 

Fig. 50. - The effect of the finite VF-boson width on the inclusive cross section. 

The above comments point out the need of going beyond the NW approximation, 
both in order to obtain reliable cross sections and to study all the details concerning the 
W-boson decay products. 
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5'2.1. Tree-level calculations. A first attempt in the direction of taking into account 
finite- width effects can be found in the approach described in ref. [253] . If one is inter- 
ested in inclusive VF-boson decays, then the total cross section can be approximatively 
factorized into two stages, namely off-shell VF-boson pair production times their decay, 
as follows: 

PS riV^-V^)^ 

(156) ^i^) = Bfih^fsh dSip{si) dS2p{S2)(To{s,Si,S2). 

Jo Jo 

In eq. (156) ao{s, si, S2) represents the total cross section for the production of two 
W bosons with squared masses si and S2, respectively (the Feynman diagrams of the 
production process at the tree level arc still given in Fig. 49). Its explicit expression can 
be found in ref. [253]. -B/^/^ is the branching ratio for the decay channel W falb- 
p{si) is the weight factor describing the inclusive decay of a W boson of squared mass 
Sj. It takes into accoimt that the W bosons are produced off-shell and hence, rigorously 
speaking, they must be described by a propagator. The explicit expression of p{s) is the 
following 

(157) p[s)='-, v^^^^) 



7r(s-M2,)2 + M2,r(s)2' 

where r(s) is the total T4^-boson decay width. Equation (156) is a very simple and 
effective approximate description of the physics of M^-pair production. Also, the main 
factorizable QED initial-state corrections to it have been considered in the literature and 
can be foimd in rcfs. [254 256]. On the other hand, its main limit is that it can provide 
at most double-differential spectra, namely the combined invariant mass distribution of 
the W bosons. If one is interested in the full description of the final-state products, a 
complete 4/ calculation is mandatory. 

A natural extension of eq. (156) to the description of a full 4/ final state can be 
obtained by considering the diagrams of Fig. 51, where the W bosons are the interme- 
diate states of the reaction and hence the full dynamical information on the final-state 
fermions is contained. The diagrams considered constitute the so called CC03 class, 
and describe the full 4/ process in the so called double-resonant approximation. The 
complete calculation of the matrix elements corresponding to them, together with the 
appropriate four-body phase space, allows one to go beyond the "off-shell production 
times decay" approximation of eq. (156), and obtain the fully differential description of 
the process. 

On the other hand, as already observed in Sect. 5'2, a given 4/ final state can 
be obtained also via other sub-processes. For instance, if one considers a semi-leptonic 
channel, i.e. a channel in which one ly-boson decays into hadrons and the other one 
into leptons, and restricts himself to those channels that do not have electrons in the 
final state, then the same final state as in Fig. 51 can also be generated by the diagrams 
of Fig. 52. These diagrams (background diagrams), together with the ones of Fig. 51, 
constitute the so called CCll class, and are characterized by having a single W boson 
in the intermediate state, so that they are single-resonant. Of course, the full matrix 
element for such a final state is the sum of all the diagrams, and the squared matrix 
element includes also the interferences between double- and single-resonant diagrams. 

In the presence of electrons in the final state, the situation is more involved. Actually, 
the full matrix element receives contributions also from additional i-channel diagrams, 
so that the total number of diagrams becomes 20 (the so called CC20 class). 
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Fig. 51. - The tree-level Feynman diagrams for the CC03 class. 



A full classification of the various Feynman diagrams contributing to a give 4/ final 
state goes beyond the aims of the present paper. The interested reader is referred to 
refs. [234,257,258], and references therein, for a more detailed account. Here it is worth 
noticing that for particular final states, for instance fully hadronic final states, also neutral 
current (NC) and QCD backgrounds appear and become relevant. 

The calculation of the scattering amplitudes for 2 — > 4 processes is, already at the 
tree level, considerably more involved than the corresponding calculation for 2 — s- 2 
processes, typical of LEPl/SLC physics. There are two main reasons for this, namely 
the fact that a single Feynman amplitude for 2 ^ 4 is algebraically more involved and 
the fact that, typically, for a given final state there are much more Feynman diagrams 
contributing. The calculational techniques adopted in the literature can be classified as 
follows: 

• helicity-amplitude techniques: in this approach, the scattering amplitude for a given 

process, and for a given helicity pattern of the initial- and final-state fermions, 
is computed analytically as a complex number by exploiting the formal prop- 
erties of the spin projection operators; the squared modulus of the amplitude 
is then computed numerically; the approach, in all its actual implementations 
(see refs. [259-262] and references therein), is particularly powerful for massless 
fermions, albeit also mass effects can be taken into account; 

• automatic calculations: these approaches adopt both standard techniques for the 

evaluation of the squared matrix element [263] and the helicity amplitude formal- 
ism for the evaluation of the scattering amplitude [264], properly interfaced with 
software packages that render the calculation of cross sections almost automatic; 

• numerical evaluation of the generating functional for the connected Green's func- 
tions: it is a new method, presented in ref. [265], in which the scattering amplitude 
is computed numerically by means of an iterative algorithm starting from the effec- 
tive action of the theory and with no use of Feynman diagrams; it becomes strongly 
competitive with respect to standard techniques as the number of final-state par- 
ticles becomes larger and larger. 

5'2.2. Gauge invariance. Being in the framework of a gauge theory, as is the case 
of the SM, means that the calculations of physical observables must be gauge invariant, 





Fig. 52. - The additional eight troc-lcvcl Fcynman diagrams for the CCll class. Notice that in 
the semi-leptonic case (e"'"e~ udfi^v^, for instance) the additional diagrams are seven; in the 
fuUy-leptonic case {e^e~ — > t'^VtH^Vp,, for instance) the additional diagrams are six. 



or at least the gauge violations have to be confined well below the required theoretical 
accuracy. In the case of 4/ final states there are two sources of problems connected with 
gauge invariance in theoretical calculations (for a review see for instance ref. [266]). 

The first one originates when the matrix element of the process is calculated consid- 
ering only a subset of the Feynman diagrams contributing to the scattering amplitude. 
For example, in the case of W^-boson production, the natural extension from the on-shell 
approximation to the more realistic off-shell production and subsequent decay would be 
the calculation of the three Feynman diagrams of the on-shell case with the two final- 
state fermionic currents attached to the virtual PF-boson lines, i.e. the double-resonant 
diagrams of Fig. 51. In ref. [250] it has been shown that from the evaluation of the double- 
resonant diagrams in an axial gauge, gauge dependent terms arise, which exhibit a single 
pole structure in either of the two invariant masses and are cancelled when also the single- 
resonant diagrams (see Fig. 52) contributing to the process are evaluated. The single- 
and the non-resonant contributions are generally suppressed by a factor Tw/Mw — 2.5% 
with respect to the double-resonant ones, but they need to be taken into account in view 
of the required theoretical accuracy. 

Even when the calculation of all the tree-level diagrams contributing to a given 4/ 
final state is performed, the problem of gauge invariance is not yet solved because of the 
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presence of singularities in the phase space due to the massive gauge boson propagators, 
which display poles at = M^, being the invariant mass of the decay products of the 
unstable bosons. These singularities arc cured by the introduction of the finite widths of 
the gauge bosons, which shift the poles away from the real axis. However, in field theory 
the widths arise from the imaginary parts of higher-order diagrams describing the gauge 
boson sclf-cncrgics, rcsummcd to all orders. So the trcc-lcvcl amplitude is supplemented 
with only a subset of higher-order contributions and this can destroy the gauge invariance 
of the calculation. 

A pragmatic approach, suitable for the construction of a lowest-order event gen- 
erator, is the so called "fixed- width scheme", where the propagators — M^) are 
systematically replaced with — -t-iFM) also for space- like momenta. This pro- 
cedure has no physical motivation, because the propagator for space-like momenta does 
not develop any imaginary part. An improvement on this is the "running- width scheme" , 
where the width is a function of p^, equal to zero for space-like momenta, in agreement 
with the calculation of the imaginary part of the gauge boson sclf-cncrgy. 

Since the resonant diagrams are not gauge invariant by themselves, both of the above 
mentioned schemes violate gauge invariance, and their reliability needs to be checked by 
a truly gauge invariant scheme, even if it cannot be uniquely defined. In the literature 
several gauge restoration schemes have been studied. The simplest one is the "fudge- 
factor scheme" [267, 268] , which amounts to calculate the matrix element without any 
width (so that the calculation is manifestly gauge invariant) and introduce an overall 
factor of the form {p^-M^)/{p'^-M'^+iTM) for every singularity, in order to transform it 
in a resonance. The problem connected with this scheme is that when the double-resonant 
diagrams arc not dominant, as is the case of charged current processes at energies below 
and at the H^-pair production threshold, it can lead to large deviations [269]. 

Another possibility is the "pole scheme" [198,270-272], where the complete ampli- 
tude is decomposed according to the pole structure in gauge invariant subsets of double-, 
single- and non-pole terms. Introducing the widths in the pole factors does respect gauge 
invariance. At present there is some debate about the correct way of implementing this 
scheme and about its validity in the vicinity of thresholds. Since the pole scheme is a 
gauge invariant decomposition of the amplitude according to its degree of resonance, it 
can be used as a starting point for the evaluation of higher-order corrections, as will be 
discussed in Sect. 5'2.3. 

The most theoretically appealing solution to the gauge invariance problem is the 
"fermion loop scheme" [269,273-275], where the minimal set of one-loop Feynman dia- 
grams necessary for compensating the gauge violation caused by the self-energy graphs is 
included in the calculation. Since the lowest order decay widths of the gauge bosons are 
given by the imaginary parts of the fermion loops in the one-loop self-energies, accord- 
ing to this scheme the imaginary parts of all the other possible one-particle-irreducible 
fermionic one-loop corrections must be included. For the process e+e^ 4/, after the 
resummation of the vector boson self-energies, the only left out contributions are given 
by the imaginary parts of the fermionic corrections to the triple gauge-boson vertex, 
which have been calculated in ref. [274] in the limit of massless fermions. The calcula- 
tion has been extended in ref. [275] to take into account the complete fermionic one-loop 
corrections, including real and imaginary parts, and all contributions of the massive top 
quark. 

Due to the complexity of the fermion loop scheme, the CPU time needed for the 

calculation of the cross section increases considerably with respect to the one performed 
within the fixed or running width scheme, based on a naive treatment of the bosonic 
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widths. For this reason, from a phenomenological point of view, it is important to quan- 
tify the amount of gauge violating effects in 4/ processes of interest at LEP2. As a 
result of the detailed investigations pursued in refs. [274,275], the violations connected 
with the gauge group SU (2) are not relevant in the LEP2 energy range. They become 
important for energies reached at the next generation of linear colliders. The case of 
U{V)e.m. gauge invariancc violation is different, since the effects arc enhanced by a factor 
of 0{s/ml) in processes with an electron or positron in the final state, almost collinear 
to the incoming particle [268,269,273,274]. This happens, for instance, for the CC20 
process e'^e~ e~i>f.ud, which is particularly important for studying triple gauge-boson 
couplings. In ref. [274] it has been shown that the running width scheme for this pro- 
cess gives a totally wrong result when the minimum scattering angle of the electron is 
set to zero, unless it is improved with the fcrmionic corrections to the trilinear gauge 
boson coupling. On the other hand, the fixed width scheme, which does not violate the 
electromagnetic current conservation, gives reliable results. It is worth noticing that the 
results for the total cross section, obtained by means of different methods respecting 
current conservation, agree to 0{T^/M^). In the limit 0, where represents 

the space-like momentum transfer of the electron through the photon, the fermion loop 
correction to the standard Yang-Mills vertex can be written in the following factorized 
form [274]: 



and p± are the squared momenta flowing in the W^='=-boson propagators. This approxi- 
mate factorized prescription can be easily implemented in a Monte Carlo event generator. 
However, the validity of the above formula is limited to the LEP2 energy range. 

As a last comment, it is worth noticing that the problems connected with J7(l)e.m. 
gauge invariancc in the CC20 processes are rendered much more mild when a cut on the 
scattering angle of the outgoing electron is imposed. For a realistic cut of 10° the results 
obtained with the fixed width scheme are compatible with the ones obtained with the 
fermion loop corrections [274]. 

5'2.3. Radiative corrections. As already discussed in Sect. 5'2, the main physical 
motivations of the LEP2 experiments are a precise determination of the VF-boson mass 
and the study of the triple gauge-boson couplings. Since in the LEP2 energy range the 
integrated cross section is not very sensitive to the presence of anomalies in the Yang-Mill 
sector, the angular distributions are the best suited observables for the determination of 
the anomalous couplings. For this reasons the radiative corrections to both total and 
differential cross sections should be under control. Moreover, as will be illustrated in 
Sect. 5'3.1, one of the methods used for the VF-boson mass measurement (the direct 
reconstruction method) relies upon the knowledge of the energy loss due to ISR defined 
as 



Cfi 




where 




0, q^ < 0, 



(158) 
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which requires the best possible control of photonic radiative corrections. The other 
method (the threshold method) requires a good knowledge of the total cross section 
in the threshold region, i.e. a total cross section including all the relevant radiative 
corrections. However, it should also be kept in mind that, given the limited statistics of 
the LEP2 experiments with respect to LEPl, a theoretical relative accuracy of the order 
of 1% at the VF-pair production threshold and of 0.5% at higher energy will be enough. 

Before discussing various aspects of radiative corrections to 4/ processes, it is worth 
recalling the role played by the PF-boson mass at LEP2 with respect to the situation of 
LEPl. As already pointed out, in this last case it is an observable calculated by means of 
the input parameters a, and Mz, and, with the inclusion of the radiative corrections 
to the n decay, Mw shows a parametric dependence on mt, tuh and Ug, as results from 
eq. (61). In the case of LEP2 it is more natural to consider Mw as a free parameter to be 
fitted by the experiments. As a consequence, eq. (61) can be solved with respect to rritop 
as a parametric function of and Mh, to be compared with the direct measurement 
from the TEVATRON [8,276], allowing for a stringent determination of the Higgs-boson 
mass. 

Turning now the attention to LEP2 processes, the ideal situation would be the knowl- 
edge of the complete set of 0{a) radiative corrections to processes with four fcrmions 
in the final state, but at present the calculation, due to its complexity, is not available. 
In fact, such a task would require the calculation of 0(10^ 10^) one-loop diagrams. It 
is important to notice that, due to the presence of charged vector bosons in the Born 
amplitude, a separation between virtual photonic and pure weak corrections respecting 
the SU{2) gauge invariance is not possible. Anyway a class of dominant corrections im- 
portant for LEP2 physics can be singled out, and the effects of non-leading terms can be 
estimated relying upon the knowledge of the full one-loop electroweak corrections for the 
reaction e+e~ — » W~^W~ , where the charged gauge bosons are considered in the on-shell 
limit [251,252], including the bremsstrahlung process e+e~ W^W^'f [277]. Since 
the exact 0{a) calculation, already for on-shell W^s, leads to complicated and lengthy 
expressions in terms of twelve s- and t-dependent effective couplings, the possibility of 
constructing an Improved Born Approximation in the LEP2 energy range has been stud- 
ied [278], with the aim of providing simple and transparent formulae to be eventually 
implemented in Monte Carlo's for off-shell W^-pair production. 

Considering the realistic case of 4/ production, in the following a brief account 
is given on the state-of-the-art about the three main classes of radiative corrections of 
interest for LEP2 physics: weak, QCD and QED corrections. 

As already pointed out, at LEP2 the W^-boson mass Mw should be a fitting param- 
eter. By simply adding it to the three data points used as input parameters at LEPl (see 
Sect. 3'2.2) leads to a four data points scheme, based on a, G^, Mz and Mw, that is 
overcomplete, so that one of the typical LEPl input parameters has to be replaced with 
Mw- Considering eq. (61) with Ar = 0, at least two tree-level schemes are possible: 

= aV2G^M^, 

and 
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Fig. 53. - The contribution of various radiative corrections. The left window shows the 

cross section for the channel e^e^ fi'Uf^ud at the tree level (dashed line) and includ- 
ing ISR (dash-dotted line), the coulombic correction (dotted line) and naive QCD correc- 
tions (continuous line). The right window shows the relative effect of each correction. 
Aa/a is defined as (aisR — cr Born) /cr Born for the dashed line (the relative effect of ISR), 
[ffiSR+couL.—oisR) /oisR for the dotted line (the relative effect of the Coulomb correction) and 
{oiSR+couL.+NQCD -~ (Ti s R+cou L.) / (Ti s R+cou L. for the dash-dottcd line (the relative effect of 
naive QCD corrections). Numerical results produced by WWGENPV [279]. 




However, for practical calculations, the leading universal weak effects can be absorbed 
in the tree-level calculations by defining the input parameters in terms of and a{s). 
This choice, applied to the second scheme, is the one adopted in [258] . 

Going beyond this phenomenological attitude, it is necessary to talk about a com- 
plete one-loop calculation. By adopting the fermion loop scheme, all one-particle irre- 
ducible fermionic one-loop corrections are included, and a theoretically appealing strategy 
in this direction would be to calculate the remaining bosonic corrections, at least within a 
suitable approximation. The pole scheme [198,270,271] offers such a framework, provided 
that the cm. energy is sufficiently far away from the threshold region (a few H^-boson 
widths above). Actually, as stated in Sect. 5"2.2, this scheme is gauge invariant, contains 
the corrections for the production and decay of on-shell bosons as building blocks, and 
the number of Feynman diagrams to be evaluated is considerably reduced with respect to 
the complete set, because the potentially relevant corrections in the LEP2 energy region 
are given by the ones affecting the double-resonant diagrams. However, a simpler ap- 
proximate procedure to take into account the leading bosonic corrections in the on-shell 
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W's limit has been recently suggested [280], which amounts to perform the substitution 



with Aj/^j*^ = 11.1-10 ^. This replacement is equivalent to use the SU (2) gauge coupling 
g{M^) at the high energy scale of M\v, defined by the theoretical value of the radia- 
tively corrected leptonic W^-boson width, ^^(M^) — iSnT]^ /Mw-C) The essence of the 
approach is to use Tf^ instead of as input parameter, and it has been checked to 
work well in the case of on-shell VF-pair production, providing results with an accuracy 
better than 1% with respect to the full one-loop results. Moreover, in rcf. [280] it has 
been shown that strong cancellations occur between fermionic and bosonic corrections, 
at least at the level of on-shell VF-pair production. Since the same behaviour can be 
reasonably expected also in the case of off-shell VF-pair production [275], these results 
confirm that the fermionic loop contributions to e+e~ — > 4/ need to be completed with 
some approximation of the bosonic corrections. 

As far as the ly-boson mass determination by the threshold method is concerned, it is 
worth mentioning a specific bosonic correction arising from the exchange of a relatively 
light Higgs boson between the slowly moving W bosons. Given the presently allowed 
range of variation for the Higgs-boson mass, it induces a theoretical uncertainty on the 
total cross section near the threshold of VF-pair production. In the limit iriH <C Mw the 
leading behaviour in the on-shell case is given by [283] 



In ref. [284] this correction has been computed for off-shell W-boson production neglect- 
ing the s-channel vertex contributions, due to the dominance of the u exchange diagram 
at threshold. For an Higgs-boson mass of 60 GeV the correction amounts to about 

0. 8 — 0.9% and decreases for higher Higgs-boson mass values (being less than 0.1% for 
niH = 300 GeV). This effect translates into an uncertainty on the determination of Mw 
of 15 MeV from the LEP2 threshold run. 

Concerning QCD corrections, when at least two quarks are present in the final state, 

1. e. in semileptonic and fully hadronic channels, also the strong radiative corrections to 
the 4/ processes have to be considered. In the case of an inclusive set-up and in the 
CC03 approximation, the QCD corrections amount to use the QCD corrected values of 
the total and partial W^-boson widths in the cross section. These have been calculated in 
the literature [285] and can be expressed in factorized form by the following replacements: 



The impact of QCD corrections on the angular distribution of the decay products of a 
W boson and their application to the on-shell VF-pair production has been discussed 
in ref. [286]. The QCD corrections for the inclusive set up and in the double-resonant 



G^-G^/(l + A2/fj^) 



a Mw 



^Born • 




(^) The full one-loop clcctrowcak and QCD corrections to the decay width of the W boson have been 
calculated in ref. [281]. An Improved Born Approximation absorbing the bulk of the radiative corrections 
can be obtained by writing the lowest order width in terms of and M-^ [282]. 
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approximations can be used as estimates of the corrections for more realistic cases. For 
this reason the above replacements are referred to in the literature as naive QCD correc- 
tions. The quality of the approximation when cuts on the jet directions arc imposed has 
been discussed in ref. [287]. The effect of naive QCD corrections on the cross section of 
a typical 4/ process is shown in Fig. 53. 

Recently the exact 0{as) corrections to scmilcptonic and hadronic processes have 
been calculated [288] and compared with the naive approximation both in a fully inclu- 
sive situation and in the presence of ADLO-TH cuts [258]. The differences are below 
the foreseen experimental error at LEP2 for the scmilcptonic channels, while for fully 
hadronic final states the naive approach can give unreliable results. The analysis has 
been performed also at a cm. energy of 500 GeV, typical of future e+e" linear colliders, 
where the naive approximation deteriorates at the percent level for NLC/TH cuts [289]. 

Concerning at last QED corrections, in principle they should consist of the virtual 
photonic corrections and the real-photon bremsstrahlung in order to have an infrared 
safe result. A separation between initial- and final-state radiation respecting the U{1) 
electromagnetic gauge invariance, as in the case of LEPl s-channel processes, is not 
possible because the i/-exchange diagram involves a non-conserved charge flow in the 
initial state. The problem is circumvented if the treatment of electromagnetic radiation 
is restricted at the LL level. The first theoretical calculations taking into account such 
effects in 4/ productions can be found in refs. [290-296]. In this case the enhanced 
logarithmic terms coming from final-state radiation (FSR) cancel by virtue of the KLN 
theorem (when sufficiently inclusive cuts are considered) and the initial-state corrections 
reveal the presence of infra-red and collinear singularities, which are process-independent. 
Being phenomenologically relevant, these can be resummcd to all orders by means of one 
of the methods described in Appendix A. Given the steep slope of the born cross section 
in the threshold region, the numerical effect of ISR is to reduce the cross section by 
several percent (sec Fig. 53), depending on the cm. energy [254]. By means of various 
procedures (for instance the variation of the scale in the LL theoretical predictions) it 
is possible to state that the overall accuracy of the LL QED corrections is of the order 
of 1%. Another possible way of estimating non- leading QED effects is the one followed 
in ref. [297] , where a gauge invariant definition of ISR is given according to the current 
splitting technique. The electrically neutral neutrino of the f-channel diagram is split 
into two oppositely flowing charged Icptons. One of them is attributed to the initial state 
to build a continuous charge flow, while the other one is attributed to the flnal state. In 
this way the 0{a) ISR form factor for s-channel processes used at LEPl (see eq. (106)) 
is recovered, plus non-factorizing terms numerically small in the LEP2 energy range 

However, when aiming to study the impact of exclusive photon radiation on partic- 
ular observables, such as the photon spectrum or the differential angular distributions of 
flnal-state particles in the presence of realistic experimental cuts, some improvements on 
the LL approximation have to be introduced. For example, in ref. [279] the transverse 
degrees of freedom of emitted photons have been introduced by means of the pr depen- 
dent SF approach; in order to take into account the kinematical effects of the photonic 
radiation from all external charged particles. Alternatively, in ref. [298] the exclusive 
photonic contributions, including also multiple photon radiation from the W bosons, 
are implemented by means of a 0{a) YES exponentiation. In this formulation the 0{a.) 
corrections are implemented using the results for l^-pair production in the on-shell limit. 

In addition to the LL effects, a particularly sizeable photonic correction in the W- 
pair production region is the so-called Coulomb singularity (see Fig. 54). It originates 
from the electromagnetic interaction between the slowly moving W's. This effect is 




known since a long time [299] for stable particles to give a correction factor Sacoui of the 
form 

na 



where (3 = y/l — My^ / E"^ represents the velocity of the W bosons. The off-shellness and 
the width of the W bosons modify radically the Coulomb interaction because they act 
as an effective cut off on the range of the electromagnetic interaction. Only for high 
enough energies, such that the typical interaction time between the VF's is smaller than 
the W^-boson lifetime, the effect derived for stable particles is substantially unchanged. 
In terms of Feynman diagrams the Coulomb singularity arises from three- and four-point 
scalar functions contained in the graphs with a photon line connecting the two virtual 
W bosons (see Fig. 54). Being the coefficients of these scalar functions gauge invariant, 
they can be worked out to give the following correction factor to the lowest order cross 
section resulting from the double-resonant diagrams [300]: 



0(^Coul = <^Born ^ 



■ 2 f\f3M-P\ 
1 arctan ' 



with 



p = J^s2 - 2s(4 + ) + (ki - e 



4M^ 

pM = \jl —, M^=M^- iMwTw - is- 

where k+,k- are the four- momenta of the and W~, respectively. By performing 

the on-shcll limit on this formula the correction factor for stable particles is recovered. 
Higher-order effects due to the Coulomb correction are unimportant [300]. Numerically 
the Coulomb singularity amounts to a correction of the order of 6% at threshold and 
decreases smoothly to about 2% at ^/s = 190 GcV (sec Fig. 53). 

The QED corrections discussed above are phenomenologically relevant, gauge-inva- 
riant and include, besides the exponentiation of LL contributions, part of the full 0{a) 
clcctroweak corrections. The natural development is to consider the remaining non- 
leading 0{a) corrections. As stated previously, a step towards the calculation of the 
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bosonic contributions to the 0{a) radiative corrections to 4/ processes can be performed 
within the pole scheme [198,270,271], even if this gives rehable results only a few F^^'s 
above the VK-pair production threshold. In the threshold scan the required theoretical 
accuracy is, on the other hand, more relaxed with respect to the high energy points. 
Given the dominance of the double-resonant diagrams in the LEP2 energy region, a good 
approximation is to consider the radiative corrections in the double pole approximation. 
These can be divided in factorizable and non- f actor iz able corrections (see Fig. 55 for an 
example of non- factorizable virtual QED corrections). 




Fig. 55. - Examples of Feynman diagrams for non-factorizable virtual QED corrections. 

The former contain manifestly two resonant VF-boson propagators and can be dis- 
tinguished between corrections to VF-pair production and decay, while the latter do not. 
At present, the factorizable corrections to e+e~ — > 4/ can be taken into account by using 
the known 0{a) corrections for on-shcU M^-pair production and decay. This is a first 
reasonable approximation to the full one-loop calculation, since the neglected terms are 
of the order of aVw/T^Mw This approach has been already followed in ref. [298]. Re- 
cently, the effects of non-factorizable QED diagrams in the soft-photon limit have been 
investigated [301-305]. As a result of this calculations, the non-factorizable corrections 
vanish in the case of initial-final state interference, and in all cases when the integrations 
over both invariant masses of the virtual W bosons arc performed. Hence the W-boson 
production angle is insensitive to these corrections, so that the studies of anomalous triple 
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gauge-boson couplings at LEP2 are not aflFected by them. The case of the VF-boson in- 
variant mass distribution is different, since it receives a relative contribution of the order 
of 1%. The same corrections vanish for ZZ-mediated and ZiJ-mediated 4/ processes. 

5'2.4. Computational tools. As done in Sects. 2"4 and 3'5, for small-angle Bhabha 

scattering and physics, respectively, the basic features of the computer codes available 
for the calculation of observables of 4/ production processes are briefly summarized. The 
aim of the discussion is to show how the theoretical results described in Sects. 5'2.1 and 
5'2.3 are in practice implemented in computational tools used for the experimental analy- 
sis of Ty-pair production at LEP2. Since, as it will be discussed in Sect. 5'4, Higgs-boson 
physics at LEP requires the study of 4/ final states, some of the programs developed for 
M^-boson physics at LEP2 can provide also predictions for Higgs-boson production and 
decay. The goal of this discussion is to make an inventory of the theoretical approaches 
to the problem and of their realizations in the form of FORTRAN codes, rather then to give 
an exhaustive description of the programs. For a more detailed account of the programs, 
the reader is referred to refs. [258,306,307], and to the authors of the packages for the 
most recent developments. 

The most of the 4/ codes developed for LEP2 arc Monte Carlo programs. Indeed, 
since the phase-space integration over the 4/ final state involves seven dimensions, it can 
be more efficiently performed, especially in the presence of arbitrary kinematical cuts, 
by means of standard or adaptive Monte Carlo techniques. However, other numerical 
algorithms are available in the literature. For instance, semi-analytical programs, where 
a number of phase-space variables is integrated analytically, ending with two integrations 
over invariant masses (sec eq. (156)) to be; pc^rfornied numerically, and deterministic tools 
have been developed. In particular, semi-analytical codes are useful as benchmark for 
the Monte Carlo's and for fitting purposes. Independent of the particular numerical 
algorithm employed, variance reduction techniques (such as importance sampling) are 
demanded and used in practice in order to take care of the peaking behaviour of the 
integrand. 

Although different in many theoretical and technical aspects, most of the programs 
include the following common ingredients, necessary to match the LEP2 experimental 
precision (see the discussion in Sect. 5"2.3): 

• the full or the numerically relevant set of Feynman diagrams for the processes 
e+e" 4 fermions; 

• ISR in the coUinear approximation, usually in the form of QED SF's; 

• Coulomb correction; 

• QCD corrections, in the naive realization, to the total and partial W^-boson widths; 

• possibility of anomalous couplings (see Sect. 5 '3.2); 

• interface to hadronization packages. 

Some codes implement additionally QED FSR from fcrmionic legs. A few programs 
generate pr-carrying photons, at the level of ISR and/or FSR, while others are able 
to calculate the matrix element for the radiative process e+e~ — > 4/ -|- 7. Finally, 
some programs include one or more of the gauge restoration mechanisms discussed in 
Sect. 5-2.2. 
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ALPHA [265] — This program is based upon an original theoretical algorithm recently 
proposed in the literature [265] , and particularly powerful for the treatment of processes 
with a high number of final-state particles. Differently from all the other 4/ codes, ALPHA 
computes automatically the tree-level 2 — * 4 scattering amplitudes without making use 
of Feynman diagrams, but computing iteratively the saddle point of the path integral for 
given external momenta. All 4/ final states can be treated. It can be used to predict the 
rate for the radiative process e+e~ ^ 4/ -|- 7 [308]. 

CompHEP 3.0 [263] — It is a package of symbolical and numerical modules giving in 
output cross sections and distributions for processes with up to five particles in the final 
state, with a high level of automation. It uses the BASES&SPRING package for adaptive 
Monte Carlo integration and unweighted event generation. All 4/ and 4 fermions -|- 
1 photon final states can be treated. ISR is implemented using SF's in the coUinear 
approximation. 

ERATD [309] — It gives results for any 4/ final state, using the "E-vector" formalism for 

the calculation of the 2^4 hclicity amplitudes [259] . ISR in the coUincar approximation 
and QCD corrections are allowed. A qualifying feature of ERATO is the incorporation of 
all CP conserving anomalous couplings. The fermion-loop scheme as gauge restoration 
mechanism is active. The program is a Monte Carlo that can run as an event generator 
and as an integrator. It has been interfaced to JETSET/HERWIG. 

EXCALIBUR [310] — It is a Monte Carlo integrator that can give predictions for all possible 
4/ final states (excluding Higgs-boson exchange), with the possibility of selecting the 
contributions of subsets of diagrams and of particular spin configurations. The helicity 
amplitudes are computed using the Weyl-van der Waerden formalism [260]. ISR in the 
coUinear approximation. Coulomb and QCD corrections, as well as anomalous couplings, 
are available. 

GENTLE/4fan 3.0 [311] — It is a semi-analytical package designed to compute selected 

4/ production cross sections and invariant mass distributions for CC and NC mediated 
processes. SM Higgs-boson production in the NC case is included. All final states that 
do not contain identical particles, electrons or electron neutrinos can be studied. ISR in 
the coUinear approximation is available both in the radiator and SF approach. Moreover, 
non-universal QED corrections to some classes of diagrams can be included. Coulomb 
and QCD corrections, as well as anomalous couplings, are implemented. 
grc4f [312] It is a Monte Carlo generator for all 4/ final states automatically generated 
by the package GRACE [313]. All relevant radiative corrections and anomalous couplings 
are supplied. In particular, QED radiation is implemented by means of SF's for ISR, but 
a QED PS (QEDPS) [314] is also available both for ISR and FSR. Fermion masses can be 
kept nonzero everywhere. The numerical Monte Carlo integration and event generation 
is performed with the help of the package BASES&SPRING. Interface to hadronization is 
provided. 

KORALW 1 . 03 [315] — It is a Monte Carlo that can treat any 4/ final state via an in- 
terface to the GRACE library. All the phenomenologically relevant radiative corrections, 
anomalous couplings and interface to JETSET arc provided. In particular. ISR is for- 
mulated according to the YFS exclusive exponentiation, thereby including the effect of 
transverse photon momenta, while FSR is treated within the SF approach. Recently, 
0{a) non- leading QED corrections have been implemented [298] using the known results 
for on-shell W^-pair production. 

LEPWW [316] — This generator contains the CC03 and NC02 tree-level diagrams, supple- 
mented with 0{a) ISR and FSR treated according to the package PHOTOS [317]. QCD 
corrections, anomalous couplings and interface to JETSET are provided. The program 
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aims at a 1-2% precision. 

LPWW02 [318] — It is a Monte Carlo code containing the Feynman diagrams with two 
resonant VF's and Z's, initial- and final-state radiation with SF's in the coUinear approx- 
imation. Coulomb singularity and QCD corrections are incorporated. It is interfaced to 
JETSET. 

PYTHI A/ JETSET [319] — It is a general-purpose event generator for a multitude of pro- 
cesses in e^er , ep and pp collisions. It is especially designed for detailed modelling of 
hadronic final states, rather than precision electroweak studies. It includes the CC03 
diagrams for the Born matrix element and a hybrid SF-t-PS treatment of ISR. FSR 
is implemented according to a PS description. The Coulomb correction is available. 
Hadronization is built-in. 

WOPPER 1.5 [320] — It is a Monte Carlo generator for unweighted 4/ events. It includes 

the ceil set of diagrams, ISR according to a PS algorithm, with finite pr for photons 
generated according to the l/{p ■ k) pole. Coulomb and QCD corrections, anomalous 
couplings as well as interface to hadronization are provided. 

WPHACT 1 . [321] — This Monte Carlo program can compute all processes with four 
fermions in the final state. In particular, by virtue of the helicity formalism adopted 
that allows a semi-automatic calculation of the matrix elements [322,323], finite 6-quark 
mass effects are fully taken into account for NC processes. Higgs-boson production is 
included. ISR is included using SF in the coUinear approximation. The Coulomb term 
and anomalous couplings are incorporated. Interface to hadronization is present. 
WTO [324] — It is a deterministic code implenicniting the largest part of the 4/ Feynman 
diagrams (including Higgs-boson signal) obtained with the helicity method described 
in [261]. ISR using SF's in the coUinear approximation, anomalous couplings, QCD 
and Coulomb corrections are available. The fermion-loop scheme as gauge restoration 
mechanism is active. Some options for the estimate of the theoretical error are provided. 
Interface to hadronization is present. 

WWF 2 ■ 2 [325] — This Monte Carlo generator is a kind of merger of an explicit e^e~ — > 
4/ + 7 matrix element with SF's for the description of ISR beyond the LL approximation. 
However, non-leading corrections are implemented in an approximate way. The CC20 set 
of trec-lc!vel diagrams is included. FSR is implemented by using an explicit one-photon 
matrix element. All the other relevant radiative corrections, as well anomalous couplings 
and JETSET interface, are provided. 

WWGENPV 2 ■ 1/HIGGSPV [279] — They are Monte Carlo generators computing the largest 
part of 4/ tree-level matrix elements (including Higgs-boson production) derived ac- 
cording to the helicity formalism of ref. [261]. ISR and FSR are treated within the 
SF approach, including pt/pl effects. Coulomb and QCD corrections are implemented. 
Anomalous couplings and interface to JETSET are also available. 

Before the start of the LEP2 operations, an extensive comparison of all the Monte 
Carlo and semi-analytical programs available for the analysis of 4/ processes was carried 
out in the context of the working group "Event generators for WW physics" of the CERN 
workshop "Physics at LEP2" [258]. The predictions for both total cross sections and more 
exclusive observables of primary importance for the measurement of the properties of 
the W boson (such as VF-boson production angle, invariant masses, radiative energy loss 
etc.) were compared in detail. "Tuned" comparisons generally showed very satisfactory 
agreement between the codes (at the level of 0.1%), pointing out a correct implementation 
of the advertised features and a high technical precision. Other highlights of the study 
can be summarized as follows 
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• the contribution of background processes is in general not negligible and therefore 
dedicated 4/ codes are more suitable for precision measurements of the W^-boson 
properties than general-purpose programs; 

• for several observables, the effect of photon pt is important on both ISR and FSR; 

• the available estimates of the present theoretical error show that it is not much 
smaller than the expected experimental uncertainty; as a consequence, a more 
reliable estimate of the complete one-loop corrections to e+e~ — > 4/ could give 
some relevant piece of information. 

5'3. The W-hoson Mass and Anomalous Couplings. - As discussed in Sect. 5'2, 

the basic aims of the study of W^-pair production in 4/ final states at LEP2 arc the 
precise direct measurement of the mass of the W boson, Mw, and the detection of 
possible anomalies in the jWW and ZWW vertices {anomalous couplings). All the 
decay channels of VT-pair production are used in the experimental analysis: the hadronic 
channel WW qqqq (45.6% decay fraction), the semi-leptonic channel WW Ivqq 
(43.8%) and the leptonic channel WW lulu (10.6%). 

5'3.1. The W^-boson mass. Two different methods are employed to measure the mass 
of the W boson at LEP2 [248]: 

• the threshold method; 

• the direct reconstruction method. 

The former determination exploits the strong sensitivity of the W^W~ production cross 
section to the W^-boson mass near the nominal threshold; in the latter method, the 
Breit-Wigner resonant shape of the invariant mass distribution of the W^-boson decay 
products is directly reconstructed, thus yielding M\y. The statistically most precise 
determination comes from the direct reconstruction method, since most of the LEP2 
data is and will be collected at energies well above the threshold. Furthermore, very 
different systematic errors affect the two strategics, that therefore can be considered as 
complementary methods able to provide an internal cross-check on the M^-boson mass 
measurement at LEP2. 

As can be seen from Fig. 50, the cross section for VK-pair production increases very 
rapidly near the nominal kinematic threshold at -^i = 2Mw^ although the finite VF-boson 
width and ISR significantly smear out the sharp rise of the (unphysical) cross section for 
producing two on-shcU W^s (see Fig. 53). This means that for a given ^/s near threshold 
the H^-pair production cross section is very sensitive to the M^-boson mass, so that a 
measurement of the cross section in this energy region directly yields a measurement 
of Mw- In the threshold determination, the error on the W-boson mass due to signal 
statistics is given by 

1 

where aww is the total W^-pair production cross section, £ww is the overall signal 
efficiency and L is the integrated luminosity. By studying the sensitivity factor given by 
sJ(Jww \dM\Y / daww\ as a function of y/s — 2Mw, one finds that there is a minimum 
located at 



(159) AMvr = 

acTww 



(160) 



~ 2Mw + 0.5 GeV, 
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corresponding to around 161 GeV as the cm. energy providing the optimal sensitivity to 
Mw (see Fig. 56). Since the statistical uncertainty is the most relevant source of error for 
the threshold method, the optimal strategy for data taking consists therefore in operating 
at ~ 161 GeV in order to minimize the statistical error. The systematic errors on 
Mw in the threshold method are due to the uncertainties in the PF-pair production 
cross section: the luminosity error, unknown corrections to the theoretical cross section, 
as those discussed in Sect. 5'2.3, the uncertainty in the beam energy, the background 
subtraction and others. However, as already stressed, the overall error in the threshold 
measurement is largely dominated by the limited statistics. The present LEP average 
for Mw from the WW production cross section at ~ 161 GeV is [8] 

(161) Mw = 80.40 ±0.22GeV. 

The direct reconstruction method consists in reconstructing the peak in the invariant 
mass distributions of the VF-boson decay products, by using either the semi-leptonic 
decay mode (W^W li'qq) or the hadronic one {WW — > qqqq). The measurements 
are carried out at energies above the threshold (y^ > 170 GeV) in order to reduce the 
statistical error thanks to the larger WW cross section. The accuracy of the jet energy 
measurement is too poor to allow a precision measurement of Mw- Therefore, kinematic 
fit techniques, using the constraints of energy and momentum conservation, together with 
the equality of the two VF-boson masses in an event, are used by the LEP Collaborations 
to improve the mass resolution [248]. This implies that a good knowledge of the cm. 
energy is essential. As already pointed out, systematics affecting this method are largely 
independent of those present in the threshold method. In particular, the qqqq channel 
turns out to be more problematic than the luqq one: indeed, besides a larger background 
contamination, the phenomena of colour reconnection and Bose-EinsteAn correlations (sec 
below) introduce additional systematic uncertainties that are absent in the semi-leptonic 
channel. Systematic errors that are common to both channels are 

• error from the LEP beam energy; 

• errors from the theoretical description (ISR, background, 4/ diagrams and frag- 
mentation); 

• errors from the detector. 

In order to estimate the uncertainties due to 4/ diagrams and ISR, the theoretical in- 
gredients discussed in Sects. 5'2.1 and 5'2.3, and the corresponding computational tools 
based on them, are very useful tools. Indeed, the effect of 4/ subprocesses other than 
the ones belonging to the CC03 class, induces an uncertainty of 20 MeV and 25 MeV on 
Mw for the qqqq and Ivqq channel, respectively, and is evaluated by using MC samples 
generated with one or more of the 4/ generators described in Sect. 5'2.4. Moreover, pre- 
cision calculations of the radiative energy loss (see eq. (158)) in 4/ production at LEP2 
are mandatory. Actually, the relatively large average energy (< >~ 1.1, 2.1, 3.2 GeV 
at = 175, 190, 205 GeV, respectively) carried away by the radiated photons leads to a 
significant mass-shift if it is not taken into account in rescaling or constraining the ener- 
gies of the final-state VF-boson decay products to the beam energy. Indeed, the average 
mass-shift caused by ISR is of the order oi < E^ > Mwl\fs, i-e. about 500 MeV at 
\/s = 175 GeV. However, since a fit to the mass distribution gives more weight to the 
peak, the actual mass shift at y/s = 175 GeV is about 200 MeV, anyway important in 
view of the envisaged experimental precision on Mw- Since, as discussed in Sect. 5'2.3, 



PRECISION PHYSICS AT LEP 



117 



the exact treatment of 0{a) corrections to off-shell VF-pair production is not yet avail- 
able, the theoretical predictions for < >, including the presently under control LL 
and Coulomb corrections, arc affected by an intrinsic uncertainty, that translates into 
an error of 10-15 MeV on the VF-boson mass. This theoretical uncertainty is due to 
non-leading photonic terms, for instance associated with the radiation of photons off the 
intermediate W bosons, and should be reduced soon in the light of the recent progress 
in the calculation of non-leading QED corrections discussed in Sect. 5'2.3. 

Colour reconnection and Bose-Einstein effects [248] , that are systematics specific to 
the qqqq channel, arc interconnection phenomena that may obscure the separate iden- 
tities of the two W bosons, so that the 4/ final state may no longer be considered as 
consisting of two separate VF-boson decays. Indeed, in the LEP2 energy regime the aver- 
age distance between the and decay vertices is smaller than 0.1 fm, i.e. less than 
a typical hadronic size. Therefore the fragmentation of the and the W~ may not be 
independent and may seriously affect the mass reconstruction. The colour-reconnection 
problem arises from the fact that the T4^-boson pair, decaying into quark- antiquark pairs 
qq and QQ, respectively, can either fragment into two strings stretched between qq and 
QQ, or into two strings stretched between qQ and Qq. Owing to the colour charges, 
the probability and the properties of the two possible configurations are very differ- 
ent. Of the two phases characterizing the fragmentation of the initial quark-antiquark 
pairs into hadrons, i.e. perturbative parton cascade and non-perturbative hadronization, 
the perturbative QCD interconnection effects have been shown not to be very impor- 
tant [326,327], whereas the infiuence of non-perturbative fragmentation can be only 
estimated by comparing difii'erent models. The interesting aspect is that quite different 
approaches give imcertainties of the same order, i.e. about 100 MeV, conservatively. Also 
Bose-Einstein correlations are due to the overlapping of the and W~ hadronization 
regions: consequently, low-momentum bosons coming from different Ws can experience 
coherence effects. Their influence on the VF-boson mass-shift can be estimated by using 
theoretical models. However, it is worth noticing that experimental studies at LEP2 
show no evidence yet for Bose-Einstein correlations and that recent theoretical investi- 
gations suggest that they may not be a serious systematic (sec for instance [328]). The 
present experimental attitude consists in trying to measure interconnection effects in 
data through detailed comparison of qqqq and Ivqq channels. 



The present LEP average for Mw from the direct reconstruction method at ^/s = 

172 GeV is [8] 



It is worth noticing that good consistency between experiments and between qqqq and 
Ivqq channels is observed. 

By combining the i/s = 161 GeV (threshold) and the ^/s = 172 GeV (direct recon- 
struction) results, the present Mw value from LEP2 is [8] 



The LEP2 result for Mw can be combined with the measurement by CDF/DO at the 
TEVATRON, in order to compile a world average. The M^-boson mass measurements at 
the TEVATRON are based on fits to the transverse mass distribution of li' produced in 
the process qq ^ W lv{l — e^ji). By combining the CDF/DO results with the UA2 
measurement, the average Mw from hadron machines is [8] 



(162) 



Mw = 80.53 ±0.18 GeV. 



(163) 



Mw = 80.48 ±0.14 GeV. 



(164) 



Mw = 80.41 ±0.09 GeV. 
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The present world average Myy from LEP2 and pp colliders is [8] 
(165) Mw = 80.43 ± 0.08 GeV . 

By the end of LEP2 and TEVATRON run 2, both experiments expect to achieve an 
error on Mw of around 0.03-0.04 GeV, comparable with the present error from indirect 
determination by means of precision data. 



W"^W~ cross section at LEP 

PRELIMINARY 




Q I I I I I I I I I I I I I I I I I 

160 170 180 190 

Vs [GeV] 

Fig. 56. - The W-pair production cross section (from [8]). 



5 '3. 2. Anomalous couplings. In addition to the VF-boson mass, the second impor- 
tant piece of information that VF-pair production at LEP2 can provide is the structure 
of the triple gauge-boson couplings (TGC). Contrary to LEPl physics, where TGC only 
enter through loop corrections to two-fermion production, at LEP2 these couplings are 
responsible of the behaviour of the tree-level W^-pair cross section. Since the specific 
form for the TGC in the SM is a consequence of the Yang-Mills character of the theory, 
the study of VF-pair production at LEP2 allows to test directly the non-abelian nature 
of the SU{2) X U{1) gauge theory. Indeed, in spite of the precision tests of the elec- 
troweak interaction at LEPl and SLC, the non-abelian sector of the SM remains poorly 
measured to date, even if the very existence of non-abelian coupling can be considered 
experimentally proved by the recent LEP2 data (see Fig. 56). 

WW Z and WWj couplings different from the form predicted by the SM are called 
anomalous couplings. The latter will in general lead to a bad high-energy behaviour of 
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cross sections violating unitarity. However, since at LEP2 the cm. energy is too low to 
point out such effects, the presence of anomalous couplings can be eventually established 
through the study of the angular distributions of the W bosons and their decay prod- 
ucts. As for the W-boson mass measurement, the control of radiative corrections and 
background contributions is necessarily required, since they induce deviations from the 
tree-level angular distributions. 

Possible anomalies in the WWV {V = Z ov 7) vertex can be parameterized in terms 
of a purely phenomenological effective Lagrangian [329,330] 



icYfY'' = 9WWV [gYv^ {w-,w+'' - w+^w--) + kvW+W-V^'' + 



^^^^^ 



that provides the most general Lorentz invariant WWV vertex, observable in processes 
where the vector bosons couple to effectively massless fermions. The definitions in 
eq. (166) read as follows 

gwwi = e gwwz = e cot 

(167) W^^ = d^W^- d^W^ V^^ = - d^V^ . 

Therefore, 2x7 independent parameters are needed to describe the most general Lorentz 
invariant WWV vertex. Within the SM and at the tree level the couplings in eq. (166) 
are given hy = = kz = k-y = I, with all other couplings vanishing. It is worth 
noticing that gY , ky and Xy conserve C and P separately, while gX violates C and P 
but conserves CP. Furthermore gy, ky and Ay parameterize a possible CP violation in 
the bosonic sector. In particular, the C and P conserving parameters in C^j^^ can be 
linked to the static e.m. moments of the W^ boson as follows [331] 

charge Qw = egj , 

2Mw 

(168) electric quadrupole moment qw = — , ,2 (^7 ~ ^1) ■ 



magnetic dipole moment fj,w = 7777 — {gj + k^ + X~^) , 



In practice it is impossible to set limits on all the above couplings, so that a number 
of assumptions have to be made in order to reduce the number of parameters. In the 
literature, two different theoretical strategies are followed. On the one side, the number 
of parameters can be reduced by advocating symmetry arguments [330,332]. On the 
other hand, the dimensional analysis of the operators involved in the parameterization 
allows to establish a hierarchy of the operators themselves according to a scaling law of 
the kind {E/KY-^, where E is the energy of interest, A is the typical energy scale of non- 
standard physics and d is the dimension of the operators. The latter procedure allows 
to single out a reduced set of parameters able to produce phenomenologically relevant 
effects at the energy scale of interest [333] . 
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For instance, if C, P and electromagnetic gauge invariance are imposed, the number 
of parameters reduces to five, namely \-y,Xz,k-y,kz and gf . Furthermore, since the 
effective Lagrangian of cq. (166) contains as a particular case the triple gauge-bosons 
Lagrangian of the electroweak theory, it is possible and convenient to introduce deviations 
from the SM predictions for the five parameters as 

(169) Agf=.gf-1, A^,z = fc7,z-1, A,, \z . 

A set of parameters widely used by the LEP Collaborations is given by the following 
three combinations, which do not aff'ect the tree-level gauge boson propagators and are 
not yet constrained by low energy and peak data [249] 

Afc^, - Agf cos^ '&W = otB^,, 

(170) \z = \ = aw, 

with the constraint Afcz = Agif — Ak-^t&v^'dw The above combinations are all zero 
according to the SM. 

The experimental strategy followed at LEP2 for the determination of anomalous cou- 
plings relies upon the measurement, in 4/ production processes, of the differential cross 
sections d"cr = d"f7(9, dl^.^pl^^a.) with n— 1, 3, 5, where 6 is the polar production angle 
of the W~ boson, d*^ and ip*^ are the decay angles in W-boson rest frames. Among 
the VK-pair decay channels, the strongest information comes from the semi-leptonic chan- 
nel because it allows unambiguous charge assignment. Also single VF-boson and single 7 
production processes, that are sensitive to the non-abelian WW-f coupling are studied 
by the LEP Collaborations, even if they are less powerful than the WW 4/ reactions. 
Present combined results from the four LEP experiments are [8] 

aw4. = 0.02l°;J^ - 0.28 < aw^, < 0.33 (95% CL limit) 

aw =0.15+5];^^ - 0.37 < aw < 0.68 (95% CL limit) 

(171) aB<i, = OA^to lr - 0.8K < 1-50 (95% CL limit) 

Therefore no discrepancy with the SM is seen in the data. Moreover, an important result 
of the anomalous couplings analysis is that the existence of the ZWW vertex is firmly 
established, namely gwwz = is excluded at > 95% CL (see Fig. 56). 

5'4. Higgs-hoson Searches. As shown in Sect. 4'2.2, the electroweak precision data 
can be conveniently used to obtain indirect determinations of the Higgs-boson mass 
through radiative corrections to the physical observables. Although the weak dependence 
of the theoretical predictions on the Higgs-boson mass as well as the their intrinsic 
uncertainties prevent to derive stringent predictions for uih, there is a slight preference 
in the precision measurements for a relatively light SM Higgs boson, with an upper 
limit mn < 420 GeV @ 95% CL [7]. Future more precise direct measurements of the top- 
quark mass at the TEVATRON and of the VF-boson mass at LEP2 and the TEVATRON, 
associated with recent progress in the calculation of two- loop electroweak graphs, should 
improve the present constraints on the mass of this yet elusive particle. Indeed, the 
search for the Higgs boson certainly constitutes one of the main tasks of present-day 
experiments in high-energy physics. The discovery of this neutral scalar particle would 
establish the validity of the electroweak theory in its standard formulation, and shed 
light on the mechanism of electroweak symmetry breaking and mass generation. 
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Fig. 57. - Tree-level Feynman diagram for Higgs-boson production in the channel e^e 
fi^H'bb (Higgs-strahlung). 



By searching directly for the SM Higgs boson, the LEP experiments have set a 
lower bound on the Higgs-boson mass, namely ttih > 77GeV @ 95% CL (four LEP 
experiments combined) [182, 183]. This lower limit comes from negative searches of 
Higgs-boson production in electron-positron collisions around and above the peak. 
The dominant production mechanism for the SM Higgs boson in the LEP energy range 
is the s-channel Higgs-strahlung process e+e~ ZH, where the intial states annihilate 
into a virtual Z^ boson that converts into a Z° and a Higgs boson. Other production 
mechanisms, i.e. the fusion processes, where the Higgs boson is formed in WW.ZZ t- 
channel collisions, with the W, Z's radiated off the incoming electron and positrons, have 
smaller cross section at LEP energies. However, it is worth keeping in mind that for 
fe^e and e+e" in the final state, the two amplitudes for the Higgs-strahlung and fusion 
process both contribute and interfere. 




Fig. 58. - Tree-level Higgs-boson production via WW fusion. 

In view of the LEP2 physics potential, the reachable Higgs-boson mass limit is 
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roughly given by ~ ^/s — 100 GeV. Even assuming the highest cm. energy in the 
LEP2 operation = 205 GeV), the discovery potential of the LEP2 experiments will 
not exceed a Higgs-boson mass of around 100 GcV. For a SM Higgs boson with a mass 
in the range between 77 (present lower limit) and 100 GeV, the dominant decay mode 
is by far H — > hh. Branching ratios for the decays H — > t~^t~ ,cc,gg are suppressed by 
an order of magnitude or more. The corresponding total decay width is predicted in the 
electroweak theory to be very narrow, i.e. less than a few MeV for niH less than 100 GeV. 
As a consequence of the fact that the Higgs boson predominantly decays into bb pairs 
when TTiH < 100 GcV, the observed final states in Higgs-boson searches at LEP consists 
of four fermions. In the case of the Higgs-boson signal these 4/ final states are mainly 
achieved through the Higgs-strahlung process e+e~ — > ZH, with the subsequent decays 
Z f f and H bb. In general, whenever excluding Higgs-boson decay channels 
other than H ^ bb (such as H ^ t~^t~ that is actually considered in Higgs-boson 
searches at LEP), the interesting final states for the physics of the Higgs particle at LEP 
arc bb^i'^ fi~ , bbe'^e~ (leptonic channels), bbvD (missing energy channel) and bbqq (four-jet 
channel). Therefore, Higgs-boson physics at LEP2 requires, as for VF-pair production, the 
calculation of 2 — > 4 fermions scattering amplitudes, including signal and backgrounds in 
the SM. The Feynman diagram for the Higgs-boson signal in the channel e+e~ ii^ ^i~bb 
is depicted in Fig. 57. The additional Higgs-boson production mechanism in the case of 
Uf.Pebb final state can be found in Fig. 58; an analogous diagram with ZZ fusion is 
present for the e^e~bb final state. For the channel fi'^ii~bb, that provides the cleanest 
event signature, background events are generated, among others, by double vector-boson 
production e+e~ — > ZZ, Z'y and 77 with the virtual Z° and 7 decaying to and 
bb pairs (see Fig. 59). Actually, just to give an idea of the degree of complexity of the 
full calculation of 2 ^ 4 amplitudes for Higgs-boson physics, 25 diagrams are present in 
the ji^fjrhb channel, 50 diagrams in e~^e~bb and 68 in the bbbb channel (unitary gauge). 
Full 4/ calculations for Higgs-boson searches are certainly less demanded with respect 
to the case of precision studies of W^W~ production. Indeed, given the small Higgs- 
boson production cross sections (of the order of tens of fb, two-three orders of magnitude 
smaller than the e^e~ W'^W~ cross section) and the low statistics, precise predictions 
are not strictly necessary for discovery physics. However, the theoretical and technical 
expertise developed for 4/ production in W^-boson physics and described in the previous 
sections allows to apply the corresponding "machinery" to a comprehensive analysis of 
Higgs-boson production as well. The main motivations for complete 4/ calculations for 
Higgs-boson physics at LEP is that no approximations are introduced in the tree-level 
matrix element, differential distributions for the final-state products arc directly available 
and the background can be precisely evaluated. Therefore, a full calculation allows to test 
the degree of validity of the usual approximation of computing production cross section 
X branching ratios and of summing squared amplitudes incoherently, thus quantifying 
off-shellness and interference effects. Furthermore, once the Higgs boson is discovered, 
complete 4/ calculations and relative generators could be of great help in the extraction 
of the Higgs-boson properties (mass, spin, parity, etc.) that require a good control of 
exclusive distributions. 

Complete 4/ calculations for Higgs-boson physics at LEP have been carried out in 
the last few years by several groups [258,306,334-338]. The typical calculation strat- 
egy consists in computing all tree-level diagrams for a given channel and including the 
best available set of radiative corrections, i.e. running electroweak couplings, initial-state 
QED radiation, naive QCD final-state corrections and running heavy-quark masses. This 
framework is implemented, although with some minor differences, in 4/ programs such as 
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Fig. 59. - Tree-level Feynman diagrams for the main background processes in the channel 
e~ — > fi^ n~bb. 



compHEP, GENTLE/4f an, HIGGSPV, HZGEN [335], WPHACT and WTD, that have been described 
in Sect. 5'2.4. The programs compHEP, GENTLE/4f an, HZGEN and WPHACT can evaluate 
matrix elements taking into account the finite b quark mass, while HIGGSPV and WTO 
consider masslcss final-state fcrmions. WPHACT includes also the contribution of SUSY 
neutral Higgs-boson production. Some of the codes provide unweighted events with the 
four-momenta of all final-state particles, that can be used to process the event through 
the detector and to apply cuts for data analysis. The above 4/ generators arc described 
in more detail in [306], where also reference to other programs for Higgs-boson physics 
at LEP, such as PHYTIA and HZHA, based on production times decay approximation and 
with an incomplete treatment of the SM amplitudes, can be found. In ref. [306] some 
comparisons between the predictions of different programs are shown for a sample of 4/ 
final states, with the conclusions that the agreement between the dedicated 4/ generators 
is systematically at the level of 1% or better, pointing out the high technical precision 
achieved, and that there are situations, such as the bbv^Ve final state, where the difference 
between 4/ codes and standard computational tools for Higgs-boson searches at LEP, 
such as HZHA, not including the full set of SM diagrams, can reach 20-40%, depending on 
the cm. energy and the Higgs-boson mass value. Thanks to the availability of complete 
4/ calculations, some topics previously not investigated but of interest for Higgs-boson 
searches, such as the amount of signal-background interference and the effect of finite 
fo-quark mass corrections, received the due attention. In the case of massless final-state 
fermions, the interference between signal and background diagrams is exactly vanishing, 
because massless fcrmions arc coupled to spin-vectors in 7, Z decays and to spin-scalar 
in Higgs-boson decay, so that the two amplitudes do not interfere as a consequence of the 
different helicity pattern. If the mass of the h quark is kept different from zero in the full 
matrix element, then a finite interference does develop. However, it was found that this 
effect never exceeds the percent level. Indeed, since for niH < 100 GeV the Higgs-boson 
width is of the order of a few MeV, it can be expected a priori that the non-vanishing 
signal-background interference is highly suppressed at LEP, as demonstrated by explicit 
calculations. Therefore, signal-background interference may be neglected in Higgs-boson 
search experiments at LEP. Other conclusions that comes from these analyses is that 
ISR is large, varying between 10-20% on the total cross section, and must be included in 
Higgs-boson event generators, especially in the light of the distortions induced on exclu- 
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sive distributions that are relevant for the determination of the properties of the Higgs 
particle. 
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Fig. 60. - The cross section for Higgs-boson production in the channels e — > bbn jjT and 
e^e~ —> hhvV. The dotted and dash-dotted lines show tree-level predictions; the solid lines 
include ISR. The background cross section is also shown. Numerical results by HIGGSPV [279]. 



The predictions obtained with the help of the 4/ generator HIGGSPV for the Higgs- 
boson production cross section in the channels e+e~ — > bbfji^ jjT and e+e" — > hbv9 can 
be seen in Fig. 60, for two representative values of the Higgs-boson mass. The effect 
of backgrounds and ISR are shown as a function of the cm. energy up to 0.5 TeV. It 
is worth observing that in the LEP2 energy range the size of the total cross sections 
in both channels can give by itself clear evidence of Higgs-boson production. In the 
case of the hhvv final state the rise with the cm. energy of the cross section is due 
to the logarithmic enhancement introduced by the f-channel WW fusion processes that 
arc absent in the bbi.j,~^ fj~ final state. As said above, in order to extract information 
on the quantum numbers assignment of the Higgs particle, the angular distributions are 
known to be the most sensitive observables. An example of such quantities is given in 
Fig. 61, showing the distribution, obtained by means of HIGGSPV, of the 6-quark scattering 
angle in the laboratory frame, at ^/s = 192 GeV and for three Higgs-boson mass values 
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and background alone (e+e" bhvD channel). As can be seen, in the presence of 
the Higgs-boson signal as dominant contribution to the hbvv cross section, a clearly 
isotropic spin zero behaviour is present, whereas it disappears whenever considering 
the background only. Since this differential distribution can be meaningfully and more 
extensively analyzed by means of complete 4/ calculations, the above example should 
clarify the usefulness of dedicated precision tools for the measurement of the Higgs-boson 
properties. 




Fig. 61. - The 6-quark scattering angle distribution in the laboratory frame at = 192 GeV 
for the process e^e~ — » bbi'i>. Numerical results by HIGGSPV [279]. 



To conclude the discussion on Higgs-boson searches in electron-positron collisions, it 

is worth noticing that, if LEP does not discover the Higgs boson and its mass lies in the 
intermediate range (140 GeV < uih < "^Mz), the search for the Higgs particle at future 
e+e~ colliders will require the analysis of six-fermion production processes, due to the 
decay chain H W~^W~ 4/. Also in this case, in view of the expected experimental 
precision and particularly for a precise determination of the Higgs-boson properties, 
complete six-fermion calculations of Higgs-boson production and relative computational 
tools should be desirable. Actually, the first results along this direction very recently 
appeared in the literature [339,340]. 



126 



G. MONTAGNA, O. NICROSINI and F. PICCININI 



6. — Conclusions 

Since 1989 the four experiments ALEPH, DELPHI, L3 and OPAL at LEP and 

the SLD Collaboration at SLC collected a large amount of precision data by studying 
the processes of two-fermion production in electron-positron collisions at centre of mass 
energies close to the mass of the Z° boson. About 16 millions of Z'^ events have been 
recorded. The LEP phase around the resonance (LEPl) was terminated in 1995, 
while SLC is continuing in data taking on the Z^ peak. During 1996, the LEP energy 
was increased in order to allow the production of W-boson pairs for the first time in 
high-energy er collisions. The LEP2 phase is in progress and will continue imtil 1999. 
Meanwhile, almost final results of the analysis of the LEPl/SLC precision data have 
become available, although the data analyses, including the final LEP energy calibrations, 
are not yet finished. 

Looking at the high accuracy reached by the most recent measurements of the elec- 
troweak observables, the program of precision physics at LEP/SLC can be considered as 
a complete success, thanks to a combined experimental and theoretical effort. Actually, 
besides the very successful performance of the LEP machine and important experimen- 
tal/technological achievements, the effort undertaken on the theoretical side in precision 
calculations for Z^ physics significantly contributed to the progress of precision tests of 
the Standard Model (SM). 

A key ingredient of the striking success of precision physics at LEP/SLC is undoubt- 
edly the high precision determination of the machine luminosity. This has been possible 
thanks to the experimental achievements on luminometers, and to the high-precision 
calculations of the small-angle Bhabha process, today accurate at the level of 0.1%. 

The extremely accurate theoretical predictions for observables of e+e~ collisions 
into two-fermion final states at large scattering angles provided the necessary theoretical 
background for the precise determination of the electroweak parameters performed by 
the high-precision experiments. 

Concerning the comparison between theory and experiment, the electroweak pre- 
cision data confirm the validity of the SM with an impressive accuracy. The standard 
theoretical framework is able to accomodate all experimental facts, providing predictions 
in agreement with precision measurements. Just to mention some of the major achieve- 
ments, the number of light neutrino species has been unambiguously determined to be 
equal to 3; the Z-boson mass is at present known with a relative precision of 2 x 10~^, 
comparable with the precision of the muon decay constant G^; lepton universality has 
been tested with an imprecedented precision; the stringent constraints on the neutral 
current couplings allow the determination of the weak effective mixing angle with an 
absolute error at the level of 3 x 10~^; the W-boson mass is indirectly determined with 
a relative error below 0.1%. 

By virtue of their high accuracy, the precision measurements show a clear evidence 
of pure weak radiative corrections and they can be used to infer valuable information 
about the fundamental parameters of the SM. Actually, the measurements are sensitive, 
via the virtual effects of radiative corrections, to the top-quark mass mj, the Higgs-boson 
mass tuh and the strong coupling constant as- The indirect constraints obtained from 
precision data on the mass of the particles not energetically accessible clearly illustrate 
the role of quantum loops in raising precision physics to the level of discovery physics. 

From all the available data (excluding the direct determination of mt and Mw from 
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the TEVATRON and LEP2) the best current estimate of the top-quark mass is [7, 8] 

mt = 157tl° GeV. 

This indirect determination is in beautiful agreement with the measurement of the top- 
quark mass from direct production at the TEVATRON [168], mt = 175.6 ± 5.5 GeV, 
illustrating the constraining power of the precision measurements. 

The value of as{Mz) derived from a fit to all electroweak data is [7,8] 

as{Mz) = 0.120 ±0.003, 

in very good agreement with the world average and of similar precision. This indirect 
determination, when associated with the ag measurements from other processes and at 
different energy scales, significantly proves the running of as predicted by the non- 
abelian structure of QCD. 

The Higgs-boson is the still missing block of the SM. Although the constraints in- 
ferred on the mass of the Higgs boson from precision data are not conclusive, as a conse- 
quence of the weak logarithmic dependence of radiative corrections on mn , electroweak 
measurements imply an indicative mass window, thus providing valuable information in 
view of the planned searches at the LHC and future e+e" colliders. A fit to all data 
gives [7,8] 

niH = 115t^f GeV, 

mii- < 420 GeV (95% CL), 

where also the theoretical uncertainty due to missing higher-order corrections is taken 
into account in the one-sided 95% CL upper limit. A general caveat is in order here. The 
upper limit quoted above depends heavily on the presence of the left-right asymmetry 
from SLC in the data set. Therefore, the only reliable conclusion is that a very heavy 
Higgs boson, that would be rather problematic for future direct searches, is excluded 
by precision data. Future improved measurements of the W-boson mass at LEP2 and 
the TEVATRON and of the iop-quark mass at the TEVATRON, together with recent 
progress in the calculation of higher-order radiative corrections, are expected to sharpen 
the existing constraints on m/f . 

The first experimental results of LEP2 are available. As in the LEPl case, the 
theoretical tools, noticeably the accurate predictions concerning reactions of the kind 
e'^e~ 4-fcrmion final states, provide the necessary theoretical scheme for the data- 
theory comparison. The recent LEP2 measurements of the W pair production cross sec- 
tion already show clear evidence for the existence of the gauge bosons self-interactions, 
thus testing a crucial prediction of the non-abclian structure of the theory. Future mea- 
surements of the VF-boson mass are expected to reach an accuracy comparable with the 
one obtained in the direct determination performed at the TEVATRON and the indirect 
determination at LEPl/SLC as calculated via radiative corrections. 

Looking beyond the SM, the possible scenarios of new physics are strongly con- 
strained by the delicate agreement between the SM and the precision data. Since no 
significant anomalies arc present, the natural candidates are represented by models that 
do not modify the structure of the SM significantly. For this reason, theories with funda- 
mental Higgs particles, such as supersymmetry, are generally considered more favourite 
than models implying composite Higgs bosons and new strong interactions. Given the 
large value of the iop-quark mass, a fundamental Higgs-boson with a relatively low mass. 
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as inferred by precision electroweak data, can be nicely accommodated in a supersymmet- 
ric theory (such as the Minimal Supersymmetric Standard Model) broken at some high 
energy scale. Supcrsymmetry provides also a viable scenario to realize the unification of 
the coupling constants in a grand unified theory (GUT). Actually, the determination of 
sin^ "ffeff and aa{Mz) as obtained from precision data is compatible with gauge coupling 
unification at a large energy scale, provided the mass spectrum of the GUT includes 
supersymmetric particles. 

Although new elementary particles have not been directly discovered by experiments 
at LEP and SLC, it can be said that precision physics provided important pieces of 
information to our present knowledge in particle physics and contributed significantly to 
obtain hints on how new physics beyond the SM is realized in Nature. 
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A. — Universal Photonic Corrections 

A common set of radiative corrections to s-channel annihilation and f-channel scat- 
tering processes is represented by QED radiative corrections. In the so called leading 
logarithmic (LL) approximation, they are dominated by long-distance contributions, and 
hence are process independent. Going beyond the LL approximation, it is necessary 
to take into account process-dependent corrections, which are no more universal and 
can be computed by means of standard diagrammatic techniques. In the following, an 
overview of the most popular algorithms developed for the computation of QED correc- 
tions in the LL approximation, namely the Structure Function (SF) method, the Parton 
Shower (PS) method and Yennie-Frautschi-Suura (YFS) exclusive exponentiation, will be 
given. (*) Just for definiteness, only s-channel processes will be considered, the general- 
ization to t-channcl processes being almost straightforward. It is worth noticing that for 
two of these methods, namely the SF method and YFS exclusive exponentiation, proper 
procedures for matching the all-orders LL results with finite-order exact diagrammatic 
results have been developed. 



(*) It is worth noticing that at least two additional frameworks for the computation of QED corrections 
in the LL approximation have been also developed, namely the Coherent States approach [341] and the 
Unitary Method [342]. The interested reader is referred to the original literature. 
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A'l. The Structure Function Method. - Let us consider the annihilation process 
e~e+ — > X, where X is some given final state, and let (Jo(s) be its lowest order cross sec- 
tion, possibly including all the short-distance process dependent corrections, as those dis- 
cussed in Sect. 3. Initial-state (IS) QED radiative corrections can be described according 
to the following picture. Before arriving at the annihilation point, the incoming electron 
(positron) of four-momentum radiates real and virtual photons (brcmsstrahlung) . 

These photons, due to the dynamical features of QED, are mainly radiated along the 
direction of motion of the radiating particles, and their effect is mainly to reduce the 
original four-momentum of the incoming electron (positron) to xii^2)P~(+)- After this 
"pre-emission" , the "hard" subprocess e~ {x\p-)e'^ (x2P+) — ^ X takes place, at a re- 
duced centre of mass energy squared s = x\X2S. The resulting cross section, corrected 
for IS QED radiation, can be represented as follows [343]: 



(172) 




dxidx2D{xi, s)D{x2, s)ao{x\X2s)Q{c\its), 



where D{x, s) are the electron structure functions, representing the probability that an 
incoming electron (positron) radiates a coUinear photon, retaining a fraction x of its 
original momentum, and 0(cuts) represent the rejection algorithm implemented in order 
to take care of experimental cuts (see Fig. 62 for a graphical representation). 



fi ^2 f„ 




Fig. 62. - Graphical representation of eq. (172). 

Equation (172) takes into account only IS photonic radiation, for the sake of simplic- 
ity. It has however to be noticed that in the literature more acciiratc LL results can be 
found, taking into account also kinematical effects due to IS radiation and/or the effect 
of final-state radiation for the small-angle Bhabha process [21, 122]. 

Actually, when considering photonic radiation only the non-singlet part of the SF 
is of interest. If the running of the QED coupling constant a{s) is neglected (which 
corresponds to neglecting the corrections due to additional pair production, that can be 
taken into accoimt properly by moans of dedicated formulae [121]), the non-singlet part 
of the SF is the solution of the following Renormalization Group (RG) equation: 



(173) 
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where P+{z) is the so called Altarelli-Parisi (AP) splitting function [344], given by 

P+{z) = P{z)-5{l-z) [ dxP{x), 

Jo 

(174) P{z) = 

1 — z 

Equation (173) can be also transformed into an integral equation, subject to the boundary 
condition D{x, ml) = 5{1 — x): 

(175) D(..s) ^ ^(1 - «) + ^ f /; f p.(.)i. (f ^ 

Equation (175) can be solved exactly by making use of numerical methods. One of 
these is the so called inverse Mellin transform method, which is briefly described in the 
following. The Mellin transformation is defined as 

(176) F{n) = [ da;a;"-V(a;), 



r+'°° dn 

' -y—ioo 



(177) fix) = j^^ ^^""^W- 



By taking the Mellin transform of eq. (172), switching to the differential form of the 
evolution equation and solving for the Mellin moments, one finds 

(178) D(n, s) = exp |c(n) , 

where C{n) are the Mellin moments of the AP splitting function 

(179) C(n) = / dzz^'-^P+iz) 

Jo 

and 

(180) r? = — In ' * 



By computing the moments of eq. (179), inserting them into (178), and performing nu- 
merically the anti-transformation given by (177), the exact numerical solution of (175) 
is obtained. There is only a technical remark at this point: the integrand of the anti- 
transformation is strongly oscillating in the asymptotic region, so that a proper regular- 
ization procedure must be adopted in order to obtain numerical convergence. 

For practical purposes, it is also of great interest to obtain approximate analytical 
representations of the solution of the evolution equation. A first analytical solution can 
be obtained in the so called soft photon approximation, i.e. in the limit a; ~ 1. In such a 
limit, the dominant contribution to eq. (177) comes from the large n region, so that one 
can approximate the moments of the AP splitting function by means of their asymptotic 
expansion, namely 

(181) C(n) ~^-27b-2 Inn -F0(- 



PRECISION PHYSICS AT LEP 



131 



7b being the Euler constant. By inserting eq. (181) into eq. (178) and computing the 
anti-transformation one obtains 



(182) 



Dgl{x,s) = 



_ exp [^?7 (I - je)] 1 



r(i 



2^(1 



X) 



r being the Euler gamma-function. The solution shown in eq. (182) is known as the 
Gribov-Lipatov (GL) approximation [345]. Its main feature is that it exponentiates at 
all perturbative orders the large logarithmic contributions of the kind 77ln(l — x). Its 
main drawback is that it is valid only in the soft limit, i.e. it does not take into account 
properly hard-photon effects. 

The evolution equation (175) can also be solved iteratively. At the n-th step of the 
iteration, one obtains the 0{a") contribution to the structure function. The iterative 
solution of the evolution equation is, in some sense, complementary to the exponenti- 
ated GL solution. Actually, on the one hand the GL solution is exponentiated at all 
perturbative orders, whereas the iterative solution must be truncated at a given finite 
perturbative order. On the other hand, the GL solution is valid in the limit a; ~ 1, 
whereas any given iterative contribution can be computed exactly. 

In order to go beyond the soft-photon approximation, the following general strategy 
has been adopted. Given an iterative solution up to a given perturbative order, for each 
perturbative contribution it is possible to isolate the part which is contained in the GL 
solution. Then, by combining the GL solution with the iterative one, in which that part 
has been eliminated in order to avoid double counting, one can build a hybrid solution 
of the evolution equation, which exploits all the positive features of the two kinds of 
solutions and docs not present anymore the limitations intrinsic to each of them. As 
a matter of fact, in the literature it is possible to find two classes of hybrid solutions, 
namely the additive and factorized ones. The additive solutions have the following general 
structure: they are built by simply adding to the GL solution the finite order terms 
computed by means of the iterative solutions. The factorized solutions, on the contrary, 
are obtained by multiplying the GL solution by finite order terms, in such a way that, 
order by order, the iterative contributions are exactly recovered. In the following, the 
third order exponentiated additive and factorized solutions are reported. The additive 
one reads [239]: 



Da{x,s 

S^\x,s 

d^2\x,s 
d^-* {x,s 



= ^rf^^(a;,s), 

i=0 

= Dgl{x,s), 
= -7V{^ + x), 



1 2 

= 32^ 



(1 + x) (-41n(l -x)+3lnx)-A 



In: 



1 - X 



b - X 



= {{i + x) [18C(2) - 6U2{x) - 12 In'il - x)] 



+ 



1-x 



--(1 + 8x + 3a;2) Ina; - 6{x + 5)(1 - x) ln(l - x) 



-12(1 + x^) Ina; ln(l -x)^-{\^ Ix^) In^ x 
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--(39- 24a;- 15a;2) 
4 



where C is the Riemann ^-function. The factorized one reads [119,238]: 



Df{x,s) = DGL{x,s)^d^p' , 

i=0 



(0) 



^ 42 



--(l + 3x2)lna;- (l-a;)2 



j(2) 



8 V2/ 



(l-a;)2 + -(3a;2-4a; + l)lna; 



(184) 



+ — (1 + 7a;2) In^ a; + (1 - a;2)Li2(l - a;) 



It is worth noting that also higher order solutions are known, numerically and/or ana- 
lytically, for both the classes of solution. However, LEP phenomenology is not sensitive 
to hard photonic terms beyond the third order. 

If only cuts on the invariant mass of the event after ISR are considered, eq. (172) 
can be simplified by performing explicitly one of the integrations. Actually, for s' > xqs, 
it can be shown that the QED corrected cross section can be written as 



(185) 



'^{^) — / dxH{x, s)ao {{1 — x)s) , 
Jo 



where H{x, s) is the so called radiator, or flux function, defined as 



dz 



D{z,s)D 



l-x 



\-X 



(186) H{x, s) 
By defining 

(187) K{x, s) = H{1 - X, s) 

and taking the Mellin moments of eq. (187), the following identity can be shown: 



K(n, s) = / dxx''-^K{x, s) = D^{n, s) = exp \^C{n) 



(188) 



where C(n) are the Mellin moments of the AP splitting function, given by eq. (181). In 
virtue of eq. (188), the following identity holds: 



(189) 



Hr,{x,s) = £>2^(1 - X,S). 



It is worth noticing that in QED there are some non-leading corrections that behave 
much like the leading logarithmic ones. Actually, it can be verified by inspection that 
after the substitution 



(190) 



2a 

TT 



In 



mi 



the 0{a) distribution of the invariant mass of the event after ISR as extracted from 
eq. (185) reproduces exactly the standard diagrammatic calculation. 
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Moreover, it has to be noticed that it is possible to match the all-orders LL results of 
the SF method with exact finite-order diagrammatic results. For the line shape, the 
match has been performed in the soft-photon approximation up to O(a^), according to 
the procedure described in Sect. 3'4. Also a more general procedure for implementing the 
matching beyond the soft-photon approximation, and hence taking into account hard- 
photon effects, is known and can be found in ref. [23] . 

A' 2. The Parton Shower Method. - The Parton Shower method is substantially a 
method for providing a Monte Carlo iterative solution of the evolution equation, at the 
same time generating the four-momenta of the electron and photon at a given step of 
the iteration. It has been studied within the context of QCD (see for instance rcfs. [346] 
and [347]) and subsequently developed also for QED (see refs. [348], [349], [350] and 
references therein). In order to implement the algorithm, it is first necessary to assume 
the existence of an upper limit for x, in such a way that the AP splitting is regularized 
as follows: ^ 

(191) P+{z) = e{x+ - z)P{z) -6{l-z) [ ^ dxP{x). 

Jo 

Of course, in the limit x+ ^ 1 cq. (191) recovers the usual definition of the AP splitting 
function given in eq. (174). By inserting the modified AP vertex into eq. (173), one 
obtains 

(192) s^^D{x, s) = ^j ^P{z)D (-, s) - ^D{x, s) j dzP{z). 
Separating now the variables and introducing the Sudakov form factor 



(193) n(si,s2)=exp 



which is the probability that the electron evolves from virtuality — S2 to — Si without 
emitting photons of energy fraction larger than 1 — a;+, eq. (192) can be recast into 
integral form as follows: 

(194) D{x,s) = Jl{s,ml)D{x,ml) + ^ £^ ^n(s,s') ^P{z)D . 

The formal iterative solution of eq. (194) can be represented by the following infinite 
series: 



(195) n(s„, mi)D , to; 

\Zl - --Zn 

The particular form of (195) allows to use a Monte Carlo method for building the solution 
iteratively. The algorithm is standard, and can be described as follows: 

1 - set = TOg, and fix a; = 1 according to the boundary condition D{x,ml) = 
(5(1 -x); 



2 - generate a random number ^ in the interval [0, 1]; 
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3 - if ^ < n(s, Q^) stop the evolution; otherwise 

4 - compute Q'^ as a solution of the equation ^ = Il{Q'^, Q^); 

5 - generate a random number z according to the probability density P{z) in the 

interval [0, x+J; 

6 - substitute x xz and Q'^; go to 2. 

An important feature of this algorithm is that it can be used to generate exclusive events 
containing the complete information on the four-momenta of the particles. This fact is 
well known for the QCD analog of the algorithm, and relics upon the kincmatical rules 
of the so called jet calculus. For the case under consideration, in which at the i-th step of 
the iteration a virtual parent electron branches into a virtual electron plus a real photon, 
the space- like kinematics apply. In particular, one has 

(196) eip)^e{p')+^{q), 

where = —K^, p'^ = —K'"^ and = 0. Assuming that the parent electron moves 
along the z axis, the four-momenta of the particles can be parameterized as follows: 

p={E,0,p,), 
p' = {E',p±,zp,), 

(197) q={E^,-p^,{l-z)p,). 

From the kinematics of eq. (197), in the limit Pz — > oo, the following relations can be 
derived: 

E^Pz-—, 

^Pz 

E zpz ' 



2zpz 2zpz 

(198) ^;^.(i_,),, + _^^ 
from which, by imposing energy conservation, one obtains 

(199) -K^ = -^+ 



z z{l — z)' 

Given K^, K'^ and z, eq. (199) allows one to compute the absolute value of the transverse 
momentum at the branching. 

A'3. YFS Exponentiation. - The Yennie-Prautschi-Suura exponentiation procedure 

is essentially a technique for summing up all the infra-red (IR) singularities present in 
any given process accompanied by photonic radiation [351]. It is inherently exclusive, 
i.e. all the summations of the IR-singular contributions are done before any phase-space 
integration over the virtual- or real-photon four- momenta arc performed [352]. In the 
following, the general idea underlying the procedure will be shortly described (see for 
instance [353] for a more detailed analysis). 

Let us consider the process e+(pi)e"(p2) ^ fiili) ■ ■ ■ .f«(g«), where /i(gi) • • • /„(?«) 
represents a given arbitrary final state, and let ^Ao be its matrix element. By using 
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Fig. 63. - The Feynman diagrams for the initial-state 0(a) bremsstrahlung correction. 



standard Feynman-diagram techniques, it is possible to show that the same process, 
when accompanied by / additional real photons radiated by the initial-state particles, 
and under the assumption that the I additional photons are soft, i.e. their energy is 
much smaller that any energy scale involved in the process, can be described by the 
factorized matrix element built up by the lowest order one, Ado, times the product of I 
eikonal currents, namely 



' r 



(200) 



ki-p2 



ki-pi 



where e is the electron charge, fc, are the momenta of the photons and £i{ki) their po- 
larization vectors (see Fig. 63 for a representation of the Feynman diagrams relative to 
the initial-state 0{a) bremsstrahlung correction). By taking the square of the matrix 
element in eq. (200), multiplying for the proper flux factor and the Lorentz-invariant 
phase space volume element, neglecting the four-momenta of the radiated photons in 
the overall energy-momentum conservation, summing over the final state photons po- 
larizations and combining properly all the factors, the cross section for the process 
e+(j3i)e~(p2) hill) • • • fn{Qn) + l real photons can be written as 



£i{ki)-P2 £i{ki)-Pi 



(201)da« =dao^n hdhd cos -dd^.-^Y^e^ f^JMl^ 
l\ 2{2-Kf'^ \ ki-p2 

1 — i L 



ki-pi 



By summing now on the number of final-state photons, one obtains the cross section for 
the original process accompanied by an arbitrary number of real photons, namely 



(202) 



1=0 

= dao exp 



kdkdcosM^^ye^ ( 4^ - 

^2(27r)3^ \ k-p2 k-pi 



Equation (202), being limited to real radiation only, is IR divergent once the phase space 

integrations are performed down to zero photonic energy. This problem, as well known, 
finds its solution in the matching between real and virtual photonic radiation. At any 
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Fig. 64. - The Feynman diagram for the photonic vertex-hke correction. 



rate, cq. (202) already shows the key feature of exclusive exponentiation, i.e. summing 
up all the perturbative contributions before performing any phase space integration. 

In order to get meaningful radiative corrections, besides initial-state real photon 
corrections it is necessary to consider also initial-state virtual photon corrections, i.e. the 
corrections due to additional internal photon lines connecting the initial-state electron 
and positron (see Fig. 64). For a vertex- type amplitude, the result can be written as 



' (27r)4 J + ie'^^P^'^ 2pi • + fc2 + is 



where T stands for the Dirac structure competing to the lowest order process, in such a 
way that Mq = v{pi)Tu{p2)- The soft-photon part of the amplitude can be extracted by 
taking A;'' ~ in all the numerators. In this approximation, the amplitude of eq. (203) 
becomes 

Mvi = Mo X V 
(204) v=^JdH> ^ ^P'-P' 



(27r)3 J h?- +ie {2pi ■ k + k'^ + ie){2p2 ■ k + k'^ + ie) 

Some comments are in order here. First, as in the real case, the IR virtual correction 
factorizes off the lowest order matrix element, so that it is universal, i.e. independent of 
the details of the lowest order process under consideration. Moreover, it is also free of 
ultraviolet (UV) infinities; this, of course, does not mean that the UV behaviour of the 
complete amplitude is irrelevant, but that the IR part can be treated independently of 
renormalization problems. Last, as in the real case, it is divergent in the IR portion of 
the phase space. 

The correction given by n soft virtual photons can be seen to factorize with an 
additional 1/n! factor [351], namely 



(205) 



Mv^=Mo X ^F", 
n! 
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so that by summing over all the additional soft virtual photons one obtains 
(206) Mv=Moxe^[V]. 

As already noticed, both the real and virtual factors are IR divergent. In order to 

obtain meaningful expressions, one has to adopt some rcgularization procedure. One 
possible regularization procedure is to give the photon a (small) mass A and modifying 
eqs. (201) and (204) accordingly. Once all the expressions are properly regularized, 
one can write down the YFS master formula, which takes into account real and virtual 
photonic corrections to the lowest order process. In virtue of the factorization properties 
discussed above, the master formula can be obtained from eq. (202) with the substitution 
dao dao \ exp(F)|^, i.e. 



{207)da = dao\ exp(l/)p exp 



kdkd cos 'dd(fi 



2(27r)= 



e{k)-p2 e{k)-pi 



k-p2 



k-pi 



As a last step, it is possible to perform analytically the IR cancellation between virtual 
and very soft real photons. Actually, since very soft real photons do not affect the 
kinematics of the process, the real photon exponent can be split into a contribution 
coming from photons with energy less than a cutoff kmin plus a contribution coming from 
photons with energy above the same cutoff. The first contribution can be integrated over 
all its phase space and then combined with the virtual exponent. The physical meaning 
of this procedure is to sum over all the degenerate final states: in fact, a very soft real 
photon produces a signature that is indistinguishable from the signature typical of the 
elastic event. After this step it is possible to remove the regularizing photon mass by 
taking the limit A ^ 0, so that eq. (207) becomes 



da = dao exp(y) exp 
(208) 

where Y is given by 



kdkdQ{k - kmin) cos ^'^'/^ 2(27r)3 E ' 



2 f sik) ■P2 e{k) -pi 



(209)y = 2V+ [ kdkde{kmin - k) cosMip } V ( 
J 2(27r) 



k-p2 



k-pi 



e{k)-p2 e{k)-pi 



k-p2 



k-pi 



The explicit form of Y can be derived by performing all the details of the calculation, 
and reads 



y = /31n- 



E 



+ <^YFS, 



1 ^ a f n 



(210) 



/3 = 



2a 



In I 



2 _ 1 

2 

1 



As in the SF method, the method of YFS exclusive exponentiation has been refined 
in order to take into account all-orders corrections on top of finite-order exact results 
(see for instance ref. [28] and references therein). 



138 



G. MONTAGNA, O. NICROSINI and F. PICCININI 




k 



V 



k-q 



Fig. 65. - Tlie fermionic vacuum polarization Feynman graph. 



B. — Vacuum Polarization Corrections 

The input parameter adopted for the analysis of the precision electroweak data 
collected at LEP are the QED coupling constant at zero momentum transfer a, the 
Fermi constant and the .Z-boson mass Mz- Given these input parameters, all the 
other unknown (or poorly known) parameters, such as the top-quark mass, the Higgs- 
boson mass and the QCD coupling constant as, are determined by means of best fits. 

Actually, the QED coupling constant a is known with great precision at zero mo- 
mentum transfer from measurements such as the electron/muon g — 2, the Lamb shift, 
the muonium hyperfine splitting, the neutron compton wavelength, the quantum Hall 
and the Josephson effects [354,355]; but for precision physics at LEP it must be evolved 
up to the Z-boson mass scale, so that the relevant parameter is rather the running QED 
coupling constant a{s). As is well known, the running of a is largely dominated by the 
contribution of fermionic vacuum polarization diagrams, which represent by themselves a 
gauge-invariant and universal subset of radiative corrections for two-formion production 
processes (see Fig. 65). The bosonic contribution to vacuum polarization is, on the con- 
trary, gauge-dependent. Its contribution is small, at least around the resonance, and 
can be taken into ac;count togetlica' with other gauge- and process-dependent radiative 
corrections at a given perturbative order. From now on, only the fermionic contribution 
to vacuum polarization will be discussed. 

The QED running coupling constant is given by 



The vacuum polarization contribution can be naturally split into the contributions 
coming from the leptons, the light hadrons and the top quark according to the following 
relation: 




(213) 



Aa{s) = Aai{s) + Aa]^> (s) + Aat{s). 




Fig. 66. - The integration contour C to be used for the integral of eq. (217) in order to obtain 
the dispersion relation of eq. (218). 



The contributions from leptons and top-quark loops can be safely computed by 
means of ordinary perturbation theory. In particular, the leptonic contributions can be 
very well approximated by the expression of the loop- integrals for « s, and read 

On the other hand, the top-quark contribution, due to the decoupling properties of QED, 

can be represented by 

In particular, eqs. (214) and (215), when evaluated at s = Mf give the values 

Aa;(M|) = 0.03142, 
(216) Aat(M|) = -0.6 x lO"'^. 



For the %/i^-quark contributions, on the contrary, the perturbative expression of 
eq. (214) cannot be used. In fact, due to the ambiguities inherent in the definition 
of the %/it-quark masses, the answer provided by eq. (214) is affected by very large 
uncertainties, proportional to Sniq/mq. Moreover, one can a priori expect very large 
non-pcrturbativc QCD corrections in the region of the hadronic resonances. The stan- 
dard procedure adopted in order to circumvent this problem is to exploit the analyticity 
properties of the vacuum polarization amplitudes by making use of dispersion relations 
(DR) techniques (sec for instance [357]). 

Given a complex- valued function F of complex argument s, if is olomorphic in a 
region TZ and on its boundary C, the value of F at any point inside TZ can be computed, 
according to the Cauchy's integral representation, as a contour integral: 



Under the assumption that 




Fig. 67. - The total cross section for the process e e — > hadrons, normalized to the total cross 
section for the process e+e~ ^ M^A*^ ^^w centre of mass energy (from [356]). 



• F{s) is real for real s, up to a threshold AP, 

• F{s) has a branch cut for real s > M^, 

• F(s) is olomorphic except along the branch cut, 

and taking as integration contour the one shown in Fig. 66, one can derive the following 
once-subtracted DR 

g2-q2 /-oo luiF{s) 



(218) F(g2)^F(g2) + ^_*I / 



A/2 s~(gs~q^ - is 

Equation (218) allows one to compute the value of F at any point with the only 
knowledge of the imaginary part of F along the branch cut. The optical theorem provides 
the link between the absorptive part of the hadronic vacuum polarization amplitude and 
the total cross section for the reaction e^e~ — > hadrons, namely 

(219) a{s) = ^h i^n^[s). 
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luminosity measurement 



Burkhardt, Pietrzyk '95 



contribution 
in magnitude 



in uncertainty 




'.05. 



narrow resonances 



1.05- 2.5 GeV 




narrow resonances 



Fig. 68. - The relative contributions to Aa{t = —1.424 GeV^) in magnitude and uncertainty 
(from [356]). 



Burkhardt, Pietrzyk '95 



> 12.GeV 





5 - 7 GeV 



1.05 -2.5 GeV 



7 -12 GeV 



2.5 - 5 GeV 



Fig. 69. - The relative contributions to Aq(s = M|) in magnitude and uncertainty (from [356]). 



Hence, from the experimental data for e+e — s- hadrons (see Fig. 67) one can compute 
the hadronic contribution of the hadronic vacuum polarization as follows: 



(220) 



47r2( 



-Re 



ds'- 



4mi 



s' — s — le 



Several determinations of the hadronic contribution to the running of the QED 
couphng constant a have been performed [356,358-367]. The differences between the 
various determinations can be traced back to the procedures adopted for fitting the 
data, treating the experimental errors and performing the numerical integrations. The 
differences are also determined by the different thresholds chosen to start the application 
of perturbative QCD at large s, and to the value used for as{Mz)- Figures 68 and 69 
show the relative contributions to Aa in magnitude and uncertainty at a scale typical of 
the luminosity measurement (Fig. 68) and at the peak (Fig. 69). In Fig. 70 several 
determinations of Aa(Mz) and a~'^{Mz) are shown, together with the value used by 
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>,.cqcrlehner 1 986 
Lynn et al. 1987 
Burkhardt et al. 1989 
...eqer lehner 1 99 1 
Nevzorov et al. 1 994 
Adel et al. 1995 
Burkhardt et al. 1995 
Eidelman et al. 1 995 
Martin et al. 1 995 
Swartz 1996 
Alemany et al. 1 997 



0.026 0.028 0.03 128.7 128.8 128.9 '29 

Aa,„, (M,) (M,) 

Fig. 70. - The determinations of A.a'^\Mz) and a~^{Mz). The dashed vertical Une selects the 
value used by the LEP Collaborations. 

the LEP Collaborations. (*^) 

For practical purposes, the running of the QED coupling constant for large time- 
like momenta can be computed by means of the following effective formula, in terms of 
fermion masses and effective light-quavk masses: 

(221) ^„(,).£|)^05«^,(,„_|_|j. 

In cq. (221), the sum is extended to all light fermions (chargc;d leptons and light quarks), 
Qf is the fermion charge in units of the electron charge and Ncf is 1 for leptons and 3 
for quarks. The effective light-qaaik masses = m„ = 47 MeV, rus = 150 MeV, nic = 
1.55 GeV and m;, = 4.7 GeV insure that the hadronic vacuum polarization contribution 
of refs. [356, 364] is reproduced. 

C. — One- loop Integrals and Dimensional Regularization 

The evaluation of one-loop diagrams generally leads to the problem of ultraviolet di- 
vergences, so that a regularization procedure has to be adopted in order to deal with well 



Two new determinations of a{Mz) have recently appeared [368,369] giving a ^(Mz) = 128.923 ± 
0.036 and a^^(Mz) = 128.98 it 0.06, respectively. After the completion of this work, a further analysis 
of a{Mz) has been performed in [370], yelding the result a-^{Mz) = 128.928 ± 0.023. 
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defined objects before implementing a renormalization program. As said in Sect. 3'2.1 
the method commonly used in the case of gauge theories is the dimensional regulariza- 
tion [45 47] , where the singularities arise as poles of the form 1/(4 — n) , n being the 
number of space time dimensions where the integrals are convergent. In the literature 
a general method has been developed [48], which allows to write every loop diagram in 
terms of certain combinations of scalar form factors [49]. In this appendix some technical 
details concerning this procedure are provided for the simple case of the two-point form 
factors. In the evaluation of a self-energy diagram the following integrals are encountered: 

(222) Bo;B^;B^^{p^,mi,m2) = A / l^i^rl^\2 ^ 2v 

iTT^ J [q^ + ml){{q + pY + m^) 

The vector and tensor integrals can be decomposed into Lorentz covariants and scalar 

coefficients: 

(223) B^{p'^,mi,m2) = p^Bi(p^.mi, TO2), 



(224) Bij,^{p'^,mi,m2) = PnPvB2i{p'^ ,mi,m2) + 5^,,B22{p'^ ■,mi,m2). 

The functions Si, B21 and ,622 are algebraically related to the two-point scalar integral 
Bq and to the one-point scalar integral A{m), defined by 

(225) A{m) = ^ / 

The expression for Bi can be obtained multiplying eq. (223) by Pij, and taking into 
account of the identity 

(226) q-p=\[{q + pf +ml- {q^ + ml) - / - + mf] . 
It follows that 

(227) ^^^^ ~ ^^""2) - (P' +ml- mi)Bo] , 

where the arguments have been omitted for simplicity of notation. The expressions for 
B21 and B22 are analogously obtained multiplying eq. (224) by P,mPv and S/^v- 

(228) B22 = ^ l-mfBo + ^[^(ms) - (m? + ml + i/) - (m? - -p^)^,] 



(229) B21 = , „ 



1 

372 



A{m2) + 2(m? - - p^)Bi + hml+ml + ^p^) + m\Bo 



In so doing only the basic scalar integrals A{m) and Bq{p^ , mi, 7722) need to be calculated. 
Let us consider the one-point integral A{m) in four dimensions: 

(230) A{m) = / dqo j d^'q '^^ . . 

Jo J q^ + m-^ - IS 

In the complex qq plane the integrand has poles for 

(231) go = ± [VWT^ - ze) . 
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The introduction of the Wick rotation allows to transform the Minkowski space in a 
Euclidean space and perform the integration by means of polar coordinates. 

Working in Euclidean space in the case of n dimensions, a set of generalized polar 
coordinates can be introduced [47, 371] 



''^'&n-ld'&n-l- 



(Pq= uj"-^dLj / A9i / sini?2«it?2... / sin" 
Jo Jo Jq Jo 

Since the integrand of A(m) depends only on q^, the angular integrations are easily 
performed by means of the relation 

(233) / 

Jo 

yielding 



Y ( ™+i ] 



(234) A{m) = / u:^-'du;- 







(235) 



2^t-2 1 r(f)r(i 



r(f) 2 r(i)m2(i-t) 



where r(x) is the Euler T function with the property a;r(a;) = r(a; + 1). Since the 
physical results are recovered for n approaching 4, it is convenient to introduce the 
notation n = 4 — e, and expand the expression of A{m) around e = 0: 



(237) A{m) = m\~^m 



r(|) 



-1 + ^ 



(238) ~ (-A + lnm^ - 1) , 
where 

(239) A=— 7-ln7r 

4 — n 

contains the ultraviolet divergence, and the expansion of the F function around the zero 
of its argument has been used 

(240) T{x) = --^ + 0{x), 

X 

with 7 indicating the Eulcr constant. 

The calculation of the two-point scalar integral Bq is more involved than A(rn) 
because of the presence of two factors in the denominator. As a first step it is necessary 
to use the Feynman parameterization to combine the denominators: 

1 /"^ 1 

(241) —= dx 



ah Jq [ax + 6(1 — x)]"^ ' 
where 

b = + ml, 

a = {q+pf + ^1- 
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With a shift in the integration variable q q + xp the term hnear in the integration 
variable disappears, yielding 



(242) Bo(/,mi,m2) = / dx [ d""^ 



(?^ + X - ^e)^ ' 



with X = —p^x^ + [p^ + m\ — m\)x + m\. The angular integration is now straightforward 
and the last integral is again worked out in terms of the F function: 

(243) Bo{p\mi,m2) = ^^^T,r^. 

X r(2) 

By means of an expansion of n around 4 as for the case of the one-point scalar integral, 
the expression for Bq reads: 

(244) Bo{p^,mi,m2)=A- [ dxlnx- 

Jo 

For arbitrary values of momentum and masses the integral in cq. (244) can be worked 
out by means of the method outlined in ref. [50] . The integration becomes straightforward 
for particular values of the arguments. For example, to obtain the asymptotic expressions 
of the fermionic self-energies written in Sect. 3'2.1, the following asymptotic expressions 
for the Bo scalar function are useful: 

So(p^,m,m)|p2|^„2 ~ A-ln(-p2) + 2 + iTT, 

2 4 

Bo(p^m,m)|p2|«„2 ~ A-lnm^ - ^ + 

Bo(p^ 0, m)|p2|»„2 ~ A - ln(-p2) + 2 + 

So(p^,0,m)|p2|«„2 ~ A - Inm^ + 1 - 

Concerning the vertex correction diagrams, the following three-point integrals are 
encountered: 



(1)(2)(3)' 



(245) Co;Ci^;Cfj,^{pi,p2,mi,m2,ms) = ^ ( d'^-q 

ITT^ J 

where 

{l)=q^+ml 

(2) = {q + p^f + ml 

(3) = {q+Pi+P2f +ml 

The vector and tensor integrals can be decomposed into Lorentz covariants and 
scalar coefficients in the following way: 

C^,=p'^Cn+P^C,2, 

where the arguments in the functions Cij have been neglected for simplicity of notations. 
According to the procedure illustrated in ref. [48], the three-point scalar form factors 
can be expressed as linear combinations of two-point functions and of the fundamental 
three-point scalar integral Co. It is worth noticing that only the function C2A contains 
the UV divergences, while the other form factors are UV finite. 
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